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AN EXPERIMENTAL ANALYSIS OF CRITICAL FACTORS 
INVOLVED IN THE BREAKDOWN PROCESS 
OF LEADING EDGE VORTEX FLOWS 


Abstract 

by 

Kenneth D. Visser 


Experimental crosswire measurements of the flowfield above a 70° and 75° 
flat plate delta wing were performed at a Reynolds number of 250,000. Survey grids 
were taken normal to the planform at a series of chordwise locations for angles of 
attack of 20° and 30°. Axial and azimuthal vorticity distributions were derived from 
the velocity fields. The dependence of circulation on distance from the vortex core as 
well as on chordwise location was examined. The effects of nondimensionalization 
in comparison with other experimental data was made. 

The circulation distribution scales with the local semispan and grows 
approximately linearly in the chordwise direction. For regions of the flow outside of 
the vortex subcore, the circulation at any chordwise station was observed to vary 
logarithmically with distance from the vortex axis. The circulation was also found to 
increase linearly with angle of incidence at a given chordwise station. 

A reduction in the local circulation about the vortex axis occurred at 
breakdown. The spanwise distribution of axial vorticity was severely altered through 
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the breakdown region and the spanwise distribution of axial vorticity present 
appeared to reach a maximum immediately preceding breakdown. The local 
concentration of axial vorticity about the vortex axis was reduced while the magnitude 
of the azimuthal vorticity decreased throughout the breakdown zone. The axial 
vorticity components with a negative sense, found in the secondary vortex, remained 
unaffected by changes in wing sweep or angle of attack, in direct contrast to the 
positive components. The inclusion of the local wing geometry into a previously 
derived correlation parameter indicated the circulation of growing leading edge vortex 
flows to be similar at corresponding radii from the vortex axis. 

It was concluded that the flow over a delta wing, upstream of the breakdown 
regions and away from the apex and trailing edge regions, is conical. In addition, the 
dominating factors leading to the onset of breakdown are felt to be the local 
circulation of the vortex and the accompanying pressure field. 
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CHAPTER I 


INTRODUCTION AND BACKGROUND 


1.1 Synopsis 

A vortex structure is a truly fascinating physical phenomenon. The ability of a rotating 
fluid to maintain a cohesive structure that would surely seem to tear itself apart at a moments notice 
is of itself worth noting. Yet, vortices are present in many situations, from giant funnel clouds to 
the stirring of coffee in a cup. They can persist in a regular fashion in the wake of a flat plate or 
may form and disperse quickly behind blunt bodies. Common to all of these structures, however, 
is that their behavior can be generalized to that of a mass of fluid rotating about a common center. 

Certain types of vortex structures will exhibit a transition in their state from a coherent, 
relatively high rotational speed flow to a much slower form. Almost instantly, a rapid change will 
occur in the flowfield including an expansion of the central core region and an associated rise in 
pressure along the vortex axis, which drastically alters the downstream characteristics. Leading 
edge vortex flows, such as that generated by delta wing planforms, also experience a severe 
decrease in the core's axial and circumferential velocity and large scale turbulent dissipation as a 
direct result of this 'vortex breakdown' or 'bursting'. The underlying physics of this apparently 
stable flowfield combined with an seemingly unpredictable disruptive process has interested 
aerodynamic researchers for over 35 years. 

More specifically, the process by which breakdown occurs above a delta wing planform is 
of great concern from a flight performance point of view. The nature of a typical delta wing vortex 
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is illustrated in Figure 1.1. At high angles of attack, the boundary layer on the lower surface of the 
delta wing flows outboard. The fluid separates at the sharp leading edge, forming a free shear 
layer which curves upward and rolls into a core of high vorticity above the top side of the wing. 
Each of the two counter-rotating vortices also contain axial flow components in the central core 
regions, around and along which the fluid spirals. 



Figure 1.1 Delta Wing Vortex Flow Schematic 

Additional spanwise outflow is induced on the upper surface beneath the coiled vortex sheet. This 
spanwise flow over the surface separates as it approaches the leading edge to form the secondary 
vortex structures which have an opposite sense of rotation. The main effect of the secondary 
vortex is to displace the primary vortex upwards and inwards. The size and strength of the 
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primary vortex increases with angle of incidence. It becomes the dominant steady flow feature 
through a wide range of practical flight attitudes. The pressure field in these vortices results in an 
incremental lift termed the "vortex" or "nonlinear" lift 

Delta wing performance is limited by vortex breakdown. As the angle of incidence is 
increased, breakdown moves toward the apex. The large suction pressures created by the leading 
edge vortices are reduced, aft of the breakdown region. This increase in surface pressure results in 
a lower lift contribution by the vortex behind the breakdown region. Upon reaching the apex, a 
further increase in the angle of attack causes total separation of the leeward flow or stall. The effect 
of breakdown on the gross features of the flow can be seen in Figure 2a. A laser light sheet is 
used to visualize planes normal to the wing surface. The vortex on the far side of this delta wing 
has broken down, in contrast to the unburst vortex in the foreground. A photograph of the 'spiral' 
type of breakdown is presented in Figure 1.2b. Smoke used to trace the flow is placed into the 
central core of the vortex and vividly demonstrates the abruptness of this transition phenomena. 



Figure 1.2a Vortex Breakdown, Laser Sheet (Payne, 1987) 
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Figure 1.2b Vortex Breakdown, Spiral Breakdown 


As a further example of the significant effects of breakdown on the the vortex flowfield the 
total pressure distributions above a 70° sweep wing (Payne, 1987) at a chordwise station of x/c - 
0.5 are depicted in Figure 1 .3. At an angle of attack of 30 , the breakdown is aft of this station, 
while at ct = 40° the breakdown region has moved forward of this location. The post-breakdown 
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flow in Figure 1.3b has the scale of the pressure axis expanded to illustrate the local features. The 
pressure increase in the core, aft of breakdown, is substantial, and the expansion of the pre- 
breakdown, narrow core region is also evident. 

Early investigators observed the onset of breakdown in wind and water tunnels under 
various test conditions. Others immediately began theoretical and numerical investigations in an 
effort to not only understand what is occurring, but to predict the onset of this drastic change in the 
flowfield. Efforts have included variation of all manner of possible parameters, including wing 
geometry, angle of attack, and Reynolds numbers. Most investigators approach the study of delta 
wing vortical flows by asking such questions as: "Why does vortex breakdown occur?" and 
"What factors influence breakdown?" A different perspective on the problem can be gained by 
posing the inverse question: "Why does the vortex manage to maintain a cohesive flow structure 
upstream of the transition/breakdown region?", especially in light of rotational rates on the order of 
50,000 rpm. 

It is important to understand the physics of the vortex flowfield, both in general and that 
specific to delta wings. The following section describes the behavior of vortex flows, including an 
early model, and the characteristics of the specific case of the delta wing planform are overviewed. 
The remaining portion of the chapter is devoted to presentation of possible mechanisms for vortex 
breakdown. These are examined separately in light of the past theoretical, numerical and 
experimental explanations, which then leads to the present experimental study. 


1.2 Phenomenological Aspects of the Vortex Flowfield 


In order to gain a better understanding of the delta wing flowfield and its important 
characteristics, an overview of some basic concepts is given. This is followed by a review of an 



analytic model derived by Stewartson and Hall in the early 1960's which uses simplifying 
assumptions on the flowfield behavior as a means for obtaining a solution. Finally, additional 
concepts used to describe the vortex flow will be presented 
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1.2.1 Fundamental Vortex Behavior 

The vast majority of vortex structures resemble the classical two dimensional Rankine 
vortex as shown in Figure 1 .4. A strongly rotational central core region is surrounded by an 
irrotational region. The central region possesses a tangential velocity distribution which grows 
with increasing radius, such as V 0 = cor. This is termed a rotational or 'forced' vortex. At a 

specified radius, the velocity begins to decrease with increasing distance from the vortex axis in 
the form of a V 0 ~ 1/r distribution and the irrotational or 'free' vortex flow continues to infinity. 



Figure 1.4 Two Dimensional Rankine Vortex Structure 
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The presence of viscosity, however, tends to merge these regions such that a region exists 
between the purely irrotational flow and that of the core region. This region possesses rotational 
properties, and is described by a profile that combines aspects of both. An exact solution to the 
Navier Stokes equations for the decay of a vortex under the action of viscosity is presented by 
Batchelor (1967) with a tangential profile of the form 


V 0 (r,t) = v 



-r^/4vt 
e ) 


( 1 . 1 ) 


This profile is presented in Figure 1.5 along with V 0 (O,O) which remains undefined at r = 0. 

For t > 0, the vortex takes the form of an external irrotational flow matched to an inner core 
through an annular rotational, viscous region. Since vt is greater than zero for a timescale, t, 
based on the ratio of the root chord of the model to a freestream velocity, the resulting profile 
can be scaled to resemble that measured in real fluids. 



Figure 1.5 Viscous Vortex Structure (Batchelor, 1967) 
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Actual vortex structures occurring in nature are seen to be more complex than their two 
dimensional counterparts. Tornadoes possess large updraft velocities along their axis. The 
addition of some rotation to a fluid draining through a hole at the bottom of a container causes it 
to drain much more quickly leading to the conclusion that the interaction of the tangential and 
axial velocity component in a vortex is very significant 

The primary vortex structure above a delta wing, as well as that found in swirling tube 
flows, exhibit an additional velocity component. An axial, jetting type of profile exists normal 
to the plane of the vortex, reaching a maximum along the core axis as illustrated in Figure 1.6. 
Axial velocities have been measured in the vortex cores of delta wings to over three times that of 
the freestream. The core region also exhibits large tangential velocity gradients near its center. 
For the sake of a reasonable analysis, however, simplifying assumptions are generally made as 
will be shown. 


Z 



Figure 1.6 Vortex Velocity Field with Axial Component 


9 


Delta wing vortices are not symmetric and the influence of the wing surface can not be 
simply neglected. Leibovich (1984) summarized several aspects regarded as common to delta 
wing flows. The flow separating from the leading edge of a delta wing is generally regarded as 
a sheet of vorticity that spirals inward as it progresses down the wing. A core region is formed 
near the center part of this spiral. Vorticity is continually fed into the vortex structure from the 
boundary layer on the underside of the wing at the leading edge. The core region continues to 
enlarge in a conical fashion in the chordwise direction, as does the entire field, which results in 
an almost linear increase in the circulation with downstream distance. The axial velocity 
increases as a direct result of the drop in centerline pressure with distance from the apex. 

The definition of the core region is usually taken to be the axisymmetric center of the vortex 
which behaves in a rotational fashion. Five hole probe measurements by Eamshaw (1961) 
above a unity aspect ratio delta wing at a = 14.9° led to his division of the vortex core into three 

regions. He estimated the core to be approximately 30% of the local semispan in diameter and 
essentially conical in nature. Along the axis of this vortex core is the Viscous subcore', which 
measured about 5% local semispan in diameter and contains high velocity and large pressure 
gradients. The subcore exhibits the appearance of an axisymmetric rotating flow, rather than 
conical, along with some solid body rotation properties. Outside of the core lies a region 
containing the trace of the vortex sheet. A schematic in Figure 1.7 depicts these defined flow 
regions. It will be shown later that the subcore contains the majority of the axial vorticity. 
Leibovich (1984) noted that the Reynolds number has no effect on the size of this vortex core 
region and only serves to decrease the size of the subcore region as the Reynolds number 
increases. The axis of the vortex is generally located on the order of a core diameter above the 
wing. The vertical core locations, as measured by Payne (1987) for a wing sweep of 75° at a 
= 10° to 40° are presented in Figure 1.8. The core position maintains an approximately constant 
height, relative to the local span, above the wing at all chord stations for a given angle of 
incidence. This indicates, at the very least, that the vortex trajectory reflects a conical behavior 



VERTICAL LOCATION (2/S) 
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Figure 1.7 Three Regions of a Delta Wing Vortex Structure 

in the flow. In addition, the angle between the vortex trajectory and the planform surface can be 
seen to increase in an approximately linear fashion with angle of attack at the lower values of a. 



CHOROWISE LOCATION CX/C) 


Figure 1.8 Vortex Core Locations (Payne 1987) 
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Erickson (1981) mentions a core position variation for laminar versus turbulent boundary 
layer flows on the wing surface. The primary core moves upward and inward for the laminar 
case, relative to the turbulent one, because of secondary vortex displacement effects. Higher 
suction pressure peaks are evident in the turbulent case, but the integrated pressures are the 
same. He discusses in detail the effect of Reynolds number on vortex flow development. At 
higher Reynolds numbers, based on the chord, c, the ratio of scales of the transport process of 
diffusion to convection is on the order of 

1 

Re 0 5 


which is also a measure of the boundary layer thickness. It appears that large scale vortex 
structures are determined primarily by convective transport mechanisms, implying an 
independence of Reynolds number. Vortex core axial velocities, however, have been shown to 
be Reynolds number dependent Erickson concludes that the majority of the phenomena 
observed in the delta wing flow field is dominated by potential flow effects associated with the 
external field, that is the external pressure gradient. 


1.2.2 An Analytic Delta Wing Vortex Model 

One of the earliest analytical explanations for the flow observed in the primary vortex of a 
delta wing was the incompressible solution proposed by Stewartson and Hall (1963) based on a 
simplified model introduced by Hall (1961). This model has been found to compare quite 
favorably with experimental data by several researchers. The flow was approximated to follow 
a steady, axisymmetric behavior which simplified the cylindrical Navier-Stokes equations to the 
form below utilizing the geometry of Figure 1 .9. 



V 
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Stewartson and Hall assume the diffusion of vorticity is confined to a geometrically slender 

av r av r 

core. The direct result of this simplification is the elimination of the terms V r -^r and V x -^-. 

The flow contains no vortex sheets, and thus must be continuous and rotational to allow 
vorticity convection. The outer core pressure and velocity fields are considered to be conical, 
that is, maintaining constant values along rays emanating from the apex. An inviscid outer core 
solution is determined to specify the boundary conditions for the viscous inner core solution. 
The resulting solution for the inviscid flow regime is 


and 

where 


The outer radius of the core is defined by r = ax. The boldface type in the equations above 
denote the values of the flow velocities at the outer radius of the core. As r— >0 this solution 
breaks down because the ignored viscous terms become important. An order of magnitude 
analysis was then performed by Stewartson and Hall, similar to that of a boundary layer 
approximation. For a sufficiently small value of v, the viscous core will be slender and changes 
in the radial direction will be much greater than those axially. The original NS equations thus 
become 

Vp 2 1 SP 

= pSF 


' >• oln (s) 


v r ra 
V x = -2 x 


(1.3a), (1.3b) 
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(1.6a) 
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av x 

+ v x^ 


1_8P 
rt <5x 



7S( rV r) + l( V x)=° 


(1.6b) 

(1.6c) 

(1.6d) 


The introduction of suitable independent variables and asymptotic expansions are used to 
arrive at a set of ordinary differential equations solved by the Runge Kutta method. 
Interestingly, it was found that for suitably small values of the kinematic viscosity, the axial 
velocity component was nearly constant across the viscous region. The solution is presented in 
Figure 1.10a from Stewartson and Hall (1963). 



Figure 1.10a Vortex Model, Analytical Results (Stewartson and Hall, 1963) 
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Laser Doppler Velocimetry (LDV) measurements by Pagan and Solignac (1986) on a 
vortex in the wake of a 75° sweep delta wing at a = 19.3° are compared with Hall and 
Stewartson (1963) in Figure 1.10b. Good agreement is demonstrated, although the 
experimental subcore region is seen to be larger. The profiles of Pagan and Solignac also 
indicate three distinct vortex regions 

i) an internal, viscous, solid body rotation motion 

ii) an external, irrotational area with a potential velocity law of the 1/r form 

iii) in intermediate "Euler" region where the flow is inviscid and rotational. 



Figure 1.10b Vortex Model, Experimental comparison (Pagan and Solignac, 1986) 

Verhaagen and van Ransbeeck (1990) have compared their data to Stewartson and Hall by 
matching the numerical inner core solution at each chordwise location to the edge of their local 
rotational core. This was based on pressure and velocity measurements above a 76° sweep half 
span model. Verhaagen and van Ransbeeck defined their rotational core to be the region inside 
of which the vorticity is continuously distributed and no shear layer can be detected as shown in 
Figure 1.1 la. They found the velocity and pressure distributions along the core axis to be 
strongly over estimated by the theory and attributed this largely to the fact that the subcore size 
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and axial velocity variations were much greater than that assumed in the theoretical derivation. 
Correlation with the theory improved with chordwise distance from the apex. The distribution 
of the circumferential velocity in the inviscid region matched quite well. The core total pressure, 
however, was found to be much lower in the experiment than what the theory predicted. 

Overall , the theory agreed reasonably well with the experimental data in the range of x/c = 
0.3 to 0.7. An example of the axial velocity is given in Figure 1.1 lb at a chordwise station of 
x/c = 0.5. The good agreement should not come as a surprise. Stewartson and Hall constructed 
a theory using an isolated axisymmetric vortex. In the wind tunnel, as the vortex moves 
rearward on the wing, and subsequendy away from the planform surface, the influence of the 
wall decreases. In addition, near the apex, the wing span is on the order of the thickness of the 
model, resembling more of a blunt body than a slender wing. The vortex core size is also on the 



a) 


b) 


Figure 1.11 Validation of Flow Conicality (Verhaagen and van Ransbeeck, 1990) 
a) Rotational Core b) Axial Velocity 
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order of the wing thickness in this region which would invalidate the assumptions of a slender 
wing approach. Hence, it would be expected that the experimental data would correspond more 
closely to the predicted theory away from the apex and the trailing edge regions. 

The model outlined above represents one type of approximation used to predict vortex 
flows, that being an axially symmetric core with a continuously distributed vorticity. The other 
type of flow model employs a sheet of vorticity, emanating from the leading edge and rolling 
into a spiral as it winds into the center of the vortex. Mangier and Weber (1967) present a 
conical analysis using this type of reasoning. Experimental observations have indicated that the 
vorticity is confined to this sheet during the initial phase of roll up. Viscous effects are ignored 
by Mangier and Weber and the vorticity is assumed to be concentrated on a thin sheet, the flow 
between the sheets being irrotational. They also stare that that the total head for the core region 
remains constant, which seems questionable in light of the total pressure drop at the vortex 
center. Yet the leading terms of their asymptotic solution indicate that their solution for a 
potential flow with sheets of vorticity is identical to the solution by Hall (1961) of an 
axisymmetric flow with distributed vorticity. 


1.2.3 Further Concepts on Leading Edge Flows 

The analysis of the leading edge vortex flow and its effect on a delta wing planform has 
also been undertaken from an aerodynamic force point of view, as an alternate approach to 
predicting the physics of the fluid interactions using the equations of fluid motion. 
Polhamus(197 1) devised a method, referred to as the leading edge suction analogy, to predict 
the lift on a delta wing by separating the normal force into potential lift and vortex lift 
components. This can be written as 
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C L = K p sinacos 2 a + K y sin 2 acosa = C Lp + C Lv (1.7) 

and is illustrated in Figure 1.12. The potential term is based on lifting surface theory where K p 

is the potential lift coefficient slope at zero angle of attack. The lift due to the vortex term results 
from the suction generated by the equivalent attached flow about the edge and K y is estimated 

from the potential flow leading edge suction calculations. 



Figure 1.12 Vortex Lift Definitions (Hemsch, 1990) 

No knowledge of the shape, strength or position of the vortex is required. Lift measurements 
by Wentz and Kohlman (1971) on 70° and 75° delta wings are in excellent agreement with this 
analogy up to breakdown as shown in Figure 1.13. Surface pressure measurements by Er-El 
and Yitzhak (1988) on a 60 and 75 degree sweep delta wings also indicate that the analogy 
provided good predictions of normal force, potential and vortical components and normal force 
loadings up to a = 20°. 
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Figure 1.13 Vortex Lift Distribution (Wentz and Kohlman,1971) 


Hemsch (1990) makes a relevant point that can easily be overlooked when examining the 

resulting lift on a delta wing. The vortex lift increment is often referred to as the nonlinear lift 
increment, C^, when in fact the nonlinear increment is defined as 

%-C L -Kpa (1.8) 

and illustrated in Figure 1.12. His analysis has shown that although the nonlinear lift 
component increases with increasing sweep, the vortex lift component decreases as does the 
vortex strength as noted in Figure 1.14. 
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Figure 1.14 Vortex Lift Coefficient Behavior (Hemsch, 1990) 


The correlation of experimental data about highly swept delta wings, using similarity 
parameters, is an additional means of gaining insight into the behavior of the leading edge flow. 
Sychev (1960) introduced parameters for hypersonic flows about slender bodies which were 
extended by Barnwell (1987) to the inner region of subsonic/transonic flow. The flow is 
implied to have smooth, small axial gradients at high angles of attack. Sychev analyzed the 
Euler equations using a geometrical slenderness assumption 

g maximum semispan 

~ body length or wing chord <<: 


His results involved only the parameter kj = 8 cota. Hemsch (1988, 1989) has utilized this in 
the form below which exhibits a more appropriate asymptotic behavior for small a, 


k3 = 


tana 

Aspect Ratio 


constant 

kl 


(1.9) 
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to correlate data accumulated over a range of thin delta wings. Supersonic flows over sharp- 
edged wings were found to correlate very well for a wide range of flows, up to aspect ratios of 
four, thus extending the similarities proposed by Sychev. Normal force coefficients and 
circulation values generated numerically on a series of thin gothic wings and delta wings with 
sweeps of 65° to 85° were found to correlate well when scaled in this manner. Pitching 
moments correlated less weU, indicating a dependency on aspect ratio with subsonic Mach 
numbers. Interestingly, Eamshaw and Lawford (1964) found the positions of the primary 
attachment line and the secondary separation line to correlate well with atanA, similar in form to 
1/ki , implying the flow is roughly conical. The present study examines extensions of these 

concepts in Section 5.2. 

It is worth noting that the final word on the effects of the leading edge vortex with regards 
to the flow above a delta wing, excluding the breakdown process for the moment, is still not 
definite. Some interesting discussions are put forth by Dixon (1989) on the physics of leading 
edge vortex lift He suggests that the primary contribution to vortex lift is due to the shear layer 
and that the vortex which results from the rollup of this shear layer has only small, indirect 
effects on the lifting surface. His major conclusions indicate that total lift is not strongly 
sensitive to the position or strength as long as the vortex is near the leading edge. In addition, 
the drag depends on the amount of flow entrained by the free vortex, which in turn is a function 
of the free vortex circulation, its rate of growth, and the distance along the vortex axis. The 
following discussion on the mechanisms associated with breakdown will further serve to 
illustrate the diversity of thought on this subject. 
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1 .3 Breakdown Mechanisms and Criteria 

The occurrence of vortex breakdown is not a random event. Certain physical conditions 
must be present in a flow that will eventually exhibit breakdown. First, the flow must be of a 
highly swirling nature. The amount of swirl present in a flow can be determined by the value of 
the swirl angle = tan ' 1 (V 0 /V x ), determined from the local values of the axial, V x , and tangential, 
V 0 , velocity components. Generally speaking, this value is determined in the vortex core and 
found to be greater than 40° for the flow directly upstream of breakdown. Secondly, it appears 
that for breakdown to occur there must also be an adverse axial pressure gradient. The occurrence, 
as well as the position, of the breakdown region can be severely influenced by the pressure 
gradient outside of the immediate vortex flow, such as a convergence of the tunnel walls or an 
object placed in the downstream wake. The removal or reduction of an adverse, axial gradient can 
delay breakdown. Finally, all breakdown phenomena show a divergence of the vortex core region 
directly upstream of breakdown. This can be attributed to the axial adverse pressure gradient, 
however as Hall (1972) mentions, any divergence of the outer streamtubes, even if the external 
adverse gradient is not sufficient to induce breakdown, can produce an adverse gradient in the 
interior of the vortex core because of continuity. 

The breakdown region appears to follow one of two forms: the axisymmetric bubble 
breakdown, where the rearward part of the bubble is open ended and irregular, as if it were 
shedding from some imaginary blunt body of revolution, and the spiral type, typified by an abrupt 
kink in the subcore which proceeds to spiral about the original axis for up to several revolutions 
before dispersing (Figure 1.2b). The axial component of the flow in the core, measured at up to 
three times ffeestream, decelerates and stagnates in the space of one or two core diameters. This is 
accompanied by a substantial rise in the core pressure. Expansion of the cores, as presented by 
Leibovich (1978), range from 1.5 to 3 times the diameter present in the upstream flow, the bubble 
form showing significantly greater expansion. Reversed flow is typically noted in bubble 
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breakdown forms. This condition of reversed flow is not specifically required in order that the 
'breakdown process’ occur, however, as has been pointed out by Cutler, Naaseri, and Bradshaw 
(1989) and Faler and Leibovich (1977). 

The breakdown phenomena is not solely a delta wing phenomena. Swirling tube flows 
also exhibit breakdown characteristics similar to delta wing vortical flows. Another phenomenon 
which has shown evidence of breakdown is the tornado. This, again is different from either a delta 
wing or a confined tube flow in that there are no solid boundaries parallel to the vortex axis. The 
complexities of these vortical flows are such that a complete analytical solution is difficult to 
obtain. Theories do exist that partially predict the occurrence of certain flow phenomena, such as 
breakdown or the sensitivity of the region to the severity of external parameters, such as pressure 
gradients. 

Excellent reviews of research dealing with vortex structures and breakdown, both 
theoretical and experimental, have been written by Hall (1972), Parker (1976), Erickson (1981), 
Wedemeyer (1982), Hall (1985), Leibovich (1984), Lee and Ho (1989), and Ng (1989). These 
deal with all aspects of vortex flows, including tube vortex flows and delta wing vortices, and 
range from stability/breakdown observations to analysis of theoretical models. Information 
pertinent to the present study is summarized below. The following explanations are detailed 
according to their theoretical, numerical and experimental nature. 


1.3.1 Theoretical Suppositions 

Several theories have been proposed to explain the occurrence of breakdown. Hall (1972) 
groups these theories into three categories which will also be used here for the purpose of 
discussion. 
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1 . Vortex breakdown is similar to the separation of a two dimensional boundary layer. 

(Gartshore 1962, Stewartson and Hall 1966 ) 

2. Hydrodynamic instability results in the transition to post breakdown flow. (Ludweig 

1962,1965) 

3. An existence of a critical state occurs which causes breakdown to occur. (Squire 1960, 

Benjamin 1962, 1967 Bossel 1967,1969) 

These will be considered in detail below and a more extensive analysis can be found in Hall 
(1972). 

The quasi-cylindrical approximation for describing the vortex field, using boundary layer 
simplifications, was detailed previously with the theoiy of Stewartson and Hall (1966). The 
idea behind using this approximation as a means for predicting breakdown is the same as that of 
determining separation points in a boundary layer flow. If the calculation shows that a 
separation of the flow may occur in a boundary layer, the corresponding real flow has been 
shown to actually separate at or near that location even though the equations themselves fail to 
model the flow at that point. In the same manner, one might expect that a quasi-cylindrical 
vortex core calculation which began to indicate appreciable axial gradients, or a failure of the 
equations at that point, would in the physical flow correspond to vortex breakdown. Gartshore 
(1962) and Bossel (1971) introduced other assumptions on velocity profiles and even included 
reversed flow regions. Unfortunately, the differences between failure of the quasi-cylindrical 
approximations and the limit of the velocity distributions to describe the flow are difficult to 
interpret. Gartshore concluded that breakdown occurred because of a viscous diffusion of 
vorticity from the core of the vortex to the outer flow. Mager (1972) associated a singularity in 
the incompressible quasi-cylindrical momentum-integral solution with breakdown. 
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The concept of a critical state of the flow serves as a basis for several theories, all of which 
interpret this concept in a slightly different fashion. The ability of the flow to sustain an 
infinitesimal stationary sinusoidal disturbance, referred to as subcritical flow, or an infinite set 
of disturbances, supercritical, is the criterion used for breakdown. If the flow is at the point 
where it can just support a disturbance, it is termed critical. Squire (1960) interpreted this 
critical state in general terms stating that if standing waves can exist, disturbances downstream 
will propagate upstream and cause breakdown. Randall & Leibovich (1973) contend that this 
wave is in fact a finite amplitude wave. 

Benjamin (1962, 1967) proposed that the critical state had an important but not exclusive 
role stating that breakdown was 'a transition between two states of axisymmetric swirling flow, 
being much the same in principle as a hydraulic jump...'. Upstream of breakdown the flow is 
supercritical and cannot support standing waves, while downstream the flow is subcritical and 
can support standing waves. He showed that the subcritical flow could be represented as a 
small but finite perturbation of the supercritical flow. In this way he countered Squire by 
proposing that the existence of these waves do not lead to breakdown, but that the leading 
standing wave is the breakdown phenomenon. Benjamin demonstrated that increasing the swirl 
of the flow causes the supercritical flow to approach the critical state and if increased enough, 
the flow will go subcritical. 

Theoretical extensions of Squire and Benjamin based mainly on linear and nonlinear wave 
models were carried out by many investigators including: Leibovich (1970), Landahl(1972), 
Bilanin (1973), Leibovich and Randall (1973), Randall and Leibovich (1973). This critical 
approach was also exploited by Bossel (1969) who constrained the flow to a state of rigid 
rotation. He then imposed boundary conditions such that a stagnation point existed and found 
the flow to diverge and form a bubble. Larger upstream swirl formed a bubble easier and the 
form of the theoretical streamlines were similar to experimental observations. Bossel used the 
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ideas of Hall to estimate that a swirl angle, defined as tan _1 (V z / V x ), of 54.8° anywhere in the 
vortex would cause flow stagnation. Escudier (1983) goes a step further by proposing a two 
stage, supercritical - supercritical / supercritical - subcntical, transition based on swirling tube 
flows. The first stage of the transition is isentropic. The flow changes state from an initial 
supercritical condition to a state that is also supercritical and can revert back to the upstream 
flow. The second stage of the transition is non-isentropic, much like a hydraulic jump or a 
shock wave, in which the flow transitions to the post breakdown sub-critical state. In addition, 
his analysis suggests breakdown occurs for a unique swirl number, T / k rcU^, which is 

equal to V2 for a Rankine vortex with an infinitesimally small core radius, r c . He does 
conclude, however, that a free breakdown must have a different character than tube flow 
breakdowns. 


Stability arguments for breakdown were proposed as early as 1916 by Lord Rayleigh. His 
analysis indicated that a cylindrical flow with no axial motion was stable provided that the 
derivative of the square of the circulation was positive, 


i acr 2 ) 

dr 


> 0 


( 1 . 10 ) 


Extensions to this concept were proposed by Howard and Gupta (1962) in the form of 


j_ar 2 
r 3 ~ar i 
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“5T 


( 1 . 11 ) 


to maintain stability to axial disturbances for cylindrical flows with axial shear components. 
Leibovich and Stewartson (1983) introduced a criterion based on Rayleigh's arguments with 
axial flow permitted for the centrifugal instability of columnar vortices in the form 
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aafaftar 



< o 


( 1 . 12 ) 


where Q = IDIT/4Q is a dimensionless parameter introduced by Sarpkaya (1971) using the 
volume flow rate Q. This inequality is suggested to indicate massive instability and transition to 
turbulence at large Reynolds numbers if exceeded. 

The hydrodynamic instability ideas put forth by Ludweig (1962, 1967) were based on a 
stability boundary he determined for inviscid flow in a narrow annulus. Accordingly he 
suggested that the vortex flow will be unstable to spiral disturbances if 

(i • ■ (' • ££*)(' - (v^j 2 ) > 0 <*•'» 

Once the flow became unstable amplification of the spiral disturbances might occur which would 
induce an asymmetry on the core, leading to stagnation and breakdown. The explanations for a 
breakdown process involving spiral disturbances are not strong, however. No evidence points 
to the ability of weak spiral disturbances to cause asymmetry of the stable, axisymmetric, 
upstream core. Leibovich (1984) notes that the expansion of the vortex core after breakdown is 
a direct result of the mixing associated with the instabilities and turbulence. Therefore, whether 
these spiral instabilities are responsible for breakdown or not, they do play a role on the post 
breakdown flowfield. 

The stability of the axial vortex structure has even been examined in light of tornado 
structures by Davies- Jones (1982). He comments on the ability of the vortex to sustain inertial 
wave motions because of its 'rotational stiffness’, unlike a non-rotating, homogeneous, 
incompressible flow. Davies- Jones considers a flow where the centrifugal forces are balanced 
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by the radial pressure gradient and the angular momentum increases in absolute value away from 
the rotation center. If a fluid particle displaced radially outward, while conserving its angular 
momentum, it will experience a net force that will restore it to its original position. He notes that 
such a displaced fluid particle will actually oscillate about its equilibrium position until damped 
by viscous forces. 

It has been observed that, when compared with experimental data, both Benjamin's 
criticality theory and Hall's quasi-cylindrical approach do predict the occurrence of vortex 
breakdown accurately in some cases. Neither is complete, however, in describing all the 
breakdown phenomena occurring in vortex tube experiments and aspects of both theories may 
be required, as well as instability arguments, for a complete model. 

A recent study on the physical mechanisms governing vortex breakdown in confined 
cylindrical flows by Brown and Lopez (1988) introduces some interesting possibilities on the 
breakdown of delta wing vortices. Several extensive theoretical and numerical discussions on 
the aspects of this type of flow have been given by Turner (1966), Bode, Leslie, and Smith 
(1975), Rotunno (1979), and Wilson and Rotunno (1986). They deal more specifically with 
tornado flows and the insights on the behavior of vorticity are very valuable. The discussion of 
Brown and Lopez (1988), outlined below, presents a unique criterion based on the relation of 
the angle of the velocity vector to the vorticity vector on a stream surface upstream of 
breakdown. Variation of the problem parameters, especially the Reynolds number, has resulted 
in an axial flow deceleration great enough to sustain a stagnant flow region containing 
recirculating fluid and termed a vortex breakdown bubble. This phenomena is very similar in 
character to that seen in cylindrical tubes and over delta wings and so will be detailed below. 
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Brown and Lopez regard the breakdown region as a transition region from a concentrated 
vortical flow to a solid body rotation. Their basic argument is that the physical mechanisms 
involved rely on the production of negative azimuthal vorticity, that is 


av r av z 

11 = ~3T - nr 


(1.14) 


which results from a tilting and stretching of the predominantly axial vorticity vector, £. They 
further state that steady, inviscid, axisymmetric swirling flow can be viewed in three ways: an 
interaction between the total head and angular momentum of the fluid (both conserved on a 
stream surface), a balance between the radial pressure gradient and the centrifugal force, or in 
terms of the generation of azimuthal vorticity. 


They addressed the question of why the strong vortical core diverges by considering the 
Euler equation of the radial momentum: 



+ V. 
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Downstream from the point where the flow is cylindrical, V r = 0, a divergence of the flowfield, 

av r 

> 0, will occur only if the centrifugal force exceeds the pressure gradient (in the absence of 

viscous stresses). This process can be examined further from a vorticity point of view. Using 
circulation (angular momentum) and the total head respectively: 

P (V r 2 +V z 2 +V 2 ) 

r = rv„ h=^+ — — I — *- 


and substituting into the above equation of motion gives: 
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( 1 . 16 ) 


Employing the azimuthal vorticity, T|, defined earlier yields: 
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( 1 . 17 ) 


Rewriting this in terms of the stream functions: 




keeping in mind that T and H are constant on y yields: 
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Further manipulation leads to: 
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( 1 . 19 ) 


where o is the radius of the stream surface r = o(z) for a curve on which the stream function is a 
constant (i.e. y(r,z) = yp. The ter^s a 0 and p 0 are defined as: 



and represent the tangents of the helix angle for the velocity and vorticity vectors rt^^vely. 
The term £ is the axial vorticity and the 0 subscript denotes a station upstream of the breako*^ 
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The essence of the above derivation is that T) is dependent on the ratio of a 0 to p 0 and a to 
o Q . Brown and Lopez suggest that if an inviscid mechanism is the dominant factor, then a 
negative value of q is required to bring the axial velocity, V z , to zero. The value of T| only 
becomes negative for a 0 > p 0 as the ratio of a to a 0 increases from unity, that is on a 

diverging stream surface. This negative q will induce a negative axial velocity leading to a 
further increase in a and consequently a more negative value of T| . It is this ’positive feedback’ 
which may be responsible for the rapid spatial divergence of the core. Diverging of the core will 
continue until the negative q grows large enough to actually turn the flow back towards the 
axis. 


l n and Lopez then comment on how this type of mechanism would apply to pipe 
flows >ossibly to freestream swirling flows like delta wing flowfields. They state that 
upstream of the breakdown the flow is of a cylindrical nature and hence: 


it _ i ap 
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( 1 . 20 ) 
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They further anticipate it is this vorticity diffusion which leads to a radial redistribution of 
T and a stretching and tilting of vortex lines due to an increase of V 0 . This is followed by a 

reduction in the initially positive azimuthal component of vorticity with axial distance and the 
subsequent beginning of an 'inviscid breakdown' process. Leibovich (1984) also supports the 

concept of breakdown being essentially an inviscid process. A criterion for this breakdown 
process is that a 0 > (3 0 and although viscosity is required to initiate a reduction of T] and V 0 and 

the initial divergence of the stream lines, the breakdown can be considered an inviscid process 

driven by the feedback mechanism described earlier. On delta wings, the radial and axial 
velocity distributions will establish a characteristic (Xq / P 0 f° r each angle of attack and be 

critical to the occurrence, location and strength of breakdown. 

Although Brown and Lopez postulate that the breakdown processes is inviscid, the 
viscous interactions of the process can not be completely dismissed. Krause has performed an 
order of magnitude analysis that indicates that the increase in pressure in the axial direction, as 
the axial core flow stagnates, is a direct result of the viscous as well as inertial forces present in 
the core. Provided that the radial velocity components are small, this becomes a solely viscous 
process. He further shows that the outer portion of the vortex cannot support a pressure 
gradient along the axis if the radial distribution of the azimuthal (tangential ) velocity component 
is equal to that of a potential vortex (i.e. a 1/r distribution). 


A final remark is included here on studies performed in the area of tornado research, which 
has developed virtually independently from the research on delta wings. Similarities exist 
between these type of flows and those above a delta wing. Vortex breakdown phenomena have 
been observed, under laboratory conditions, in tornado vortex chambers (TVC) capable of 
simulating atmospheric tornados. These type of flowfields are generally characterized by a swirl 
parameter defined as 


where r 0 and Q represent the inlet opening and flow rate of the tornado apparatus (Davies-Jones 
(1973)). The circulation, T, is determined in a plane normal to the vortex axis. Depending on 
the geometry of the apparatus, this can then be simplified further to the swirl ratio, 

mentioned earlier. An analysis by Walko and Gall (1986) indicates that breakdown of a tornado 
in a TVC is a direct result of viscous diffusion. They note that the axial pressure gradient is zero 
for highly swirling flows, but is strongly influenced by diffusion at lower swirls. In addition, 
the flow is strongly sensitive to the condition of the axial velocity component. Staley (1985) 
proposes that the principal source of kinetic energy for amplifying non-axi symmetric 
perturbations on a tornado flow is the radial shear of the axial flow, which, as noted previously, 
is the principal component of azimuthal vorticity for these types of flows. The axial vorticity 
term, which is related directly to the radial shear of the rotating fluid, then extracts energy from 
these perturbations. Further reviews of these studies can be found by Rotunno (1977), 

Howells, Rotunno, and Smith (1988) , Pauley (1989). An extensive discussion on the physical 
phenomena associated with tornado morphology is given by Davies-Jones (1982). 


1.3.2 Numerical Approaches 

The advent of super computers have enabled researchers to extensively investigate 
numerical models of the three dimensional flows above a delta wing. The majority of numerical 
codes, excluding Navier-Stokes models, may be classified as either potential or Euler codes, 
using either a continuous vorticity distribution in the field or a vortex sheet method, as outlined 
by Hoeijmakers and Rizzi (1985). The former uses vortex sheets with a discontinuous velocity 
potential across them to model the free shear layers. Isolated line vortices are used for the cores. 
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The topology of the flow must be known quite well in advance as the rotational flow regions 
must be "fitted". The method is well established and the results for flows without breakdown 
are in reasonable agreement with experiments. Euler solutions are based on the Euler equations, 
allowing rotational flow everywhere, and capturing implicitly the vortical parts of the flow as the 
solution evolves. Two schools of thought exist with respect to the validity of utilizing these 
inviscid codes. A comparison of the results from various Euler solvers by Williams, Kordulla, 
Borsi, and Hoeijmakers (1990) with experimental data concluded that Euler methods are 
unsatisfactory in representing the flow due to their failure to represent secondary separation. 
Wagner, B., Hitzel, S. M„ Schmatz, M.A., Schwarz, W., Hilgenstock, A., and Scherr, S. 
(1988) reported, however, that reasonable agreement for overall pressure and force distributions 
between experiment and Euler code simulation is possible. 

Predicting breakdown using Euler codes does appear to be possible, however. O Neil, 
Barnett, and Louie (1989) have demonstrated a vortex breakdown effect above semispan models 
of 60° and 70° delta wings at Mach = 0.2 that closely follows experimental trends. They 
conclude breakdown is primarily governed by inviscid factors. Hitzel (1988) concludes that the 
Euler calculations indicate breakdown to be triggered by adverse pressure gradients and is 
independent of viscous effects. Further discussions on Euler simulations can be found in 
Hitzel, Wagner and Leicher (1986). 

Navier-Stokes (NS) codes allow calculations to be performed in a less restricted way than 
the Euler equations. Liu and Hsu (1987) give a review of NASA contributions to three 
dimensional, incompressible, NS simulations of slender wing vortices. Krist, Thomas, Sellers, 
and Kjelgaard, for example used a thin layer approximation of the 3-D, time-dependent, 
compressible, laminar NS equations to compare with LDA measurements by Kjelgaard and 
Sellers (1990) about a 75° delta wing at an angle of attack of 20.5°, Re = 500000 and Mach = 

0.3 . Even though the equations modeled a laminar flow, comparisons at a chordwise station of 
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70% were quite favorable. The maximum predicted velocity, however, was still found to be 
less than that determined experimentally. 

Numerical breakdown results have been presented by Grabowski and Berger (1976) on an 
unconfined viscous vortex, using solutions of the full, steady, axisymmetric Navier- Stokes 
equations. They concurred with Hall (1972) in that breakdown is the result of a critical 
retardation of the flow. As well, the results showed that a vortex with sufficient swirl can be 
reduced to the critical breakdown state by diffusion of vorticity, a non-linear coupling of the 
axial and swirl velocities, flow divergence and pressure forces. 

Ekaterinaris and Schiff (1990) used a Navier-Stokes code on a 75° sweep delta wing for 
32° < a < 40° to predict breakdown. Bubble type breakdown was observed on fine and coarse 
grids with no trace of unsteadiness. Indicators of breakdown, such as reversed axial flow and 
increases in surface pressure, along with diverging and coiling streamlines substantiated their 
arguments. Other NS studies which have predicted bubble and spiral breakdowns can be found 
in Fujii and Schiff (1989), Thomas, Taylor and Anderson (1987) and Hartwich, Hsu, 

Luckring, and Liu (1988). Axisymmetric vortex filament methods by Nakamura, Leonard and 
Spalart (1985,1986), employing experimental data for upstream conditions, have also been used 
to simulate both bubble and spiral forms of breakdown such as that observed by Faler and 
Leibovich (1978) 

The use of non-dimensional parameters, in conjunction with Navier Stokes codes, has been 
shown to be very useful in predicting flow conditions and vortex breakdown. A reduced form 
of the steady, incompressible NS equations utilizing the Sychev (1960) parameter discussed 
previous, was computed by Dagan and Almosino (1989). Very good agreement with 
experimental data on delta wings with aspect ratios of 0.25 to 0.7 and a up to 30° was obtained 
provided the basic slenderness constraint on both the wing and the vortex was met. The major 
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assumption used a viscous, inner region dominated by the vorticity equation that was 
subsequently matched to an outer potential region. 

Numerical experiments by Powell and Murman (1988) found that the level of the total 
pressure loss predicted in the vortex core was independent of the Reynolds number using a 
similarity Navier-Stokes solution. Their results matched the analytical results of Hall (1961) 
quite well, not all that surprising considering they restricted the flow to an incompressible, 
axisymmetric, conical, high Reynolds number flow with a slender core. Experimental 
comparison with Eamshaw (1961) showed an overprediction of the axial and tangential 
velocities as well as the radial total pressure distribution. 

The effects of Reynolds and Rossby number parameters were investigated by Spall and 
Gatski to evaluate the NS equations formulated in terms of velocity and vorticity. The Rossby 
number is important in the study of the Coriolis force and inertial forces in large scale 
atmospheric fluid motions (Bode, Leslie and Smith (1975)). Although the Coriolis force is 
regarded as a restoring force, it can cause fluid particles to overshoot their initial locations and 
even set up waves in the resulting fluid motion termed inertial waves . The Rossby number is 
defined as 


where r* is the radius of maximum swirl, U* is the axial velocity at r*, and Q is the vortex 
rotation rate as r— >0. If this is applied as a stability criterion to the theory of Squire (1960), R 0 

must be less than 0.56 theoretically to permit the existence of axisymmetric standing waves. 
Spall and Gatski found that since the local Rossby number should decrease near breakdown, 
and that the numerical tests reveal viscous diffusion increases R 0 , it is necessary to impose an 
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external adverse pressure gradient to decrease R 0 with downstream distance, in order for 
breakdown to occur. 


A parametric numerical study by Delery, Pagan and Solignac (1987) of the NS equations 
using velocity profiles of the form 


V* = 1+(V X _ -l)e' 


(-Br 2 ) 


v axis 


v 0= 


r 0 (l- e- Br2 ) 


(1.24) ,(1.25) 


indicated that there is a maximum value of the circulation, r 0’ which will cause the flow to 
transition to breakdown denoted by a negative value of V x in the computational domain. 

Breakdown was concluded to be an inviscid process and only a function of the strength of the 
vortex, provided the Reynolds number is large enough and the other parameters of the flow are 
held constant. A large sensitivity of the onset of breakdown to the axial velocity on the vortex 
centerline was also noticed. 

Presently, it seems that none of these ideas are sufficient to accurately predict vortex 
breakdown on a delta wing over a wide range of conditions. Computational results have been 
seen to correspond to experimental data, however no theory exists which can yield the flow 
detail in the breakdown zone nor universally predict breakdown locations which consistently 
compare with experimental results. In addition, the experimental data is not always completely 
consistent across different investigations. Thus it is essential that further studies be conducted 
in an effort to uncover information which may shed light on this problem. 
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1.3.3 Experimental Studies 

Researchers have conducted a wide variety of investigations to measure the vortex 
flowfield and to observe the onset of breakdown. A summary of relevant studies and 
contributions is presented here. Discussion of results directly comparable to the present study 
will be made, where applicable, in Chapter 5. A list of specific measurements on delta wings 
and their associated geometries can be found in Appendix A. 

Vortex breakdown was first detected by Peckham and Atkinson (1957) over a highly swept 
delta wing at large angles of incidence. This phenomena was also seen to occur elsewhere as 
noted by Smith and Bessemer (1959) in aircraft trailing vortices. Since parameter variation and 
flow control could be defined much easier in confined flows, a great majority of tests have been 
conducted in swirling pipe flows. The asymmetries present in delta wing flows are removed 
and the numerical results are more comparable. The results of these studies will be detailed 

first. 

1.3.3. 1 Tube Flows 

Some of the first experiments on swirling pipe flows were conducted by Harvey (1962) 
using flow drawn radially inward through a set of vanes and into a tube. He observed that 
variation of the swirl of the flow indicated breakdown to be an intermediate stage between two 
types of flow: those that exhibit axial velocity reversal and those that do not. Harvey also 
concluded that, because the observed breakdown could revert back to the upstream conditions 
via the introduction of a favorable pressure gradient, the reversible breakdown process 
represents a division between subcritical and supercritical regimes and is not an instability 
process. An instability would grow unchecked and lead to an irreversible process. 
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Analytical profiles of the form, 

ar 2 

v x = v xl +v x2 e<_ar2> V fl = ^ r * (1-26) , (1.27) 

equivalent to that presented by Batchelor (1967) above, fit Harvey's (1962) data quite well, as 
noted by Leibovich (1984). Recall that this is the same analytical profile as that in Figure 1.5. 
The major difference between these types of flow and that above a delta wing, apart from the 
physical wall boundary, is that the core size is a function of the boundary layer shed by the 
centerbody of the apparatus. Hence, both the flow rate (Reynolds number ) and the vane angle 
contribute to the amount of vorticity in the core, unlike sharp-edged delta wings, where 
Reynolds number effects appear to play a negligible role above a critical value. 

The effect of instabilities on stationary breakdown was investigated by Sarpkaya (1971) in 
a series of experiments in a conical tube apparatus. He noted from motion pictures that a 
toroidal vortex ring, with an axis gyrating at a regular frequency about the bubble axis, appears 
to empty and replenish fluid in the breakdown bubble. Variation of swirl and Reynolds number 
produced spiral and bubble forms, including a double helix form with no stagnation point. 
Sarpkaya (1971) felt breakdown did not depend on a single mechanism, but resulted from 
hydrodynamic instabilities, depending on the combination of Reynolds number and circulation 
of the flow. He concluded (1971) that his observations are in perfect accord with Benjamin's 
theory, and no other theory can account for this itself. Further experiments on the effects of 
adverse pressure gradients by Sarpkaya (1974) indicate that the effect of increasing this 
parameter has the same result as increasing the circulation or the mean flow rate. 

Six distinct types of vortex tube breakdown modes were isolated by Faler and Leibovich 
(1977), including those categorized by Sarpkaya (1971,1974). A schematic of these forms of 
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breakdown is given in Figure 1.15. For all the breakdown cases observed, the flow was 
laminar upstream of breakdown and turbulent shortly downstream of the breakdown region. 
Increasing the Reynolds number, based on tube diameter, or the swirl angle, caused each type 
of breakdown to be seen in a repeatable succession. The highest Reynolds number (about 8000 
based on the tube diameter) and swirl angle gave rise to the bubble form of breakdown. Type 0 
refers to the so-called axisymmetric mode. A stagnation point on the axis is followed by an 
abrupt expansion around a bubble of recirculating fluid. For low circulation rates and higher 
Reynolds number, the type 0 form alternated with the type 1, which has a slight asymmetric, 
ragged nature that tends to break up into a turbulent wake. Four others types representing a 
variation on the spiral structure. The spiral type of breakdown, type 2, occurred at lower 
Reynolds numbers and is marked by an abrupt kink of the fluid along the vortex axis. At the 
lowest values of swirl and Reynolds number, the type 6 form appeared. The vortex filament 
was moved gently off the axis, expanding in radius, until it nearly reached the tube wall. The 
filament would also, at times, shear into a tape as it moved off axis. The type 5 form evolves 
directly from the type 6 when another shearing branch forms and they wind around each other in 
a double helix fashion. The type 6 could also evolve into a type 4 where the filament or tape 
would roll back onto itself up to the point of deflection. After a time fluid would be ejected from 
this recirc ulatin g zone and the type 4 flattened bubble form would be fully developed. Type 3 
was a combination of 4 and 2. It should be noted that types 4, 5 and 6 only occurred at 
relatively low tube Reynolds numbers of about 2500. Faler and Leibovich (1977) draw several 
strong conclusions from this study. First, all flows that exhibit axisymmetric, type 0 and 1, or 
spiral, type 2, breakdown are supercritical upstream in the sense of Benjamin (1962). 

Secondly, flows that subsequently undergo breakdown are all stable to infinitesimal, inviscid 
disturbances, including axisymmetric and non-axisymmetric perturbations. Finally, no 
axisymmetric disturbance patterns exist in these flows. 



41 



Type 0 


Axisymmetric forms 


Type 1 






large amplitude, off axis displacement^v^^ ^ 

Type 6 

Figure 1.15 States of Tube Breakdown 

Laser Doppler anemometry experiments provided, for the first time, precise data on the 
velocity field of the vortex and the accompanying breakdown region. Initial investigations were 
carried out by Faler (1976), Faler and Leibovich (1977,1978), Garg (1977), and Escudier et al. 
(1980). Analysis of this data enabled Leibovich (1983) to determine the vortex status at any 
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point as subcritical or supercritical and the results agreed reasonably well with vortex tube 
experiments. Leibovich (1984) argues that these two, visually different forms of breakdown, 
spiral and bubble, are indeed different because of their different expansion ratios and the 
imposed upstream conditions. Sarpkaya (1971) has shown in his conical tube apparatus that 
increasing swirl from zero causes the spiral form to occur first at a particular mean location. At 
some larger value of swirl, the form will switch to the bubble form, several core diameters 
upstream of the previous location. This discontinuity is proposed by Leibovich to be a strong 
reason for considering these two types of breakdown to be truly distinct. He also concludes 
from stability arguments (1984) that the bubble breakdown contains less stable wake and 
approach flows than the spiral form. 

Laboratory experiments have also been carried out in cylindrical tanks having a jet issuing 
fluid tangent to the outside wall to simulate vortex behavior. Velocity profiles measured by 
Escudier, Bomstein and Maxworthy (1982) were found to be qualitatively similar to those 
measured above delta wings. A profile of the form 

V - Is- (l-e-Aro 2 ) + H (1.28) 

0 2 fir 1 

was found to fit the data satisfactorily. It was noted that there was no flattening of the axial 
profile. The implication of this is that since 3Vx/3r still existed, so did the azimuthal vorticity 
component and, thus, there was no external irrotational flow present anywhere. 

Further tests by Escudier and Zehnder (1982) with their tangential jet entry device have 
shown that a simple parameter, ReQ 3 R = constant, R being a dimensionless parameter based on 
the geometry, correlates the conditions for breakdown at a fixed axial location in their apparatus. 
In the swirl vane generators, this becomes ReQ^ = constant. This formulation breaks down, as 
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pointed out by Leibovich (1984) since bubble and swirl breakdowns have different locations as 
a function of Re and co in swirl vane devices. The most interesting aspect of the results of 
Escudier and Zehnder is that they noted an opposite sense of rotation of the spiral rotation to the 
rotation of the base flow. As mentioned previously, the spiral form of the breakdown continues 
to wrap about the core axis for several revolutions. Experiments in swirl vane generators have 
shown the sense of rotation of the tracers in the flow to be the same as that in the upstream flow. 
Spiral breakdown of leading edge flows, as shown by Lamboume and Bryer (1961), have a 
sense of rotation of the core filament, in the same manner as the tangential jet entry device, that 
is, opposite to that of the upstream swirling flow. 


1. 3.3.2 Delta Wing Flows 

Many studies have been performed on delta wing planforms. Their empirical behavior is 
well documented under a range of various conditions. Based on research cited in the literature, 
Payne (1987) has reviewed the major factors which can influence the breakdown of vortices 
above delta wings. This is summarized below followed by a brief discussion on particular 
studies relevant to the present efforts 

The actual position of this breakdown is a strong function of the pressure gradient along the 
vortex, the initial axial core velocity, and the angle of sideslip, or yaw angle. Increasing the 
sweep angle or decreasing the angle of attack causes the location of the breakdown to move aft. 
The breakdown position will move forward during flow acceleration and remain so until the 
steady speed condition is reached, whereupon it returns to its normal breakdown position as 
noted by Lamboume and Bryer (1961). The reverse is true for deceleration. An increase in the 
swirl of the flow or a larger adverse pressure gradient tends to promote the onset of breakdown. 
Thicker wings, rounded leading edges, lower Reynolds numbers, and of course more complex 
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geometries can also substantially influenced the location of breakdown. Comparison of the 
seven hole probe and LDV data acquired by Payne (1986,1987) with the present study will be 
made on an individual basis in Chapter 5 . 

Payne also summarized the effects reported to have a minimal effect on breakdown. 
Reynolds number effects on the flow are small at higher Reynolds numbers, but the effects on 
position and strength of the leading edge vortices become more pronounced at low Reynolds 
numbers, that being below Rec - 100,000. This applies to sharp edged, thin delta wings at 
moderate angles of attack. A study using a flat, 63° sweep delta wing by Schrader, Reynolds, 
and Novak (1988) revealed the major influence of Reynolds number to be in the viscous 
secondary separation region and that the overall aerodynamic forces were influenced only 
slighdy. Wing stall characteristics were weakly dependent on Reynolds number, but strongly 
affected by Mach number. A slight lift decrease was noted at the higher Reynolds numbers, but 
the slope of the lift curves remained unchanged. Extensive water tunnel tests by Erickson 
(1981) have shown vortex generation, sheet and core location, as well as vortex strength are 
accurately reflected in a water tunnel due to the insensitivity of the separation point location to 
Reynolds number changes. Erickson (1982) notes, however, that correlations are best at high 
angles of attack. 


The geometry of the wing can play a key role in determining the resulting delta wing field. 

A thickening of the delta wing has been seen to reduce the strength of the leading edge vortices 
and the non-linear lift components as noted by Peckham (1958) and Squire (1967). The loss of 
lift is associated with a decrease in dCL/dct resulting from a weaker vortex system. The angle 
of incidence for a certain lift therefore increases and a smaller lift to drag ratio results. Squire 
(1967) did note an increase in the stability as the loss of lift occurred over the forward part of the 


wing. 
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Hummel and Srinivasan (1967) found that increasing the aspect ratio of delta wings causes 
the effects of breakdown on the lift and moment to occur at lower angles of incidence and that 
these effects are decreased as the aspect ratio increases. They also feel the cross sectional shape 
of even thin wings has a considerable influence on the breakdown position. Zohar and Er-El 
(1988) note that the effect of breakdown on the suction induced by the leading edge vortices is 
lower for higher aspect ratio (lower sweeps), while the lift to drag ratio is seen to increase. This 
was deduced from surface pressure measurements on delta wings of 55°-75° sweep. Eamshaw 
and Lawfoid (1964) noted a gain in the lift with a convex surface on the suction side of the 
wing. Lamboume and Bryer (1961) demonstrated how longitudinal camber can delay 
breakdown, presumably because of the beneficial pressure gradients induced . 

Conversely, Wentz and Kohlman (1971) observed that variation of the trailing edge 
geometry of 70° sweep delta, diamond and arrow planforms, had a negligible effect on vortex 
breakdown location. Previously (1969), they noted that the vortex breakdown region crossed 
the trailing edge at the same angle of attack for all these trailing edge configurations. Thompson 
(1975) noted the effect of cropping a 75° sweep delta by 40% resulted in breakdown location 
differences of ± 2.5% x/c in a water tunnel. 

The leading edge shape of the delta wing has also been seen to have a considerable effect 
on the position of the vortex breakdown region. As early as 1955, Bartlett and Vidal 
determined that bevelled edge wings produced a higher value of CL a than rounded edges. 

Squire (1967) determined that a blunter leading edge causes a decrease in the lift curve slope on 
sharp edged delta wings with an aspect ratio of 4/3 (A = 71.6°). He found the loss of lift was 
large and occurred near the apex as indicated by the change in pitching moment. Squire also 
comments that the scale of the vortex flow at the apex is on the order of the wing thickness and 
that the blunting retards the initial vortex development, whereas further aft the vortex is too large 
to be influenced by leading edge details. Kegelman and Roos (1989) investigated leading edge 
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geometries, ranging from blunt rounded edges to leeward and windward bevels, and found the 
breakdown location to be strongly affected by these geometric variations. Sharp edged wings 
with a positive leading edge camber were seen to give the greatest vortex lift and the strongest 
post stall behavior. The breakdown position does seem to be independent of the turbulent 
breakdown of the shear layer near the leading edge. Attempts at altering the leading edge by 
Lamboume and Bryer (1961), including a trip wire on the upper surface had no significant 
effect. 


This breakdown of the vortex structure is seen to occur not so much at a point, but over a 
region. Lamboume and Bryer (1962) describe three stages of vortex breakdown: flow 
deceleration, spiral deflection, and breakdown to full scale turbulence. Measurements by 
Payne, Ng, and Nelson (1987) and others have found the vortex core region to indeed 
transform from a jet-like to a wake-like flow over a spatial region when breakdown occurs. The 
actual location of the breakdown region also tends to vary, within 5% of the chord. Faler and 
Leibovich (1977), in vortex tube experiments, describe the location of the breakdown as only 
quasi- steady in the axial location. The breakdown almost continually drifted back and forth 
along the core axis. The axial extent and velocity of the drift was seen to increase with 
increasing Reynolds and circulation numbers. 

The pressure field above the delta wing and the external pressure gradients imposed by the 
test section have been an area of extensive testing. Pressure measurements in the vortex core by 
Lamboume and Bryer (1961) for the 65° sweep wing at ft = 15 indicated that while both the 
static and total pressure fell as distance aft of the apex was increased, and most sharply at the 
apex, the total pressure was seen to have a nearly constant value along the length of the vortex at 
roughly Cpt = -5.0. Lamboume and Bryer state that the pressure distribution along the axis 

depends on three factors 

i) increasing vortex strength tends to provide a falling axial pressure 
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ii) diffusion of vorticity in the viscous core provides a rising pressure 

iii) deceleration of the longitudinal velocity component in the irrotational flow region tends to 

cause a rise in the axial pressure. The trailing edge can cause such an effect because of 
the pressure recovery experienced at the trailing edge. They suggest the pressure 
changes sensed in the core are amplified above what is experienced in the outer flow. 

In contrast to the data of Lamboume and Bryer (1961), Naarding and Verhaagen (1988) 
reported a drop in the axial total pressure with increasing x/c as well as a sharp drop near the 
apex. They comment that the difference is because their probe was on the order of the viscous 
core size and Lamboume and Bryer's probe was an order of magnitude smaller, due to the 
larger wing used by Lamboume and Bryer. The better resolution of Lamboume and Bryer 
would explain their lower pressure values, but does not explain the reduction of values as 
Naarding and Verhaagen move closer to the apex. Smearing of the pressure profile due to a 
locally large probe size, would result in higher pressures, not lower. 

Kegelman and Roos (1990) noted the total pressure loss in the core of 60° and 70° flat plate 
delta wings with a 25° bevelled windward edge to be the same at all chordwise locations for a 
constant angle of attack. The peak loss value was seen to increase with angle of attack and the 
60° vortex structure was concluded to be stronger than the 70° because of a higher peak 
vorticity. Surface pressure measurements also indicated a stronger suction peak for the 60° 
wing until the occurrence of breakdown. The effect of Reynolds number was examined and 
found to have almost no effect on core trajectory, breakdown location or aerodynamic loading. 
At high a, C Lmax is seen to decrease slightly. Below this, the net loading does not change for 

increasing Reynolds number despite an outboard shift in the secondary separation point. 

LDA measurements were made by Pagan and Solignac (1986) on a vortex generated by a 
75° sweep delta wing at a = 19.3° and allowed to enter a two dimensional variable pressure 
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duct. They noted large amplitude low frequency oscillation in the region of flow just in front of 
the breakdown bubble. They also determined that the flow behaves in an Euler form outside 
and inside the recirculation bubble by computing the difference between the pressure gradient 
and the cross product of the velocity and the vorticity vectors, the only terms in the Euler 
approximation of the equations of motion. A difference of zero indicated the flow to be 
behaving in an Eulerian sense. This condition was not verified, however, along the relatively 
thin interface layer at the outer boundary of the recirculation zone. This interface boundary also 
exhibited very high axial velocity fluctuations. Further tests to examine adverse pressure 
gradient effects were made by Delery, Pagan and Solignac (1987) with the same apparatus at oc 
= 27.5°. A close correspondence was found between the pressure gradient and the vortex 
strength as given by the local swirl velocity in Figure 1.16. 



Figure 1.16 Effect of Adverse Pressure Gradient and Swirl on Breakdown 

(Delery, Pagan and Solignac, 1987) 


A limit to the vortex strength was found, beyond which breakdown occurred even in the 
absence of an adverse gradient. Further analysis of the data at 19.3° revealed a breakdown 
oscillation of about 8 Hz. Earnshaw and Lawtord (1964) observed low frequency force 
fluctuations for sweep angles less than 65° at low angles ot incidence. Higher sweep angles 
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reduced these fluctuations and if 65° < A < 76° no fluctuations were seen until maximum lift was 
achieved. 

The stability concepts put forth by Ludweig (1961) were concluded to be experimentally 
supported by Engler (1988). He states that Ludweig's stability theory is correct based on an 
opto-acoustic technique of measuring the stability parameter in the flow field above a 68° sweep 
delta wing. Row visualization by Lowson (1988) on 70° and 80° sweep wings at Reynolds 
numbers of 3000 to 20,000 points to quasi-two dimensional instabilities existing in the shear 
layer leaving the leading edge. A second local streamwise instability was also seen to be present 
and both instabilities were inhibited by the vortex stretching and wrapping process. The 
frequency of this leading edge shedding followed a Re 0 '* 5 law in the same form observed by 
Gad-el-Hak and Blackwelder (1985). Ng (1989) points out, however, that the boundary layer 
thickness is varying continuously along the leading edge. Since the shear layer thickness is 
dependent on the boundary layer thickness just prior to separation, and the characteristic velocity 
is constant along the leading edge, any shedding frequency should vary according to position. 

Studies concerning the secondary vortex have not been nearly as extensive as those 
centered on the primary structure. The influence of the secondary vortex structure and its 
separation location on the overall flow has been investigated by Hummel (1979), who carried 
out extensive tests on an aspect ratio 1 (sweep = 76°) flat plate delta wing including balance, 
pressure and boundary layer measurements. He reported that the presence of the secondary 
causes a displacement of the primary vortex inwards and upwards and that the presence of the 
secondary vortex locally increases suction on the surface below it. These effects are small for 
turbulent boundary layers, but large for laminar. Also, Hummel noted that the trailing vortex 
forms in a rotation sense opposite to the leading edge vortex. 
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Kjelgaard and Sellers (1990) obtained a series of measurements above a 75° sweep wing, 
including LDV, five hole probe and surface visualization. They confirmed the transition of the 
boundary layer to occur at a Reynolds number 

Re. = 0-29) 


where is the streamwise distance from the apex to the transition point, of Hummel (1979) as 
being approximately 800,000 - 900,000. The transition of the boundary layer from laminar to 
turbulent shifts the secondary separation point outboard. Naarding and Verhaagen (1988) found 
this transition Reynolds number limit to extend to about 2 million for their biconvex wing at 
lower incidences. Carcaillet, Manie, Pagan and Solignac (1986) noted the beginning and end of 
the transition region which, from their tests lies roughly between that of Hummel and Naarding 
and Ver haa gen Naarding and Verhaagen note that for wings of aspect ratio 1 (A = 76°), if the 
Reynolds number based on the root chord is below 500,000 the entire boundary layer on the 
upper wing surface can be expected to remain laminar. 

Increasing the Reynolds number was seen to move the mean boundary layer transition 
point towards the apex as reported by Carcaillet, Pagan and Solignac (1986) from tests on a 75 
sweep flat plate delta. They comment on the non-conical nature of the flow, this being much 
more evident in the core pressures than the velocities. Measurements with a 3-D LDV indicated 
large turbulence levels in the primary, the feeding sheet, and large shear stresses in the region 
between the feeding sheet and the secondary vortex. 


Breakdown correlations are continually being offered as a means for prediction of 
breakdown. The swirl angle was one such parameter mentioned previous. Hawk, Barnett, and 
O'Neil (1990) further analyzed the data acquired by Kegelman and Roos (1990). They point 
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out that conventional models of the fields would give swirl angle distribution contours as 
concentric rings, growing in magnitude away from the vortex model, and reaching a maximum 
at the outer edge of the core. No such pattern was evident They observe most local swirl angle 
values are greater than 45° in the vortex burst region. These maximum swirl angles tend to 
occur in the shear layer and not the rotational core, maintaining constant values along the wing, 
even though the vortex strength is increasing. In addition they appear to be independent of the 
wing sweep, angle of attack, or loading. 

Cornelius (1990) has examined the effect of the Rossby number and a second Rossby 
parameter, defined as 


Rft 


Uc 

ft r v 


Ju^ dy dz 

Rmr = — (130) . (1.31) 

JV 2 + W v 2 dy dz 

A 5 


respectively, where r y is the radius of the core vorticity = y A^/k where ft is the rotation rate 
defined from the integrated axial vorticity, U c is the axial velocity, is the area encompassing 
the axial vorticity, and V v and W v are the mean velocities in vortex coordinates. He concludes 
that Rq > 0.55 and Rmr > 1.0 will lead to a growth of spiral instabilities while for Rq > 0.55 
and Rmr < 1.0 the bubble form is imminent. Note the difference in Cornelius’ definition of this 
parameter from that of Spall and Gatski. The Rossby number has also been defined for 
atmospheric vortex structures as 

^ 0 max maximum tangential velocity 

Rossby Number ^ f " average vorticity * core radius for V 0 max 

c 


where values of 10^ for tornados and 10^ for mesocyclones have been reported by Davies- 
Jones (1982). The implications of the Rossby parameter will be examined further in Chapter 2. 
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Lamboume and Bryer (1961) have conducted tests in wind and water tunnels using various 
flat and cambered plate geometries and noted a correlation of the breakdown position with the 
parameter involving the wing sweep and the angle of attack, defined as y = cos'kcosa sinA). 
This parameter indicated a correlation of the data taken on their swept wings of 55° to 70°. They 
suggested, on the basis of their tests, that breakdown is not based on the amplification of small 
upstream disturbances. They also determined that the total pressure in the core and the adverse 
pressure gradient along the axis, which can be altered by the geometry of the wing or externally 
to the field, are essential factors in causing breakdown to occur. Although no measurements 
were made, the proportionality of the core radius to \ vc/U was proposed. Eamshaw and 
Lawford (1964) noted that at moderate incidences, the location of the surface flow separation 
and attachment lines were seen to correlate in a linear fashion with atanA, indicating the flow 
behaves in a roughly conical fashion. 

Finally, a swirl type correlation for experimental data based on the parameter 
„ (dw/dr)r_o r Wmax 1.994 Hf 

s ■ V P (U2) 

which is defined on conditions in the core, was presented by Wilson (1977). The parameters k 
and X were supplied by a subsonic potential flow panel method. The dependence of this 
parameter on the angle between the wing leading edge and the ffeestream, 0 = cos' ^ (cosa sinA) 
was seen to follow a bounded linear distribution. 

Many efforts have been aimed at uncovering the physics behind the vortex breakdown 
phenomena. The ability to predict the location of vortex breakdown, theoretically, numerically or 
from experimental data with a high degree of proficiency is very desirable. As can be inferred 
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from the above discussion, investigative results, both experimental and numerical, are not always 
conclusive and have even led to contention on particular aspects of the problem. There are certain 
aspects, however, upon which every one agrees. 

The flow upstream of breakdown on a delta wing can be modeled well by a number of 
different methods, both analytically and numerically. Comparison with experimental data away 
from the apex and trailing edge regions indicates very good agreement The tangential velocity 
profile, for example, can be approximated quite well with a simple exponentially declining 1/r 
distribution. Aerodynamic forces can also be predicted very well within a certain angle of attack 
regime. Overall, investigators agree the delta wing vortex flow follows a conical fashion. The 
subcore, however, is often assumed to be cylindrical when used in computational schemes, but 
this assumption has yet to be suitably quantified empirically. 

The prediction of the breakdown phenomena and the resulting changes in the flowfield 
structure does not compare consistently with empirical data. As can be expected, breakdown 
prediction is the area where the most disagreement arises with respect to what is physically 
occurring. Numerical results are typically qualitative both in location and flow features. It would 
appear that the use of any flow approximation short of a full Navier-Stokes code is insufficient to 
model the physics of the flow behavior at breakdown. Euler codes can not account for diffusion, 
boundary layers or transition phenomena and appear capable of only predicting pre-breakdown 
flows with a high degree of confidence. The breakdown itself may behave in an inviscid manner, 
but some type of viscous interaction to approach the state needed to breakdown seems to be 
required. 

Understanding a fluid phenomena, such as vortex breakdown, demands a strong 
understanding of what is physically occurring. The relation of the strength of the vortex to the 
associated axial flow, for instance, is critical to the state of the vortex. The axial velocity is in turn 
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is affected by the external flow conditions which are not always documented in experimental 
studies. Anything that alters the local pressure gradient along the axis and therefore affects the 
axial velocity can have profound effects on breakdown. This ranges from the local planform 
camber to the tunnel pressure gradient Most of the concepts on the flow state prior to breakdown 
are theoretical in form and have not been verified experimentally, but investigators agree that the 
flow must reach some critical form which then transitions to a state that will initiate breakdown. 
Certain physical features of the breakdown have only been addressed only at an empirical level, 
such as the types of breakdown, but these features may have a minimal effect on the overall flow. 
Bubble and spiral forms are seen on a delta wing, yet the post breakdown flow behaves in the 
same manner far enough downstream of the breakdown zone. On the other hand, there could be 
more substantial physics involved. The sense of the spiral breakdown, for instance, is opposite to 
the rotation of the primary vortex on a delta wing, but carries the same sense for a breakdown in a 
tube. 


No standards for theoretical or experimental parameters exist that describe the flow state in 
a manner consistent enough to indicate the onset of breakdown, except for perhaps the swirl angle. 
Even this value is subject to dispute, varying from as much as 42 to 52° in the flow measured 
directly preceding breakdown. In addition, the swirl angle does not account for any local pressure 
field and is a point property as opposed to an integrated field effect, which may not be truly 
representative of the flow. If the combination of certain measured flow features, perhaps in the 
form of some parameter, could indicate the state of the flow at a particular chordwise station, and 
if their values changed in such a way in a downstream direction to approach a limit, breakdown 
could also be predicted. The following chapter outlines this line of reasoning. 



CHAPTER II 


A PATH FOR THIS STUDY 


2. 1 Characterizing the Onset of Breakdown 

The motivation for the majority of delta wing research is to determine why the leading edge 
vortex structure breaks down. The present study falls under this umbrella as well and poses 
several questions as a more concrete basis for direction. What is occurring in the flow that 
prevents the vortex structure from remaining in its pre-breakdown state? For that matter, what 
keeps the vortex in its pre-breakdown state in the first place? Is it possible to define the conditions 
on the flowfield in such a way as to tell when the breakdown process will occur? Are there 
physical quantities that can be measured that will indicate the state of the vortex? 

The following discussion examines these concepts. First, some general comments on the 
use of particular parameters to evaluate the vortex flowfield are made. A specific look is then taken 
at the nature of the vorticity field above a delta wing and what light may be shed on the nature of 
breakdown by a more in depth study of this property. Finally, the specific goals and objectives of 
this study are outlined. 

In broadest terms, the vortex flowfield can be regarded as a transition from one flow state 
to another which may occur as a result of a combination of certain flow parameters reaching a 
critical or unstable state. One can attempt to quantify such critical breakdown parameters in two 
forms: as a function of either the independent or the dependent variables. The former would 
involve factors such as angle of attack, sweep angle of the wing, and sideslip. An indication 
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parameter involving dependent variables would be based upon flow conditions resulting from the 
geometry of the problem such as the local condition of the adverse pressure gradient or the local 
vorticity. 

From the review in Chapter 1, the determination of such breakdown indication parameters 
could be based on several forms including 

i. A local balance of the pressure forces with the acceleration of the fluid. 

ii. A local balance between the generation and convection of vorticity. 

iii. The size of the local length scale compared to the wing geometry. 

iv. A relation based on external pressure gradient, wing geometry, angle of attack, etc. 

The conditions at breakdown can also be used to establish an empirical criterion for the 
onset of this phenomenon. Then, given the relevant parameters of the flow, either from a 
geometrical standpoint (alpha,sweep angle) or a consideration of the flowfield (vorticity, 
circulation, local pressure distribution), the position of breakdown could be determined. Any of 
these factors can also be examined in a chordwise progression to see how they vary up to, and 
beyond, breakdown. The swirl angle criteria is an example of one such a consideration. 

This development of a parameter can also be approached using similarity arguments. A 
function can be generated by considering the variables which are important to the flow field and 
then constructing a non-dimensional parameter (or set of parameters) which would indicate more 
precisely the conditions leading to breakdown. This more rigorous means has been demonstrated 
by Sychev (1960) and then Hemsch (1988) as noted earlier. 

The delta wing flowfield also contains factors which are difficult to quantify. The complete 
role of the secondary vortex, which is a direct result of the viscous nature of the flow, is not well 
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understood. Both the primary and secondary port side vortex structures are illustrated in Figure 
2.1 using a titanium tetrachloride vapor method for marking vortical flows by Visser, Nelson, and 
Ng (1988). 




Figure 2.1 Visualization of the Secondary Vortex on a 70° Delta Wing 

The coherent behavior of the primary vortex in Figure 2.1 is reflected in the secondary, but with an 
opposite sense of rotation. It appears that the secondary vortex undergoes a change, visually 
resembling the breakdown phenomena of the primary vortex, well before any such occurrence in 
the primary vortex. The secondary does, however, retain a vortical flow behavior aft of this point, 
as will be detailed in the present measurements. Reported characteristics of the secondary vortex 
also include displacing the primary vortex inwards and upwards. Its behavior is also dependent on 
the Reynolds number, especially for transition to a turbulent boundary layer. Questions still 
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remain, though, as to the effect, if any, of the secondary vortices on the magnitude of the primary 
vortex pressure peaks when compared to an Euler type solution. Even more nebulous concepts, 
such as the effect of the secondary vortices on the vorticity distribution above the wing, remain to 
be investigated. 

Although the data base of delta wing flowfield information is growing, most experimenters 
are constrained to a specific configuration at a fixed angle of attack or chord location. This 
information, helpful for the sake of comparison, provides little or no information on the changing 
flow field state in the chordwise direction. Aspects such as rate of vorticity generation or 
circulation distribution in the streamwise direction are impossible to address. 

Other subtler reasons also exist which make comparison difficult. Data compiled from 
several different investigators (Kjelgaard and Sellers (1988), Verhaagen (1990), Carcaillet, Manie, 
Pagan, and Solignac (1986), Pagan and Solignac (1986)) are compared in Table 1. The sweep 
angles and the angles of attack are roughly the same. The values of vorticity presented are the 
maximum values found in the flow field and coincide with the core axis. Typically, investigators 
use the model root chord, c, and the freestream velocity, Uoo, as nondimensionalizing scales for 

the flow. The results are seen to vary quite substantially. The local semispan, s* could also be 
considered a viable length scale, since flow visualization indicates the vortex structures scale with 
the wing geometry. Use of the local semispan would account for local geometric changes due to 
sweep angle and allow for comparison of data taken at different chord stations. The data of 
Kjelgaard and Sellers (1988), for example, indicates a lower value of vorticity than that of Payne 
for a location 20% farther from the apex. If s* is used as a scaling factor, the magnitude of the 
axial vorticity component is seen to exceed that of Payne. It is also noted that identical geometries, 
such as that of Payne and the current investigation, have produced different values of the axial 
vorticity component. These differences indicate there are other factors which have not been taken 
into account and add to the difficulty in obtaining a clear picture of the physics involved. A closer 
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examination of the grid resolution of each investigation, indicated in the last column of Table 1 , 
provides some insight. The highest derived vorticity values correspond to the finest grid resolution 


Investigator 

Angie 

of 

Attack 

(°> 

Sweep 

Angle 

o 

Chord 

(mm) 

Freestream 

Velocity 

Uoo 

(m/s) 

Chordwise 

Station 

x/c 

Maximum 

Voitidty 

n 

X 

(1/S) 

n*c 

£V* 

Measurement 

Grid 

Resolution 

(y/s) 

u°° 

u<~ 

Payne 

20.0 

75 

406.4 

10 

0.5 

8,383 

341 

46 

0.04166 

Kjelgaard 

and 

Sellers 

20.5 

75 

568.8 

12.8 

0.7 

7,113 

316 

59 

0.0323 

Visser 

and 

Nelson 

20.0 

75 

406.4 

9.7 

0.5 

12,340 

517 

69 

0.030 

20.0 

75 

406.4 

10.0 

0.5 

22,774 

925 

124 

0.015 

Verhaagen 

20.4 

76 

2220.0 

' 

25 

0.5 

17,400 

1545 

193 

0.0145 

Carcaillet 

20.0 

75 

500 

20.3 

0.6 

23,061 

568 

91 

* 

20.0 

75 

500 

119 

0.6 

114,144 

480 

77 

* 

20.0 

75 

1450 

40 

0.8 

11,034 

400 

86 

* 

Pagan 

and 

Solignac 

19.3 

75 

560 

14.5 

1.4 

(wake) 

6,732 

260 

70 



* unavailable 


Table 1 Comparison of Maximum Axial Vorticity Data 

and vice versa. Since vorticity is a measure of the smallest scales of the flow, it would only make 
sense that a finer measurement grid would be able to 'capture' the high gradients of velocity 
present in the field. The behavior of the vorticity field is now examined in more detail and 
arguments for the use of this property as an indicator of the onset of breakdown are given 
Previously obtained experimental data is also presented to illustrate the potential viability of such 
arguments. 
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2.2 The Delta Wing Vorticity Field 

The aspects of the vorticity field would seem to be of utmost importance in gaining a 
better understanding of the flowfield behavior. The vorticity or rotationality of the fluid could well 
be the decisive factor in the mechanism of vortex breakdown. Both the generation of vorticity, 
which is transported into the vortex, and the convection rate downstream of this vorticity could 
play a crucial role in determining where the breakdown of the primary vortex occurs. Lee and Ho 
(1989) state that ' a stationary leading edge vortex is achieved only when the convection of 
vorticity along the core axis balances the vorticity generation from the boundary layer of the leading 
edge'. They further conclude that a reduction in the axial convection, via the adverse pressure 
gradient at the trailing edge of the planform, results in vortex breakdown and that the swirl angle 
will indicate the vorticity balance. A similar argument, postulated by Ng (1989), is that a critical 
vorticity concentration occurs, above which the aerodynamic forces cannot maintain a stable vortex 
over the airfoil. If this concentration is exceeded the vortex transitions to another state, such as 
post breakdown, to redistribute the excess vorticity. An increase in the angle of attack leads to a 
higher rate of generation of the axial vorticity component without an accompanying increase in the 
axial velocity. The subsequent increase in the vortex strength leads to vortex breakdown . 

Ashenberg (1987) attempted to model the flow about a slender wing to avoid the limitation 
of using a suction analogy, since this analogy is not capable of predicting the flow near the surface 
of the wing. He assumed the type of breakdown does not affect the aerodynamic properties of the 
wing and that conical flow is assumed near the apex which is used as an upstream condition for the 
solution. Downstream of breakdown, two dimensional sources were distributed along the vortex 
axis. Although the lift was found to be overpredicted and pressure peaks were displaced laterally 
outboard, he noted from his mathematical expressions that the expansion of the bubble caused a 
reduction of the vorticity shed from the leading edge. 
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Consider the state of the vorticity field as an indicator of the onset of breakdown. Perhaps the 
vortex structure may be unable to exist in the cohesive pre-breakdown state if, say, a maximum 
local value of vorticity in the vortex is reached or if, perhaps, an upper limit exists on the total 
amount of distributed vorticity in the vortex being convected at a given time. The first proposal 
reflects the point type property of vorticity. Determining Q x and £2y, or in polar notation, the 
radial, £2 r . and the azimuthal, D t heta> vorticity components, over various wings at a series of 
chordwise locations, and suitable nondimensionalization of the data might point to a critical 
maximum in the flowfield. The supposition of Lopez (1988) that a change in sign of Qtheta 
causes breakdown could also be verified experimentally. 

The second hypothesis arises from the reasoning that the vortex may have a limit on the 
maximum amount of vorticity per unit area or volume. If one continues to feed vorticity into the 
system, it can only 'hold' so much before it must revert to a more stable configuration in order to 
contain or transport the increase in vorticity . Extending this further, since the breakdown position 
maintains an average mean location, it could be assumed that the vorticity being generated is 
balanced by the vorticity being convected downstream for some given set of fixed conditions, such 
as sweep angle and angle of attack. If the flow conditions are then in some way altered, so as to 
add more vorticity upstream of the breakdown without a corresponding increase in the convection 
rate downstream, the breakdown would be seen to move upstream. This would indicate that some 
sort of critical condition, based on a maximum vorticity distribution, exists causing the initial 
vortex structure to transition to the post breakdown state. An increase in the angle of attack or a 
decrease in sweep angle would momentarily cause the relative vorticity generation rate to become 
higher than the convective rate. Hence, the critical vorticity distribution would be reached earlier 
(i.e. upstream of the initial breakdown location) and the breakdown would move towards the apex 
until a stable situation is again reached. For this reason it may be advantageous to consider the 
relation of the vorticity generation to convection terms. 
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The argument of a critical vorticity distribution can be substantiated in light of the work by 
Pagan and Solignac (1987) mentioned previously in 1.2.2. Their results indicate that a maximum 
vortex strength, as given by the maximum swirl velocity ratio, is strongly dependent on the local 
freestream pressure gradient, controlled by moveable flaps in their tunnel section. This can also be 
interpreted as the maximum amount of vorticity at a given station, or circulation/vortex strength, is 
limited by the ability of the flow to move downstream, regulated by the pressure gradient. Thus, 
the state of the vortex, with regards to its breakdown potential, can be described by a ratio of 
strength generation conditions to flow transport conditions. Specific ratios could include flowfield 
properties such as the circulation at a station to convective velocity, or geometrical variables, such 
as a function of a and A to imposed tunnel gradient. 


Empirical verification of such proposals can best be investigated by utilizing parameters 
based on the state of the flow as revealed by the ratios of the relevant variables. Several examples 
of this were outlined in Chapter 1. Interestingly, if the swirl parameter 


S = 


2Q 


( 1 . 22 ) 


mentioned previously in the use of tornado studies and reduced to the form of V0/V x for delta 

wings, is examined further, several possibilities are revealed which may be useful in correlating 
delta wing flows. Rewriting the parameter in terms of the average axial velocity component, V x 


S 


roT 


2nr o 2v x avg 


2nr o v x avg 


(2.1) 


the relation can also be seen to indicate the transport of circulation in the axial direction. The 
parameter can also be expressed in terms of the average vorticity over the area normal to that 
bounded by r Q 
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S 


r o ^x 
2V Xavg 


( 2 . 2 ) 


This is simply the inverse of the Rossby number and can be interpreted as a measure of the 
convection of the vorticity. If the radius r 0 is defined to be that of the subcore of the delta wing 

vortex, which will be shown later to contain the majority of positive axial vorticity, the streamwise 
behavior of this parameter could provide an indication of whether or not breakdown was imminent. 

The flowfield can be examined on a point by point basis, using the local axial values of 
vorticity and velocity and some length scale, such as the chord or the local span, s*, 


to observe any type of local effect. Since a comparison of the data using a point property analysis 
is susceptible to effects such as the grid resolution, an integrated approach can also be used. The 
radius of the subcore could be the length scale, 

S = (2.4) 

v x 

and averaging the axial vorticity and velocity components over the region of the subcore could help 
in alleviating potential resolution discrepancies. 


As can be seen, different variables can be incorporated into these nondimensional 
parameters. If the velocity in the numerator above was V 0 instead of V x , then perhaps this would 

provide the 'critical condition parameter 1 . This may be so. This type of guessing can continue ad 



infinitum and therefore the parameters in this study are limited to those based on strong physical 
arguments. 
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The suppositions presented above were based on a examination of some existing data and 
theories in the literature. This data and the conclusions drawn from its examination are now 
presented below as a basis for the discussion of the objectives of the present study. 


Axial vorticity was calculated from LDV velocity data taken by Payne (1987) and Anders 
(1982) is shown in Figure 2.2. The axial vorticity, based on polar coordinates is given below: 


, = i »rv e ) 1 av r 

x r dr ~ r d0 


(2.5) 


Under the assumption that the azimuthal gradient of the radial velocity is negligible, axisymmetric 
flow and a negligible radial velocity , the relation becomes 


£2x 



dVp _ w dw 
dr “ r + dr 


(2.6) 


This is applied to data acquired along a traverse directly through the vortex core. The angles of 
attack of Payne and Anders were 20° and 19.3° respectively and the corresponding sweep angles 
were 70° and 68.2°. No breakdown is occurring over the wing. Both core traversals were taken 
near mid-chord (x/c = 0.5 and 0.67). The Reynolds numbers are comparable (425,000 and 
170,000) for this type of flow. Payne has a chord length of 406.4 mm while Anders has 137.5 
mm. If the spanwise variation in the vorticity (calculated using central differencing in polar 
coordinates) is nondimen sionalized using the local span length and the freestream velocity, the 
results are comparable. 
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Non-Dimensionaiized radius 

Figure 2.2 Axial Vorticity Data of Payne and Anders 

Axial vorticity was also calculated from LDA velocity data acquired by Iwanski (1988) over 
a one inch thick, 70° sweep, flat plate delta at a = 30°. The vorticity can be nondimension alized in 
two ways: by the local half span, s*, to account for the increase in the vortex length scale, or by 
the root chord, as is generally done in the literature. Figure 2.3 illustrates how the velocity data 
from Iwanski can be used to observe the spanwise variation in vorticity at different chord stations. 
The axial vorticity profiles are seen to increase in magnitude and narrow in width in the 
downstream direction up to x/c = .41 1 (Figure 2.3a). Further downstream. Figure 2.3b, the 
magnitude of the values drop off and the peaks broaden as the breakdown region is traversed. 



Non-Dimensionalized Vorticity , 
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From the above discussion, it is evident that the vorticity field undergoes large changes in 
the axial vorticity distribution as the breakdown region is encountered. Thus, the measurement and 
subsequent analysis of vorticity was regarded as an important aspect of this investigation. Before 
the objectives and goals of this study are laid out, however, a hypothesis is presented in the next 
section which deals with another aspect of the flowfield that could indicate the state of the vortex 
relative to breakdown, in terms of the tangential velocity component and the pressure field. 


2.3 A Heuristic Proposal 

A proposal is outlined below which would allow the state of the delta wing vortex, relative 
to its condition at breakdown, to be evaluated in terms of measurable flow quantities at a given 
chordwise station. The region referred to as the viscous subcore in the vortex is suggested to act 
as a solid body downstream of the apex generation region, where the wing geometry is thin 
compared to the local flow. This being so, the cross sectional flow can be treated as if the radial 
pressure gradient supplies the required centripetal acceleration to maintain a pre or post breakdown 
vortex. 

The flow is assumed to follow a steady, axisymmetric, incompressible behavior, as was proposed 
by Hall (1961), with the appropriate equations (1) repeated here 


Momentum: 



+ V 


dVr 

x0x 


V 2 






r: 


r 


(2.7a) 
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(2.7b) 


(2.7c) 


Continuity: ^ (rV r )+ (V x ) = 0 (2.7d) 

The analysis now deviates form that of Stewartson and Hall (1963) in the core region as 
described earlier by introducing a limiting assumption. Restricting attention to the viscous subcore 
region, a further simplification is introduced based on observations in the wind tunnel using flow 
visualization and measurements performed with x-wires and 5 hole probes. This simplification is 
also the basis for many solid body theory models, that being: 

The radial velocity, V r , is negligible in the viscous subcore upstream of breakdown. 

Although this assumption is incorporated by many of the studies mentioned earlier, further 
substantiation of this premise is now given based on the observations of this study and a previous 
study by Visser (1988). The negligible radial velocity appears to be confirmed experimentally 
using flow visualization. No tracer particles are seen to spiral into the subcore region when 
released either upstream of the wind tunnel test section or locally about the model, excluding the 
apex. Interestingly, Davies-Jones (1982) reports that no entrainment into the core region is seen 
for tornado structures in the flowfield away from ground level. It also appears that tracers placed 
into the subcore region remain in the 'core' and are neither transported to the outer vortex flow 
region nor ’spun out 1 to some particular radius within the subcore region. Laser light sections 
have shown a tracer filled cross section for the entire chordwise length of the vortex as reported by 
Visser, Nelson, Ng (1988). Any radial velocity gradients would ‘spin out’ the particles and cause 
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b) 

Figure 2.4 Vortex Visualization a) Sub - Core b) Laser Light Sheet 
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a dark center hole to appear at some point down the length of the vortex, creating in essence a ring 
of smoke at the outer edge of the subcore. No such phenomena is observed. Measurements 
performed using an x-wire and a five hole probe indicate that this core, based on the distance 
between the maximum tangential velocity components, remains approximately constant in diameter 
upstream of breakdown. Flow visualization by Visser, Nelson, Ng (1988) and Payne (1987), in 
Figure 2.4a and b respectively, indicates a core size that appears to be of a cylindrical nature rather 
than conical. It is quite obvious that near the breakdown region, as the core stagnates and begins 
to expand, radial velocity components can not be neglected. However if this analysis is restricted 
to the pre-breakdown state, neglecting V r is not an unreasonable approximation. 

Introducing this radial velocity approximation further simplifies the above equations to the form: 


V _ I3P 

r ‘ p dr 


(2.8a) 



(2.8b) 


(2.8c) 


(2.8d) 


and substituting the resulting continuity condition into the x direction momentum yields 


1 dP _ y _d_ ( *V X >| 
p dx " r dr dr J 


(2.9) 
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Recall that these equations are being used for the region of flowfield up to the outer boundary of 
the subcore or the inside of what is referred to as the inviscid Euler rotational region. 

Several aspects of the flowfield become immediately apparent. The change in the axial 
velocity component, 8V x /9x, is zero from continuity. Leibovich (1983) stated that the axial 

velocity continues to increase along the axis and Verhaagen and van Ransbeeck (1990) have 
measured the axial velocity and found it to increase with axial direction up to a certain downstream 
distance. The theory of Stewartson and Hall (1963) also indicate an axial component of velocity 
that increases with axial distance, again up to a point, whereupon it becomes constant as shown by 
Verhaagen and van Ransbeeck . Stewartson and Hall (1963) include a radial velocity term in their 
model, directed inwards to accommodate the increase in axial velocity. This radial term depends 
on distance and as x increases, the magnitude of the radial velocity component decreases to a 
negligible amount 

The constant axial velocity component that results from the present derivation is because of 
the zero radial velocity constraint imposed. Yet both from the theory and measurements noted 
above, the axial velocity seems to reach a constant value at some distance downstream of the apex. 
It is known that the axial velocity stagnates at breakdown and even exhibits reversed flow. Hence 
one could presume at the very least that the axial velocity would reach a maximum and then begin 
to decrease to the breakdown point. The deceleration of the flow occurs over a very short distance, 
on the older of 2 or 3 core diameters. Thus it may be that for a distance directly upstream of the 
breakdown zone, the axial velocity reaches and maintains a maximum value, similar to that shown 
by the above investigators. This being so, a zero radial velocity would be reasonable and the 
above assumptions justified on more than just the basis of flow visualization. 

It should also be noted from the x-momentum equation that the axial pressure gradient is 


not zero. Since 



av x 

~5T 
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< 0 


( 2 . 10 ) 


due to the jet like nature of the axial flow, the pressure gradient in the axial direction is also 
required to be negative for av x A)x = 0. This pressure gradient has been verified experimentally 

by Verhaagen and van Ransbeeck (1990) and by Lam bourne and Bryer (1961). Present 

measurements will be discussed later. As mentioned earlier, Stewartson and Hall ( 1963) found 
that dV x /dr = 0 for a low enough value of the kinematic viscosity (see Batchelor (1967)), which 

would suggest a constant pressure along and across the core. 

The zero radial velocity also implies a solid body type of rotational flow. As noted earlier, 

flow visualization has indicated a viscous subcore diameter that does not appear to vary in a conical 
fashion, but rather maintains a cylindrical form. This implies V 0 to be a constant at a specified 

radius for any axial direction, that is no dependence on x. Using this with the the 0 momentum 
equation gives 


a 

dr 




and a solution of V 0 = cor which agrees with the solid body implication. 


(2.11) 


This is not to say that this solid body rotation behavior actually occurs. Indeed there may 
be shear occurring in the fluid at different radii from the core center. If the rotational rates of the 
subcore region is approximated by the above 



(2.12) 


with V 0 typically on the order of 1 .3 - 1.5 and d = 6-10 E-03 m, to values up to 5000 rad/s 
or almost 800 rps (48,000 rpm) for the tests presented later in this study. Higher tunnel speed 
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would increase this rate substantially. Data from Payne (1987) based on LDV data gives rotational 
speeds approaching 60,000 rpm. 

Nonetheless, the particles appear to maintain their respective distances from the core, 
indicative of a solid body rotating flow. In addition measurements of the velocity field through the 
core indicate a linear variation with radius, also a property of a field undergoing solid body 
rotation. The total pressure drop, which has been measured at the core center, is seen to 
accompany such flows and so for the moment this assumption will still be maintained. 

The radial momentum equation provides the most interesting possibilities as to insight on 
the physics of the flowfield. It is interesting to note that the reduced radial momentum equation 
derived above by assuming a zero radial velocity is identical to that presented by Stewartson and 
Hall (1963) who had imposed a slenderness condition. 

XsL = IdP (2.13) 

r p dr 

This relation can be interpreted as a balance between the radial pressure gradient and the 
centripetal forces. A better explanation is that the radial pressure gradient is such to exert a large 
enough centripetal force to keep the vortex together. Alternatively it can be said that for a given 
pressure gradient, there is a maximum centripetal acceleration that can occur and hence a maximum 
tangential velocity. This type of flow is also referred to as 'cyclostrophic flow (Davies-Jones, 
1973). 


It could be surmised from the above that if an initially balanced system was to be upset 
such that the radial pressure gradient became less than the required acceleration (i.e. a local 
pressure rise along the axis), the particles would no longer be held to their orbital paths and would 
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attempt to travel in straight lines, i.e.. the core would expand or diverge. Conversely, if V 0 was 

reduced it could be expected that the pressure gradient would drive the core to a smaller radius with 
a subsequent radial inflow opposite to the outflow seen in the breakdown region. Or, extending 
this concept, if the centripetal acceleration (tangential velocity) was modified to remain smaller in 
magnitude than the local radial pressure gradient, breakdown could possibly be delayed. 

This argument is consistent with experimental observations of the swirl angle parameter 
defined as: 


swirl angle = tan 



(2.14) 


where V x is the axial velocity component and V 0 the tangential component. Typically the swirl 

angle increases in a chordwise direction to a value of between 40° to 50° whereupon breakdown of 
the vortex is observed. Larger magnitudes of V 0 relative to V x would seem to indicate breakdown 

was imminent. Since the increase of the acceleration term in the radial momentum equation is 
dependent on the square of this tangential velocity, V 0 , both of the above arguments point to a 

reduction of the local tangential velocity as a means to delay the breakdown process. A further 
implication is that the circulation or strength of the vortex, defined as the line integral in a plane 
normal to the vortex axis: 


T = jV«dl 


( 2 . 15 ) 


would have to be reduced. That is, the vortex would have to be weakened in order to delay 
breakdown and possibly enough weakening could move or maintain the breakdown region in the 
wake. 
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It would be helpful if equation (2.13) could be used in direct conjunction with properties 
that are measurable in the core region. Equation (2.13) could then be utilized to indicate the current 
status of a given vortex with respect to breakdown. Two approaches can be taken to using this 
relation. First, equation (2.13) can be simply interpreted as: 'What pressure gradient is required to 
maintain the maximum centripetal acceleration?'. Thus the relation can be evaluated at any point in 

the vortex using experimental data and without resorting to integration or any assumption on the 
behavior of V 0 . The most interesting spatial locations would occur where the tangential velocity is 

the largest Thus equation (2.13) can be stated as 


V ^ 
dm ax 

^dmax 


@P\ 

l ar J@r y 

dmax 


(2.16) 


The value of V 0max 


is is generally 


taken to be the edge of the 


subcore region. The pressure data at 


this location must also be available to obtain the local gradient The local ratio of the pressure 
gradient term to the acceleration could then then be compared. 


As an alternative approximation, equation (2.13) can be integrated from the vortex center, 
denoted as 1 below, to the radius of the subcore region, 2, using a V 0 = tor, solid body 

distribution 




yielding: 




(2.17) 


Denoting Pj as the pressure at the core axis and setting r, = 0 gives: 
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P p p to 2 r 2 2 _ p V 02 2 

2 r axis - 2 ~ 2 

or 

^ " P axis^ 

= 1 (2.18) 

A further implication is that the pressure gradient required is dependent on the square of the local 
swirl ratio at r 2 



^ (^2 ' p axis) 

O O = 1 (2 1<; 

pU 2 tan(swirl 2 )2 

What remains is to define where V 02 is to be measured, that is, where the boundary 2 is located. 
If the value V 02 is taken as the largest tangential velocity in the flow field, as in approach 1 , the 

value of P 2 is also required at this point. Four cases can be examined, the first three of which are 
based on the ratio 

2 (P 2 - P ax j s ) _ 2 (AP) 

P v 0 2 2 = PV 02 2 

and are listed as follows: 

i) Use the static pressure difference, AP = P^ - Paxiw noting that 

P ~ = P ~to,al ‘ PU “ 2 a " d P “ is - “ ^ - taxis' 2 

ii) Use the total pressures, AP = P f - P, 

ioo t axis 

iii) Use the pressure at r V0max> AP = P V02 - P^. This requires all the velocity 
components at to be known at ly to obtain the static pressure from the total pressure. 
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iv) Use two values of pressure to either side of r V02 . These values can be used in conjunction 

with the ratio above, but also with the approach to equation (2.16). This is the most correct 
evaluation, as the closest approximation to the local pressure gradient is obtained. 


The question now arises as to what this represents? Can this parameter be simply be 
interpreted as the ratio of the radial pressure to the centripetal acceleration of the fluid? That is, m 
order to maintain the centripetal acceleration necessary to keep the fluid on the circular trajectory 
defined by the vortex, a certain pressure gradient must be required? If the ratio is > 1, is the 
pressure gradient is large enough to maintain the coherent structure? If it is < 1 , are the required 
centripetal forces are too large and the vortex diverges? Could this possibly be used as some kind 
of breakdown criterion? Available data and the present tests will be used to examine the above 

conjectures in Chapter 5. 


Although the above derivations did not assume an inviscid behavior, viscous effects do not 
appear explicitly in the the reduced radial momentum equation, which formed the basis for the 
above arguments. The effect of viscosity should not be overlooked, however, as its effects are 
seen quite readily if the entire flow field above a delta wing is considered. The existence of the 
secondary vortex is a direct result of viscous interaction at the wing surface. Consider, for a 
moment the possible effects of this secondary vortex using the idealized situation below: 





Figure 2.5 Ideal Delta Wing Flow Vortex Representation 
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It would be expected that the secondary vortex induces a velocity field on the primary vortex. This 
would tend to cause the value of the tangential velocity, V0 , in the primary to be larger in 

magnitude on the side nearer to the secondary than on that closer to the wing centerline. This is 
observed experimentally. Thus it may be that it is the effect of the secondary which causes V0 to 

increase to the point where 


YfL > IdP 

r p dr 


and thereby initiate breakdown. 

The rotation of the secondary vortex would also cause the distance between the core of the primary 
vortex and the wing to be greater than if the secondary were not there. This would have the effect 
of reducing the measured suction peak value at the surface of the wing. The implications of this 
concept could be far reaching if it could be substantiated. Reduction or removal of the secondary 
might not only delay breakdown, but also increase the lift by moving the primary vortex nearer to 
the wing surface, possibly to the maximums calculated by current Euler codes. Earlier preliminary 
tests using small angular tabs as vortex generators were conducted by the author to investigate the 
possibility of altering the breakdown location by manipulating the flow near the surface. Tabs 
were placed near the apex, both on the upper surface and the bevelled leading edge. The 
breakdown location was found to move further aft for certain tab orientations. The effect these 
generators had quantitatively on the actual vorticity distribution of the secondary, or even the 
primary vortex, was not determined, however. 
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2.4 Goals and Objectives of Present Study 

The above research summary and that of the previous chapter can be summarized as three 
important conclusions, pointing the way for the study described hereafter and for further 
investigations. First, it is quite apparent that the external conditions which can be imposed on the 
vortex field to initiate or delay breakdown are well known and documented. These include the 
pressure field, the local swirl of the flow, based on either circulation or vorticity, and the axial 
velocity field, regardless if the vortex is confined or free. Additionally, from an experimental point 
of view, variation of a particular parameter has shown that the onset of breakdown is not caused 
foremost by a change in the swirl of the local flow, which subsequently leads to stagnation of the 
axial velocity and on through the breakdown process. Breakdown can just as easily be caused by 
an external change in the adverse pressure gradient which may then lead to a stagnation and then a 
change in swirl followed by core expansion and breakdown. Depending on the flow conditions, a 
combination of these factors may be required. 

Secondly, flowfield measurements taken to verify theoretical postulations are needed. 

Much research has been undertaken to determine aerodynamic behavior and verify existing codes 
for delta wing flows, but data acquired for the purpose of examining the flow state as it approaches 
breakdown in order to observe any significant changes is not extensive. In general, flowfield 
surveys, using pressure or LDV techniques are expensive and time consuming and thus relatively 
sparse. 


Finally, it is very important to accurately document the external conditions imposed on the 
experiment when obtaining a set of measurements. This is naturally assumed to be done by 
investigators and is most often the case. However other aspects, generally not accounted for, may 
exert appreciable influence on the results. One of the these aspects is grid resolution. This is, 
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naturally, a serious point of contention for numerical computations of a flow domain, but how 
often is it addressed by experimenters when conducting a series of flowfield measurements? 
Measurements conducted in flows with strong gradients, as is certainly the case here, are 
dependent on grid resolution as much as the associated numerical simulations. 

Based on these conclusions, a direction for the present study was derived. The lack of data 
on the chordwise progression of a vortex above a delta wing prompted a systematic investigation to 
measure the delta wing flowfield at a set of various stations. The flowfield properties including 
circulation and vorticity were felt to be of interest as it has been demonstrated that an increase the 
circulation will initiate breakdown. Angle of attack and geometrical variation, including sweep and 
chord length were varied in an attempt to determine if the flow reaches some measurable critical 
state, particularly involving aspects of the vorticity components, that would then initiate a flow 
transition to the post breakdown state. 

Specifically, the variations in flowfield characteristics were derived from velocity field data 
acquired at a series of chordwise stations using a hot wire technique presented in the following 
chapter. Acquisition of the data was performed using different grid sizes to examine questions on 
resolution. Single wire spectral data was obtained in the core upstream of breakdown to observe 
any dominating frequencies for interests sake. Additional measurements were taken in the form of 
surface pressures and total pressures to try and correlate these quantities with the flowfield 
velocities. Flow visualization, including on and off surface, was used to aid in interpretation of 
the flow. 

An indication of the state of the vortex based on these measured properties was the basis 
for examining the data. If a measured or derived property did not indicate changes suggesting 
breakdown was forthcoming or had occurred, then it could be regarded as a minimal factor in 
terms of a parameter analysis. Initially, each flowfield station could be analyzed in terms of their 
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local and overall properties, including derivative and integrated characteristics. Comparison with 
theoretical criterion presented previously could then be made. Scaling the data in various ways 
would hopefully enable better comparisons at a local scale and/or a global scale of the flow 
properties at each measured station, especially with respect to other data in the literature. Finally, 
the flow can be examined with the aid of the relevant parameters described previously to attempt to 
quantify the flow state at each chordwise station. 



CHAPTER IH 


CROSSED HOT WIRE ANEMOMETRY 


3.1 A Brief Overview of Current Techniques 

Many researchers have used two hot wire anemometers in an X-wire configuration to 
determine velocity fields and an extensive list of references can be found in Freymuth (1982). The 
measurement of X-wire voltages and their conversion to velocity can be divided into two 
categories: the table lookup method and the effective velocity method. 

The table lookup method requires that the probe be rotated through a series of angles and 
the velocity varied at each position. In this way a table of wire voltages versus angle and velocity 
is generated. Curves fitted through these points can be generated in such a way that, given a pair 
of measured voltages in an unknown flowfield, a velocity and angle can be found and the 
subsequent velocities along specific axis can be determined. Lueptow, Breuer, and Haritonidis 
(1988) give a description of this technique. This can be used for both X-wire and triple wire 
probes. The major drawback is that it is time consuming , especially if temperature compensation 
is required, and a rotatable apparatus is required. The benefits, however, include no assumptions 
on the probe geometry or variable cooling rate estimations for velocity components not normal to 
the wire. 

The effective velocity approach is the more widely used method. This technique centers on 
the concept that the wire senses an effective velocity comprised of velocity vectors normal and 
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tangential to the wire and thus cool it at different rates. J0rgensen (1971) has expressed the most 
general form of this equation as: 


U eff 2 = U N 2 + k l 2U T 2 + k 2 2lJ B 2 (3T) 

for an X-wire lying in the NT plane as depicted in Figure 3.1. 



The k values are the yaw and pitch factors which are functions of the wire and even the 
yaw and pitch angles to an extent Typical values taken are kj=.2 to .3 and k 2 = 1.08 to 1.12 

although a complete study by Jprgensen (1971) demonstrates this variation. Expressions for the 
velocity components are then derived based on the probe geometry and/or the particular flowfield 
geometry. 

To illustrate the extension of this measurement technique to an X-wire configuration, 
consider the geometry of two wires in Figure 3.2. 
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u 


Figure 3.2 Cross- wire Geometry 

An unknown velocity vector V, comprised of u and v (w = 0) and at some angle a, is to be 
measured by wires 2 and 1 at angles (3 + a and 90° to (3 + a respectively as illustrated in Figure 
3.3a It is assumed that only the normal components of the velocity on each wire are required. 

This is equivalent to a reduction of the Jprgensen equations to the form 

U eff = U N ( 3 - 2 ) 

The measured voltages are Vj and V 2 and are assumed to act at 90° to the direction of the wire. 
That is, the velocities seen by wires 1 and 2 to are derived from V j and V 2 . 



Figure 3.3a Unknown Velocity Vector with u > v. 
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Thus, from geometry, 


cos (a + P) = cosacosp - sinasinP = y yr - y = yi ( uV 1 ' vV 2) 


or solving for u, 


u 


cos (a + P) + vV 2 

vl 


(3.3a) 


Similarly, from sin (a + P), 


u 


sin (a + P) + vVj 


(3.3b) 


Equating these two expressions for u and solving for v yields. 


v = V jsin (a + P) - V 2 COS (a + P) (3.4) 


If, from geometry, a + P = 45° then 

v = X (Vl ' V 2 > 

(3.5a) 

and subsequently 

u = ^ (V 1 +V 2 ) 

(3.5b) 


Hence the u velocity component would be proportional to the sum of the measured voltages 
(velocities), while the value of v would be proportional to the difference. 
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Now, consider Figure 3.3b below where the value of v > u. 



Figure 3.3b Unknown Velocity Vector with v > u. 


In order to resolve the velocity into the known wire angles, the expression 


cos (a - p) = cosacosp + sinasinP 


is required. The resulting expressions for u and v become 


u = x ( V ! - V 2 ) 


v = IT (Vi + V 2 ) 


(3.6a) , (3.6b) 


which are the inverse of the initial derivations. 


It is readily seen that this type of derivation is useful only if the general flow direction is 
known in order that the probe may be oriented correcdy. In addition, the w (out of plane) 
component must be negligible and the effects of tangential cooling on the wire are ignored. Many 
investigators choose to use this approach, because general inclusion of all the terms proposed by 
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J0rgensen (1971) leads to a series of equations with multiple solutions which are difficult to solve 
uniquely. Flowfields are simplified or terms are ignored to accommodate simpler forms of the 
Jorgensen equations. Further derivations illustrating the increased complexity of inclusion of these 
terms in a general form can be found in Appendix D. 

The inherent complexity of these equations is evident. It appears that no one has used an 
X-wire to measure an unknown constant velocity flowfield, that being three magnitudes and three 
directions, even with successive rotations. The measurements that are taken, are performed in a 
flow where either the flow directions are known (Lofdahl (1986), Pailhas and Cousteix (1990)), or 
where there is a strongly preferred direction of mean motion in the three-dimensional shear flow 
(Mojola (1974)). Lofdahl (1986) fails to address the problem by disregarding the tangential 
cooling component effect. Browne, Antonia, and Chua (1989) in their discussion of calibration 
methods for yaw response in x-wire probes state that as with angle methods, the effect of the 
velocity component normal to the wire must be neglected. Klatt (1969) and Andreas (1978) also 
take this approach. 

Other hot-wire methods have been employed in an effort to measure the flow field and 
include configurations using up to nine wires. Janjua, McLaughlin, Jackson and Lilley (1982) 
have used a six-orientation single wire method in the axisymmetric flow of a gas turbine engine 
combustor to determine an unknown velocity field. Comparisons of the results with independent 
data have demonstrated the reliability of this method. The largest uncertainties were found to exist 
in the measurement of the turbulent shear stresses. Further tests by Jackson and Lilley (1986) 
indicate that this technique adequately measures the properties of the flowfield independent of the 
dominant flow direction except if it is aligned with the probe axis. In addition the time-mean 
velocity in the probe direction is inadequately deduced. It was also reported by He (1988) that this 
technique has the same precision as a muti-hole Pitot tube with the advantages of hot wire 


anemometry. 
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Triple wire probes are available commercially and contain a third wire which is oriented 
such to provide a means for obtaining the entire velocity vector in a single measurement. Both 
DISA and TSI manufacture triple wire probes capable of measuring all three velocity components 
and their respective directions. There is a constraint on the use of these probes, however. The 
velocity velocity being measured must lie in an acceptance cone of 70° about the axis of the probe 
as shown below in Figure 3.4 Thus the flow direction must be approximately known, as noted by 
Gaulier (1977), and depending on the derivation the resulting equations become fourth order 
transcendental, requiring interpretation of the roots for validity, as detailed by Lekakis, Adrian, and 
Jones (1989). 



Figure 3.4 Triple Wire Acceptance Cone 
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This restriction on the direction of the velocity vector, which again implies a known 
velocity direction, is too severe for the measurements required in the present study. Consider the 
velocity vector shown in Figure 3.4b, consisting of the components u, v, w, and scaled by the 
ffeestream velocity Uoo. If, for the moment, w is set to zero and the angle between the axial 

component of the velocity and that of the lateral velocity component is defined as 

9 V = tan'H~) (3-7) 

For an axial velocity of u/U«> = 1 and the maximum allowable value of 6 V constrained to 35°, 

v /Uoo can not be greater than 0.7. As illustrated in Figure 3.5, should w/Uoo attain a value 
max 

greater than zero, as is certainly the case in the present flowfield, the velocity would lie outside of 
the cone of acceptance. 



Figure 3.5 Constraint circle 


Since the angle v max /Uoo makes with u must be < 35°, and for any point on the circle in Figure 


3.5 above 


o 9 1 

v ^ = V ^ + W 

v max v 


90 


it can be seen that 


0 = tan-ir-Xmax.^ = tan -l 

v max l u ) 


V 


2 2 
v ^ + w z - 


u 


< 35° 


(3.8) 


Hence, a value of u = 1.3 requires v ^ + w ^ = 0.83. Seven hole probe data from a core traverse 
by Payne (1987) in Figure 3.6 indicates regions of the flow, such as u = 1.3, w = 1 , and v = 0.25, 
which violates this constraint Thus the use of commercial triple wire probes was deemed 
unacceptable for the present study using conventional data processing procedures. 



Figure 3.6 Seven Hole Probe Traverse Data 

This triple wire constraint has been overcome, as reported by Jacobsen (1977), for flight 
data acquired in the tip vortex wake of a T-33 aircraft. A three wire probe was mounted on the 
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nose boom of a Gates Leaijet such that the probe axis was parallel to the longitudinal axis of the 
aircraft The axial component of the jet is thus added to the tip vortex axial component and 
therefore keep the velocity vector in the cone of acceptance. 

The use of hot wire probes with more than three wires to measure instantaneous vorticity 
values has been proposed by Kovasznay in the early 1950’s. The Kovasznay type probe. Figure 
3.7a, consists of four prongs supporting four wires to form a Wheatstone bridge, when operated 
by a constant current anemometer. 



Figure 3.7 Kovasznay Type Four Wire Probes 

Kastrinakis, Eckelmann, and Willmarth (1979) found it was not possible to measure 
instantaneous vorticity values or rms values with a Kovasznay type probe in flows with large 
cross-stream velocity components. They suggested that each wire be supported by its own set ot 
prongs as in Figure 3.7b. Vukoslavevic and Wallace (1981) have built and tested this type of 
probe, but tests in a low speed boundary layer indicated errors occurred if the cross-stream velocity 
components were accounted for as they themselves could be in error by as much as 80%. They 
concluded that this four wire probe does not provide enough information to determine the 
instantaneous streamwise vorticity or cross-stream velocities. A review by Foss and Wallace 
(1988) discusses efforts and ensuing complications based on other configurational four wire 
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probes, as well more complicated five and nine wire probes, to obtain instantaneous values of 
vorticity. 

Swirl angles larger than 45° exist in the present flowfield which eliminates the use of 
commercially available triple wire probes. Due to the amount of data required to examine useful 
trends in the flow, the six-orientation technique was eliminated in deference to a procedure 
involving multiple x-wires. In addition, because v, and w are on the order of u, the flow will exert 
appreciable tangential cooling and cooling by velocity components normal to the plane of the X- 
wire. It would thus be incorrect to disregard these terms of J0rgensen's equation. The nature of 
the flow under investigation precludes any assumptions concerning flow direction or relative 
magnitudes, apart from the fact that the axial velocity maintains a streamwise sense for the flow 
upstream of the breakdown region. Due to the inherent non-linearity of the associated equations 
and the resulting multiple solutions, it became apparent that several measurements at each spatial 
location are required using different x-wire configurations. A summary of the technique used in 
this study follows. 


3.2 A Method for Unknown Three Dimensional Flows 

Minimization of the number of spatial measurements required to obtain the mean velocity 
field was felt to be essential to providing the most accurate data possible. For this reason, the table 
lookup method, detailed by Lueptow, Breuer, and Haritonidis (1988), was initially examined. The 
required probe uses a sensor plane parallel to the probe axis, such as the DISA 55P61 in Figure 
3.8b. 
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Figure 3.8 Gross-wire Probe Configuration 

Calibration involved positioning the probe axis parallel to the freestream and varying the velocity 
while recording both wire voltages. The probe was then rotated to keep the plane of the wires 
parallel to the oncoming stream and the procedure repeated. In such a way, a calibration relation of 
the wire voltage pairs as a function of the speed and position was created. By suitable inversion, 
the velocity and direction could then be determined from a measured pair of voltages in the 
unknown field. 

Although the normal and tangential effects were thus included, this method fails to account 
for velocities normal to the plane of the wire. Thus two grid sweeps of the field were required 
with the probe rotated 90° for the second sweep. In this way it was hoped that u and v would be 
measured on the first pass and u and w on the second. The results of a u - v pass for the vortex on 
the right side of the delta wing can be compared to data acquired by Payne (1987) in Figure 3.9. 
The survey plane was taken normal to the planform surface and normal to the x direction. Each 
axis has been scaled by the local semispan. Thus a y/s =1.0 corresponds to the leading edge of 
the planform. 
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0.0 y/ s 1.2 


(a) 




Figure 3.9 Table Lookup Method Results for sweep = 70° alpha = 20° x/c = 0.5 
a) u/U M b) v/U M clu/U^, Payne d) v/U^, Payne 


The velocity component normal to the plane of the wires was found to contaminate the 
measured velocities. The axial component reflects the influence of the w component, that being 
normal to the wing surface, in the presence of two lobes to the left and right of the vortex center. 
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The presence of the normal component in these regions causes the measured axial velocity to 
appear larger than it actually is. The effect on the transverse component is much less pronounced 
and can be seen from the displacement of the zero velocity contour. In a similar manner, the axial 
velocity obtained in the second sweep was contaminated with the transverse, v, component. 


Several schemes were employed to isolate the unwanted effects including iterative 
schemes. No procedure was completely successful as the resulting field would contain unresolved 
regions. As a final resort, the use of a physical shield on the probe during acquisition of the data 
was investigated. It was felt that this would possibly eliminate the unwanted effects of the velocity 
component normal to the plane of the wires. Although this concept has not been utilized on x-wire 
configurations, it has been employed in conjunction with a single wire. Giinkel, Patel, and Weber 
(1971) affixed a disk-like shield around a single wire to minimize the effects of lateral velocity 
components in measuring reversing flows from 0.3 - 10 m/s. Their conclusions indicate drastic 
improvements. A shield of the type shown in Figure 3.10 was therefore constructed. The topside 
is removed for clarity. 



Figure 3.10 Cross-wire Probe Shield 
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Essentially, the concept behind the shield is that it would allow only the velocity components in the 
plane of the wires to be detected. After several iterative designs, a configuration was reached that 
gave a qualitative comparison with the data of Payne, as is shown in Figure 3.11. 



Figure 3.11 Shielded Probe Method Results for sweep = 70° alpha = 20° x/c = 0.5 

a) u/Uoo b) v/Uqo 

Quantitatively, however, the shielded probe data still reflected an increase in the velocity 
components and the influence of such a device on the local flow is questionable. This line of 
investigation was therefore terminated in favor of a method utilizing the effective velocity concept. 
This is not to say, however, that this indicates a shielded probe technique to be invalid. 
Optimization of the shield design may yield data that represents the two desired components of the 
flow being measured. 

Since the component of velocity normal to the plane of the x-wire cannot be neglected in 
parts of the field, it would seem mandatory that the Jprgensen equation incorporating all the 




97 


cooling terms be employed. With two different configurational measurements, it is possible to 
obtain a series of four non-linear equations in four unknowns. Thus, it could be surmised that a 
solution exists, possibly in an explicit form. Depending on the configuration used, however, the 
equations will only yield information on some of the quantities, such as u,v,w with unknown 
directions. All the attempts so far at solving these equations without some type of knowledge of 
the field, such as directions or relative magnitudes of one velocity component much larger than 
another, or simplification, such as ignoring the effect of tangential cooling, have led to solutions 
which possess regions in the measured field where singularities were seen to occur. 

The method adopted for the present study is outlined by Sherif and Pletcher (1987). It is 
based on the effective velocity measurement method, and uses an X-wire array with the wires in a 
plane perpendicular to the probe axis. This configuration is depicted in Figure 3.8a. Their 
procedure requires two rotations of the x-wire and the resulting equations are solved explicitly in 
terms of u, v, and w. In addition, they state that the flow being measured should be of a three 
dimensional nature. Initial tests indicated it was possible to determine magnitudes of the velocity 
comparable with that obtained by other investigations to ui 2.3%. Topographically, this 
procedure also seemed to give the best result, with no serious gross errors in the field. 
Unfortunately, there is no way to determine the sign of v and w without further information, their 
technique only revealing three of the six unknowns for the field. 

Since the present field has five unknowns, the direction of u being known if data is taken 
upstream of breakdown, these equations were rewritten for the probe configuration of the two 
wires parallel to the oncoming stream and at ±45°, according to the geometry of Figure 3.8b. 
Details can be found in Appendix D. The resulting equations are similar to that obtained by Sherit 
and Pletcher. They are also a set of three nonlinear equations with three unknowns plus a fourth 
equation. The direction of v and w could also be derived from this geometry as will be shown 
shortly. The intent, as noted previously, was to minimize the number of surveys required. 
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Unfortunately singularity problems occurred when the values of v and w approached each 
other in magnitude, as noted in Appendix D, and resulted in areas of the flowfield which were 
grossly in error. In addition, the values of u were not topographically representative of the field, 
although the maximum magnitudes were comparable. A Newton iteration scheme relived the 
problem somewhat, however artificial constraints were required to get rid of the singularities and 
the side effects of these were unknown. 

It was therefore decided to return to the method of Sherif and Pletcher. In order to fully 
determine the three velocity components and their associated directions, however, it was necessary' 
to take four grid sweeps above the wing at each chordwise location. Probe 1 , a DIS A 55P62, had 
wires lying in a plane perpendicular to the probe axis, while the wires of probe 2, a DIS A 55P61 , 
were lying in a plane parallel to the probe axis, as depicted in Figure 3.8a and 3.8b respectively. 
The initial two grid sweeps used probe 1 with wire 1 at the reference of zero degrees and wire 2 at 
negative 90 degrees using the geometry in Figure 3.8 and in accordance with that of Sherif and 
Pletcher (1987). The probe was then rotated 45° ± 0.5° about its axis and a second sweep 
initiated.This provided enough information for the velocity magnitudes to be determined. The 
second probe was used to take two sweeps with the plane of the wires parallel to the wing and 
perpendicular to it respectively. This second set of sweeps determined the direction of the 
transverse (v) and normal (w) velocity components. The direction of u was always assumed to be 
in the positive direction, as the probe was kept in the flow forward of the breakdown region. 

The complete derivation of the equations describing the velocity components in terms of the 
measured voltages can be found in Sherif and Pletcher (1987). The major equations will be noted 
below. In addition the expressions used to determine the directions of the v and w components 
based on the probe geometry used for the third and fourth spatial passes will be detailed. The 



geometry in Figure 3. 12 will be referred to and corresponds to that of Sherif and Pletcher, 
however the notation is in accordance with the convention present here. 
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Figure 3.12 Slanted Hot-wire Geometry 

The terms present in Jprgensen's equation can be represented as follows 

U N = u sin0 - v cos0 cos0 + w sin0 cos0 
U T = - u cos0 - v cos0 sin0 - w sin0 sin0 
U B = - v sin0 + w cos0 


(3.9a) 

(3.9b) 

(3.9c) 
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At 0 = 90° and the following reduced forms of equations (3.9 a - c) are derived, 

U N = u (3.10a) 

U T = - vcosG - w sin0 (3.10b) 

U B = - v sin0 + wcos0 (3.10c) 

If these are substituted into Jqrgensen's equation, a general expression for the measured effective 
velocity can be obtained. 

U e 2 = v 2 (ki 2 cos 2 0 + k2 2 sin 2 0) + w 2 (ki 2 sin 2 0 + k 2 2 cos 2 0) + u 2 

+ v w (sin20(k 2 2 - k^ 2 )) (3.11) 

Substituting in for values of 0 corresponding to rotation angles of 0°, 90°, and 45° gives 

U e 2 (0 = O°) = v 2 k^ 2 + w 2 1^2 + u 2 (3.12a) 

U e 2 (0 = 90°)= v 2 k 2 2 + w 2 kj 2 + u 2 (3.12b) 

U e 2 (0 = 45°)= \ v 2 (k! 2 + k 2 2 ) + \ w 2 (k! 2 + k 2 2 ) + u 2 + v w (k 2 2 - k\ 2 ) (3.12c) 


These equations can now be solved for u, v and w, 

qu e 2 (90°) - u e 2 (0°)}+ [(U e 2 (90°) - u e 2 (0°)} 2 + (U e 2 (90°) + u e 2 (0°) - 2U e 2 (45°)) 2 l °- 5 ) 

2(k 2 2 -k! 2 ) 

(U 2 (90°) + U 2 (0°) - 2U 2 (45 0 )) 

w = — £ — * 2 

2 (k 2 2 - ki 2 ) v 


(3.13a) 
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In order to determine the directional sign on the v and w velocity components, two 
additional spatial sweeps were required using the Figure 3.8b probe configuration. This probe 
was rotated about its axis 90° to obtain position 2 from position 1 as depicted in Figure 3.13. 


Figure 3. 13 Directional Determination Geometry 


From the above geometry, the value of y, corresponding to the angle between the projection 
of the velocity vector on the u-w plane and the x axis, can be seen to vary ± 90°, positive being 
defined as the particular angle shown. It can be seen that if the projected vector lies anywhere in 
this region, the effective velocity sensed by wire 4 will be greater than that of wire 3. Actually, a 
unique value can only be determined for 0° < y < 45°. If y is equal to say 55°, this would give the 
same readings on wires 3 and 4 as at y = 35°. Nonetheless, a greater effective velocity on wire 4 
than wire 3 would indicate than there is a positive w velocity component and that is what is 




required. Similarly, a greater effective velocity on wire 3 than wire 4 indicates a negative w 
velocity. Mathematically, this can be represented as. 
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y = 


45° - tan-1 


f u eff 3 ^ 

u eff 4 

v 


The sign of ythus determines the sign of the w velocity component. In a similar manner, a relation 


for then sign of v can be determined 


a = 45° - tan'l 


f Ueff 0 

u eff? 

L ) 


A positive value of a indicates a positive value of v, using the coordinate system shown. As an 
example of the effectiveness of this directional determination, an a map and a y map of the field 
above a 75° delta wing at 20° angle of attack are given in Figure 3.14. 



Figure 3.14 Directional Maps a) a b) y 
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The complete procedure for acquiring a field of data therefore combines Sherif and 
Pletchers technique with two additional grid sweeps for direction. In order to check the validity of 
this method, outside of comparing data acquired over a delta wing with published data, the probe 
was tested in a known set of flowstreams. The probe was oriented to arrive at six different sets of 
velocity components listed below: 


1. U = Uoo, v = w = 0 


2. u = v 


Uoo 

V2’ 


w = 0 


3. u = w = 0, v = Uoo 


, n Uoo 

4. u = 0, v = w = — 

V2 

. U ~ it 

5. U = W = — =r, V = Uo 

V2 


6. u = v = w = 


Uoo 

V3 


The difference between what the velocity values, scaled by Uoo, should be and what the procedure 
outlined above resulted in, is shown in Figure 3.15. 



Position 


Figure 3.15 Cross Wire Positional Response 
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Obvious differences existed between the measured and actual values. The largest discrepancies for 
the w component at position 3. Sherif and Pletcher note that this method is suitable for three 
dimensional flows and position 3 contains two component which should be zero. Hence, one 
would expect the largest error to occur there. Conversely, the condition where all the components 
are of equal value would be expected to yield the lowest error. This is not the case, for position 6 
indicates errors approaching 0.2. The lowest errors occur in positions 2 and 5. Position 5 
contains values of roughly equal magnitude, but position 2 has w = 0. Note that the values of u 
and v in position 2 are equal. For positions where u = 0, that being positions 3 and 4, a negative 
square root arose from the equations given previously and therefore the points are not displayed. 

If such a condition is used to denote a value of u < 0, that is, simply a conditional check, then the 
error difference would be zero. Although only one test was performed at each orientation, it is 
evident that as the flow departs from three dimensionality, where all three components are of the 
same order, the error increases dramatically. The equations (3.13 a - c) were further examined for 
their output sensitivity to changes in the measured input voltages. Since this was done using data 
obtained from the measurement surveys, discussion will be deferred until Chapter 5. A listing of 
the sensitivity inputs and outputs is provided in Appendix C. 

It should be noted that the important matter of temperature calibration of the wires has been 
circumvented. This is a very important consideration, as a calibration of a wire is not independent 
of the ambient temperature. To avoid the necessity of temperature compensation due to ambient 
changes, the wires were calibrated individually prior to and during each run. The ambient was 
constantly monitored and did not vary by more than ± 1°C which has a negligible effect on the 
measured voltages during the course of any one run. For information on temperature 
compensation techniques the reader is referred to Bearman(1971), Drubka, Tan-atichat, and Nagib 
(1977), Machen (1986), and Manca, Mastrullo, and Mazzei (1988). 



CHAPTER IV 


EXPERIMENTAL APPARATUS AND PROCEDURES 


4.1 The Wind Tunnel 

The experiments in this study were performed in the University of Notre Dame subsonic 
wind tunnel facilities. The tunnel is of the indraft, open circuit type and consists of the 24: 1 
contraction inlet, a test section, and the diffuser section as depicted in Figure 4.1. 



The test sections are interchangeable and are typically 610 mm by 610 mm (24 in. by 24 in.) 
cross-section with a length of 1820 mm (72 in.). The diffuser section expands the flow 
downstream of the test section through a 4.2 meter (13.8 foot) length at an included angle of 4.2°. 
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The tunnel is powered by an 18.6 kw AC induction motor which drives an 8 bladed 1 .2 meter (4 
foot) fan located in the diffuser outlet . Twelve "anti-turbulence" screens are mounted across the 
entrance to the contraction inlet to reduce the flow irregularities to a scale where they rapidly 
dissipate. 

The primary benefit of such a design is that it allows for flow visualization using a wide 
variety of tracer materials without contamination of the flow as would occur in a closed circuit 
design. The major disadvantage is susceptibility to atmospheric disturbances. Any variation of the 
pressure outside causes the tunnel velocity to vary with time. To reduce the amount of 
unsteadiness in the flow due to outside gusting, a flow restricter constructed of 5 mm (0.2 in.) 
diameter plastic tubes, 200 mm (7.9 in.) long and mounted in a frame, can be inserted between the 
test section and the diffuser. Due to the required flow speeds for calibration during this 
investigation, the flow restricters were not employed. The turbulence intensity in this particular 
wind tunnel configuration has been determined by Brendel and Huber (1984) using a single wire 
hot wire anemometer to be below 0.5% for all speeds and configurations and less than 0. 1 % at all 
clean section flow speeds for disturbances with frequencies greater than 10 Hz. 


4.2 The Test Section 

The test section utilized for the majority of tests was previously designed and used by 
Payne (1987). The section dimensions were 610 mm by 610 mm by 1820 mm (24 by 24 by 76 
inches). The four degree of freedom probe traversing mechanism was integrated with the roof of 
the test section as illustrated in Figure 4.2. 
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The system was computer controlled in the spanwise, stream wise and vertical directions. The 
probe position could also be rotated about the pitch axis manually. The streamwise and spanwise 
translator motions were driven with Slo-Syn MO623-FD08 stepper motors powered by a Velmex 
amplifier/controller model 8202M1. Directional control in the vertical utilized a DISA 52B01 
Sweep Drive Unit in conjunction with a 52C01 stepper motor. The amplifiers themselves were 
controlled by a Macintosh II computer equipped with a National Instruments MIO-16 12 bit data 
acquisition board. Code was developed utilizing LABView software to maneuver the probe to any 
relative position given the desired distances. The particular setup required the use of all 8 digital 
I/O 0-5 volt lines for control. Further technical details can be found in Visser (1989). The 
minimum step sizes possible in the streamwise, spanwise, and normal directions as noted by 
Payne (1987) were 0.0064 mm, 0.0254 mm, and 0.0208 mm respectively with an overall 
positional accuracy of ± 1 mm. 
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4.3 Models 

All the models used for testing were flat plate delta wings having a windward 25° bevelled 
edge. The flowfield velocity measurement models were made of aluminum, having a centerline 
chord of 406.4 mm (16 inches) and a thickness of 6.35 mm (0.25 inches). The sweep angles were 
70° and 75°. A schematic of the coordinate system and the associated geometric details is illustrated 
in Figure 4.3. 



Figure 4.3 Delta Wing Model Geometries a) Aluminum Full Span b) Acrylic Full Span 
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The surface pressure measurements were taken using 0.25 inch thick, 70° sweep acrylic 
delta wing planforms. Three chord lengths of 304.8 mm, 406.4 mm, and 508.0 mm (12, 16, and 
20 inch) were investigated. Channels were milled into the surface, allowing 1.37 mm ID / 1.83 
mm OD (0.054 / 0.072 inch) stainless steel tubing to be embedded along rays of y/s = 0.5, 0.6, 
and 0.7 on each of the models. Pressure tap holes of an outer diameter of 1.07 mm (0.042 inch) 
were then drilled at chord stations ranging from 0.2 to 0.9 as shown in Figure 4.3b. When 
measuring the surface pressure, all the taps were covered by a piece of transparent tape, save the 
station being measured. 

Additional velocity measurements were made using a half span aluminum model of the 
same sweep and bevel as the full span aluminum models, but a root chord of 26.95 inches. This 
enabled the effect of a finer grid resolution to be examined in conjunction with a smaller relative 
probe size. A splitter plate was employed during these tests. 


4.4 Flow Visualization Techniques 

Three types of flow visualization were applied during the course of this investigation. 
Kerosene and titanium tetrachloride were used to mark the flow externally. Visualization of the 
planform surface was accomplished using a mixture of kerosene, oleic acid and titanium dioxide. 

Kerosene was vaporized on electric resistance heater strips to produce a white smoke 
which, in conjunction with the correct lighting, visualized the flow. Four of these systems were 
combined to form a smoke generator. A squirrel cage blower forced the smoke through a series of 
cooling pipes to reduce possible buoyancy effects. The smoke was exhausted through the rake of 
tubes and entered the tunnel upstream of the contraction section. An extensive review on the 
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subject, both past and present, is given by Mueller (1978) and pertinent schematics can be found in 
Visser (1988). 

Titanium tetrachloride (TiCl 4 ) is a colorless liquid that, upon exposure to water vapor, 
reacts to form a dense white smoke composed of hydrochloric gas and titanium dioxide particles. 
This smoke can be introduced into the flow field to locally visualize certain aspects of the fluid 
motion. Other applications have included dripping or painting the liquid onto a model enabling 
brief periods of visualization before the supply had to be replenished (Freymuth, Bank, and Palmer 
(1985)). In the present tests, titanium tetrachloride was contained in a specially constructed glass 
flask under an inert nitrogen atmosphere. The low vapor pressure of the TiCl 4 allows the region 

above the liquid to saturate with the vapor component and thus the gaseous space inside the 
container contains a mixture of the nitrogen and the TiCl 4 vapor. Nitrogen was fed into the 

container, displacing the vapor in the bottle out to the test section where it exited at the model. The 
tunnel probe acts in a similar manner to a water tunnel dye marker in that the vapor was 
immediately visible, marking a particular streakline. Further details are described by Visser (1988) 
and Visser, Nelson, and Ng (1988). 


Surface visualization was conducted on the 70° and 75° degree models. Each was fitted 
with a mounting pin at the trailing edge to facilitate removal and subsequent record on film. A 
mixture of 15 parts of kerosene, 5 parts of titanium dioxide and 1 part oleic acid provided an oil 
based slurry which was spread uniformly on the model at zero tunnel velocity. The tunnel was 
then run at the test condition until the liquid evaporated leaving the skin friction lines visible. The 
models were subsequently removed from the tunnel, photographed and cleaned prior to the next 
test Tests were conducted at a Reynolds number of 250,000 while the angle of attack was varied 
from 25° to 45°. Both the leeward and windward surfaces were documented. Tests with the 
models reversed, such that the bevel was on the leeward side, were also performed. 



Still and moving visual recordings were made of the flow over the models to help in 
analyzing the data. Photographs were taken using a Nikon FM2 35mm SLR camera. These were 
used to record the various surface visualization patterns at different angles of attack. Kodak Tri - 
X Pan 400 ASA black and white print film was used for the photographs. Video taping was 
accomplished using a Panasonic Digital 5000 System Camera capable of an effective frame rate of 
1/1000 of a second. A Panasonic NV - 8950 VHS recorder was used to document the events. 
Lighting was in the form of high intensity 1000 Watt lamps placed so as to maximize tracer 
visibility. 


4.5 Pressure Measurements 

Surface pressure measurements were made on the 0.25 inch thick, 70° swept acrylic delta 
wings of 12, 16, and 20 inch chord length. Transparent tape was used to cover the all the pressure 
tap holes save the one which was being measured. The Reynolds number was held to 250,000 
and the sting location was also varied. Angles of attack ranged from 25 to 45 degrees. 

Total pressure measurements were made in the core of the vortex above the 75° swept delta 
wing at angles of 20° and 30° degrees angle of attack. Various chord stations were examined, all at 
a Reynolds number of 250,000. A stainless steel tube probe of 1.07 mm ID / 1 .47 mm OD (0.042 
/ 0.058 inch) used to measure the total pressure and is shown in Figure 4.4. To obtain the total 
pressure, the probe was maneuvered to the core center, based on the value of the maximum axial 
velocity obtained from the hotwire measurements. The probe location was then manually adjusted 
in the y/s and z/s directions until a minimum value of of the total pressure was observed. 
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4.6 Cross-wire Equipment 

The majority of velocity field data was accumulated over a 75° sweep delta wing at 20 ° 
angle of attack. Measurements were made at various chordwise stations in grid planes normal to 
the upper surface. The angle of attack was then increased to 30° and chordwise stations were 
measured upstream of the probe induced breakdown of approximately x/c = 0.5. Experiments by 
Payne, Ng, and Nelson (1987) comparing LDV and seven hole probe data have shown that the 
effect of introducing a probe into the flowfield does not greatly distort the flowfield provided the 
measurements are taken upstream of the breakdown zone. Flow visualization confirmed that 
positions measured were kept upstream of the breakdown region. A 70° sweep configuration was 
utilized at 20° angle of attack in an effort to observe the effect of sweep. A final set of tests were 
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taken with the larger half span model to observe the relative effects of grid resolution changes . 
Unless specifically noted, all the data presented reflects a Reynolds number of 250,000. 

The geometry of the probe used for acquiring the hot wire data is shown in Figure 4.5. A 
DISA geared probe holder held the probe extensions out to the actual cross wire probe. As detailed 
earlier, two x-wire probe configurations were utilized and were illustrated in Figure 3.8. Probe 1 . 
a DISA 55P62, had wires lying in a plane perpendicular to the probe axis, while the wires of probe 
2, a DISA 55P61, were lying in a plane parallel to the probe axis. The probes utilized five 
micrometer diameter tungsten wires giving an length/diameter ratio of 250. Overheat ratios were 
set to 1.8. The wires were calibrated for every test to reduce the possible errors associated with 
property changes of the wires. In addition, this procedure eliminated the need for temperature 
compensation as the ambient temperature did not vary by more than ± 1° C over the course of any 
individual test.The minimum distance between the wing surface and the probe was 3.0 mm due to 
probe geometry. 



Figure 4.5 Hot Wire Probe Geometry 
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In order to fully determine the three velocity components and their associated directions it 
was necessary to take four normal grid sweeps above the wing at each chordwise location. The 
initial two grid sweeps used probe 1 with wire 1 at the reference of zero degrees and wire 2 at 
negative 90 degrees using the geometry in Figure 3.12 and in accordance with that of Sheri f and 
Pletcher (1987). The probe was then rotated 45° about its axis and a second sweep initiated.This 
provided enough information for the velocity magnitudes to be determined. The second probe was 
used to take two sweeps with the plane of the wires parallel to the wing and perpendicular to it 
respectively. This second set of sweeps determined the direction of the transverse (v) and normal 
(w) velocity components. It was assumed that the direction of u was always in the positive 
direction as the probe was kept in the flow forward of the breakdown region. 

The cross wire probes were monitored using a TSI IFA 100 Model 150 constant 
temperature anemometer system. The accuracy of the anemometer given by the manufacturer is 
such that as the probe resistance is brought to zero, the difference on the actual resistance measured 
is: 

meter Q 

Probe Resistance 0000 8.500 

0010 8.504 

Hence, if the worst case scenario is presumed with the meter capable of displaying ± 0005 or 
0.002 Q then an estimate of the desired overheat ratio can be obtained. Maximum error for an 
overheat ratio of 1.8 is estimated to be ± 0.08%. Details can be found in Appendix C. A TSI 
model 570 signal conditioner was used in conjunction with the anemometer to provide a DC 
coupled offset of 1 volt ± 0.15%. The 570 signal conditioner applied a gain of 5 ± 0.15% to the 
signals and low pass filtered at 1000 Hz ± 10% before they were sampled by the computer. 
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4.7 Data Acquisition and Reduction System 

Data was acquired using a Macintosh II computer equipped with a National Instruments 
MIO-16H-9 12 bit D/A board. All the signals were obtained in the differential mode using uniform 
gains and a -5 to 5 volt range. The available precision for these settings corresponding to one least 
significant bit of the 12 bit A/D converter was 2.44 millivolts. The sampling rate was set to 25,000 
Hz. Two channels were used to acquire the hot wire voltages and the dynamic pressure was 
obtained on the third channel with a resolution of 1.22 millivolt/bit. 

The data was reduced to velocity values using a least squares fit based on the calibration to 
determine the effective velocities. The method of Sherif and Pletcher (1987) was then applied to 
determine the magnitudes of the velocity components. Direction was determined at the same time 
using the information form the third and fourth wire configurations. 

Several reduction codes were written for the post processing of the data. These codes 
accomplished a wide variety of tasks and, for simplicity, a schematic overview is given in Figure 
4.6. 


3-D Traverser Section 
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Reduction Code 



Figure 4.6 Data Acquisition and Reduction Schematic 





CHAPTER V 


RESULTS AND DISCUSSION 


The results of this study are divided into several sections, beginning with the surface flow 
visualization. Several interesting features were noted from the visualization, which have not been 
reported by previous investigators and so are presented here. The major part of this chapter deals 
with the data obtained from the cross wire velocity measurements. First, the velocities are 
discussed and compared with previously obtained data. An indication of the unsteady nature of the 
flow is also revealed from the fluctuations in the measured wire voltages. The data is then 
analyzed in light of its derived properties such as the vorticity and the circulation. Dimensional 
scaling and various ways of interpreting the data are offered as a means for characterizing the 
behavior of the vortex. Parameter correlations are presented next, using both conventional ratios, 
such as the swirl angle, and ratios typically not examined in this context. A brief mention is made 
of some interesting spectral behavior in the vortex core upstream of breakdown and this is 
followed by a look at the core dimensions calculated from the cross wire data itself. Lastly, the 
surface and off-surface pressure data is detailed in light of the ability of these measurements to 
indicate the onset of breakdown. The total pressure measurements are combined with the velocity 
data to examine the concepts introduced at the end of the previous chapter. 


5.1 Surface Flow Visualization 

The interpretation of the 'skin friction' lines, that is the lines resulting from surface flow 
visualization technique, can indicate a great deal about the nature of the flow. Extensive 
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discussions on flow topology deciphered from delta wing surface flowlines are given by Peake 
and Tobak (1980) and Delery (1990). Investigators, such as Kjelgaard and Sellers (1990) have 
conducted tests to establish locations of the major features for comparison with numerical 
predictions. In the present study, surface flow visualization on the 70° and 75° planforms was 
performed to examine any changes on the surface flow patterns specifically due to breakdown. No 
evidence of transition to a turbulent boundary layer was observed. Since the Reynolds number 
based on chord was much less than 500,000, this was in agreement with the earlier statements of 
Naarding and Verhaagen (1988) that the entire boundary layer on the upper wing surface can be 
expected to remain laminar. Several other interesting features were noticed, however, and are 
presented below. 

Typically, investigators have found that the surface visualization will reveal the presence of 
secondary and even tertiary separation and reattachment lines on a flat plate delta wing. Carcaillet, 
Manie, Pagan, and Solignac (1986) applied topological rules to interpret their visualization 
photographs for a 75° sweep wing at a = 20. A cross sectional schematic of the interpretation of 
their visualization data is presented in Figure 5. la. The major features, including the presence of 
the primary, secondary, and tertiary vortex structures and their respective separation and 
attachment lines, were observed. Naarding and Verhaagen (1988) also observed the presence of a 
tertiary vortex system on their unit aspect ratio biconvex delta wing. They point out, though, that 
this tertiary flow structure is confined to laminar boundary layer flows. 

Carcaillet, Manie, Pagan, and Solignac further suggest that the attachment point of the secondary 
vortex is almost identical to the primary separation location. Thompson (1975) postulates that a 
vortex structure may exist in this region between these two nodes, as is detailed in Figure 5.1b. 
Based on this, Dixon (1989) concludes that the primary vortex is not fed by a continuous sheet of 
voracity, but is comprised of a stationary, shear layer vortex which feeds the main vortex. The 
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primary vortex 



leading edge 


(a) 



Figure 5. 1 Vortex Topology a) A=75° a = 20° (Carcaillet, Manie, Pagan, and Solignac, 1986) 


b) Proposed Stationary Shear Layer Vortex Thompson (1975) 



120 


present surface flow experiments indicate, however, that if some region of vortex structure exists 
between the secondary vortex attachment point and the leading edge, it would have to be rotating in 
a sense opposite to that in Figure 5.1b, based on the direction of the surface flow pattern. 

Present experimental results indicated a very complex surface pattern extending from just 
inboard of the secondary vortex separation line to the leading edge. The results of a surface flow 
visualization test at a = 30° for the A=75° delta wing planform is shown in Figure 5.2a. A close- 
up is shown in Figure 5.2b. The detailed features observed in the region between x/c = 0.3 to 0.4 
were handsketched in Figure 5.2c from the model itself. Although all the features can not be 
explained, three interesting observations can be noted. The first is the darker region on the surface 
just inboard of the secondary separation line. This is common to the flow of the majority of cases 
considered here. It is possible that this indicates a local acceleration of the flow, although one 
might expect the flow to be slowing as it approaches the point of separation. The second 
interesting area is that located just inboard of the leading edge. This too appears darker, but the 
surface lines are plainly seen. The secondary vortex attachment line appears to lie just inboard of 
this darker region and the flow is thought to be traveling outwards to rejoin the leading edge flow. 
However, the surface lines point in the downstream direction which, unless there is a reverse flow 
region, indicates the surface flow to be travelling towards the center of the wing. This would 
indicate that any rotating flow present in this area would have the same sense as the secondary 
vortex, opposite to that proposed by Thompson (1975). 

The final feature worth noting, which occurred on most of the planforms, involved an area 
of reversed flow on the surface of the wing. This occurred just outboard of the secondary 
attachment line. The area can be observed near the aft portion of the wing in Figure 5.2a. A closer 
look is given in Figure 5.3 for the 70° sweep wing at 25° angle of attack. The 'puddles left by the 




Figure 5.2a Leadinu Fdge Details tor A— 75° ex — ^0 , Plantomi 


122 


ORIGINAL PAGE lb 

OF POOR QUALITY 



Figure 5.2b Leading Edge Details for A=75° a = 30° , Closeup of Leading Edge 
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flow moving upstream are clearly evident. Earnshaw and Law lord (1964) also recorded much 
more pronounced evidence of now forward from the trailing edges for lower wing sweeps and 
greater angles of incidence. The flow patterns that are presented for their higher sweep wings 
appear to contain features similar to that in in Figure 5.3 although they do not discuss aspects of 

this specifically. 



Figure 5.3 Reversed Flow Region Details for A-70° oc - 25 

Gross changes in the surface pattern due to the movement of breakdown were not 
noticeable as shown on the 70° wing in Figure 5.4a, where no breakdown is present on the wing, 
and 5.4b, in which breakdown occurs at approximately x/c =0.4. 



Figure 5.4a Angle of Attack Effect for A=70°, a = 25°, no breakdown present 


The only real feature that appears to have changed is the region just inboard of the leading edge. At 
the higher angle of attack, this region begins to widen closer to the apex. Note also the darker 
region inboard of the secondary vortex, mentioned previously. At a = 25, the region extends to 
the trailing edge, but at a=35° it is less evident aft of about x/c = 0.7. It is not until a = 40°, 

Figure 5.4c, that the surface shows a major topological change and the effects of breakdown can 
be more dramatically seen. 
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The breakdown is now at 20% chord and all traces of the darker region on the inboard side of the 
secondary separation line have disappeared aft of about x/c = 0.4. A weak separation line is still 
evident, however the complex structure noted earlier between the leading edge and this line is all 
but disappeared. 



Figure 5.4b Angle of Attack Effect for A=70°, a = 35°, breakdown at x/c = 0.4 

Similar features were observed on the 75° wing. The dark region inboard of the secondary 
vortex was much more pronounced. Figure 5.5a is a closeup of the apex region, up to about 10% 
of the chord, at an angle of attack of 35°. The darker area is seen to extend from the apex. The 
bottom surface near the apex, as shown in Figure 5.5b, exhibited dark regions on the surface near 
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the apex as well. This darker area, on the underside of the planform, decreased in size with angle 
of attack and vice versa. There did not seem to be any local change in the direction of the surface 
lines between this region and the lighter colored areas next to it. A suggested explanation is that 
the solution for marking the flow simply dried at a faster rate in these areas, which would indicate a _ 

larger local velocity. 



Figure 5.4c Angle of Attack Effect for A=70°, a = 40°,breakdown H x/c = 0.2 
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Figure 5.5 Apex Details, x/c - 0 to 0. 1 A=75° a - 35° a) Upper Surface b) Lower Surface 
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The surface flow visualization experiments provided several interesting insights on the 
nature of the topology of the vortex flow field. As demonstrated by other investigators, the 
presence of a tertiary vortex system is apparent. Local darker surface regions on the underside of 
the wing near the apex were noted to decrease with angle of attack. Additional dark regions 
inboard of secondary vortex separation line, on the leeward side of the wing, were seen to persist 
at all angles of attack until complete separation occurred. At higher angles of attack, these dark 
regions did not appear to extend completely to the trailing edge. These regions extended and 
intersected each other at the apex. Finally, the surface flow directly inboard of the leading edge 
and outboard of secondary vortex attachment line indicates a flow direction inwards from the 
trailing edge, towards the wing centerline, and in the axial direction. If this was caused by a local 
vortex, it would have the same sense of rotation as the secondary vortex. 


5.2 Cross-wire Measurements 

Velocity field data was acquired at a series of chordwise stations using the cross wire 
technique described earlier. A sample of the surveys are presented and compared with data 
acquired by Payne (1987) using a seven hole probe. The changes in the velocity profiles with 
distance from the apex is then examined along with the unsteadiness of the measured field. The 
velocity data was subsequently manipulated to derive field properties such as vorticity and 
circulation. These are described in the sections following. Dimensional scaling is used extensively 
to examine the behavior of these properties as the chordwise location and angle of attack were 
varied. Both integrated or average field distributions as well as discrete values are detailed. 

Finally, correlation parameters, such as the Rossby number and extensions of the discussions in 
Chapter 2, are presented. A listing of the data matrix, along with as some of the measured 
properties, is given in Appendix B. 
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5.2.1 Velocity 

A typical axial velocity distribution measured by the present technique for A = 75° at an 
angle of attack of 20° and x/c = 0.5 is given in Figure 5.6a. Note that the distances in the y and z 
direction are scaled by the local semispan. Thus, a y/s = 1 .0 corresponds to the leading edge of 
the wing. The axial velocity u/U^ is the velocity normal to the measurement plane, which was 

perpendicular to the chord line, at the given chordwise station. The jetting core structure of the 
vortex is quite well defined, with the majority of the measured field maintaining a velocity above 
the freestream velocity. Measurements encompassed a z/s of 0.055 to approximately 1.0 and 
were taken from the chord centerline out to a spanwise location of y/s =1.2 beyond the right 
edge of the wing. Grid increments were set to a y/s and z/s of 0.03 at each station. 

Comparison with previous data obtained by Payne (1987) using a seven hole probe (SHP) 
was used to evaluate the performance of the cross wire measurement technique. The inherent 
complexity of the equations given in Chapter 3 make an assessment of the error difficult for each 
variable in the equation. In addition, the values of ki and k2 were taken to be constant at 0.3 and 

1.1 respectively, when they are in fact a function of the flow vector at extreme angles. Finally, 
four spatial passes were required for each survey, and although each was carefully aligned and 
each traversal was from the centerline outboard, the error in the traversing system contributes a 
cumulative error for each station. Thus a quantitative comparison with other data was used 
appraise the technique. 

Seven hole probe data for the same test conditions in Figure 5.6a is given in Figure 5.6b. The 
field of Figure 5.6a is also rendered in color in Figure 5.6c for comparison. Both sets of data 

compare favorably to each other, however several differences are observed. The cross wire data 
exhibits a peak axial velocity of u/U^ = 2.33, which is about 6% greater than that indicated in 

the SHP data. Both profiles exhibit larger velocity gradients on the lower side of the vortex. 
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The core locations, based on the position of maximum velocity, for the cross wire data, y/s = 
0.69 and zjs = 0.40, lie further outboard and away from the wing than the SHP data, y/s = 0.62 
and z/s = 0.34. This is believed to be probe induced, the cross wire probe being approximately 
three times larger than the seven hole probe. The core location for the seven hole probe data, 
based on the position of the lowest value of the total pressure coefficient, Cp t> was located at 
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Figure 5.6 Axial Velocity A=75° at a = 20° c) Cross Wire Color Image 

y/s = 0.66. The resolution of the SHP data is 0.04167 in both directions, which, if applied to 
the core location at y/s = 0.66, would account tor both the y/s = 0.62 of the SHP velocity 
contour and y/s = 0.69 of the cross wire data. 
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Directly underneath each vortex lies a regime of flow moving at below freestream velocity. 
This region is larger on the SHP profile and appears to be distorting the vortex above it. This is 
in contrast to the cross wire data. It is unknown whether this is a phenomenon of the flow or a 
direct result of the interference of the probe. Examination of the SHP total pressure profiles 
taken by Payne (1987) revealed values of Cp t > 1 .0 in this region, which is physically 

impossible. This would indicate that the interference effects associated with the nearness of the 
probe to the wall are substantial, possibly for both probes. Surface visualization has indicated 
that the flow is approaching the secondary vortex separation point in this region and there is a 
sizable component of velocity in the spanwise direction (Figure 5.4a, 5.5a). This turning of the 
local velocity vector can account for the decrease in the axial component. The same effect can be 
seen just outboard of the leading edge, y/s = 1 .0. The presence of the shear layer itself can be 
detected from the data in Figure 5.6b, at y/s = 1.0, but is not nearly as evident in the cross wire 
data in Figure 5.6a. 

At a = 30, the vortex structures are seen to be stronger as presented in Figure 5.7. Lower 
than freestream flow is again present below each vortex. A stronger evidence of the shear layer 

is seen in Figure 5.7a. The cross wire data again indicates a greater maximum axial velocity, 
u/U^ = 3.2, versus 2.98 for the SHP data. The core locations for the cross wire data now lie 

inboard of the SHP locations. The total pressure data from the SHP experiments indicate a z Js 
location of 0.43 for the maximum Cpt compared to 0.47 based on a maximum axial velocity. 
With this in mind, the core locations can be seen to coincide, allowing for the possible error of 
half the resolution of each respective grid increment. 

A third comparison of the cross wire data with the SHP data for a different wing planform 
is made in Figure 5.8. In this particular instance, A = 70° and a = 20°, the data compares very 
well. The core locations are identical, both laterally and normal to the wing. The peak velocity 
of the cross wire data exceeds the SHP data by less than 2%. In addition, the size of the low 
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axial velocity regimes underneath the vortices are seen to be very similar. The most outstanding 
difference is that the u/U^ = 1 .2 contour does not encircle the vortex in Figure 5.8b as in 5.8a. 

The in-plane velocities, v/U^ and w/U^, for the A = 75°, a = 20° and x/c = 0.5 survey are 

contoured in Figure 5.9a and 5.9b respectively. Negative values, following the convention in 
Chapter 1, are represented by the dotted lines. For both components, the greatest error can be 
seen to occur in the region where the velocity approaches zero. More correctly, it occurs where 
the velocity component experiences a change of sign. The contours do not appear to be evenly 
spaced in this region but seem to crowd the zero contour. The sensitivity of the crossflow 
components to changes in the measured flows can be severe in regions where v and w are much 
less than u. Along the core, an increase in the input voltage of wire 1 by the standard deviation 
of the measured voltage itself resulted in less than a 1% change in the calculated value of u, but 
over an 80% decrease in the magnitude of v. This was equivalent to about 23% of the freestream 
value. Regions of the flowfield containing velocity components of similar magnitude indicated 
changes in the velocites of approximately 15% for changes to the input voltages of one standard 
deviation of the measured signal. Details can be found in Appendix C. 

A final, more quantitative comparison, was made between the cross wire and the SHP data 
by comparing the cross sectional traverses through the vortex core. The axial and transverse 
velocities from this traverse, corresponding to the A = 75°, a = 20°, x/c = 0.50 data in Figure 

5.6, are plotted in Figure 5. 10a against spanwise distance. The z/s location for the traverse was 
based on the maximum value of u/U,^. Both of the cross wire velocity components can be seen 

to be shifted outboard of the SHP data. The axial peak differences are equivalent to that noted 
earlier. Inboard of the core center, the tangential velocity from the cross wire data maintains a 
larger, negative value. In contrast to this, the cross wire data indicates a smaller magnitude value 
of the axial velocity component, in the area between the core center and the chord centerline. 

This trend is also observed in the the cross sectional traverses of Figures 5.7 and 5.8 as plotted 
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Figure 5.10 Axial and Tangential Velocity Core Profiles a) sweep = 75, alpha = 20 
b) sweep = 75, alpha = 30 c) sweep = 70, alpha = 20 
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in Figures 5. 10b and 5. 10c respectively. There is no major indication of a shift in the location of 
the core for either 5.10b or c. There is a noticeable increase in the maximum magnitudes of the 
tangential component for the A = 70° planform in Figure 5.10c. One might attribute this to grid 
resolution differences, however data was acquired at similar spatial locations for each 
experiment. More likely, the differences lie in the measurement technique. Since the flowfield is 
sensitive to probe intrusion, and the cross wire technique used a larger probe, the error would 
seem to arise in this measurement. 

Based on the above comparison, the cross wire data was concluded to compare well with 
the SHP data. Payne ( 1987) also showed the SHP data to correlate well with LDV data in the 
same study for the flowfield upstream of breakdown, and thus the present velocity measurements 
represent the flowfield well. The cross wire data was then evaluated with respect to the 
repeatability of the data itself. 

Three of the chordwise stations, where data was acquired above the 75° planform, are 
shown in Figure 5.1 1. Discrepancies were observed. The largest variations were observed 
nearest the apex, as in Figure 5.11a at x/c = 0.3. Both the axial and tangential profiles indicate 
deviations in magnitude along certain portions of the profiles. Aft of this station, the axial 
profiles appear very repeatable. Peak values in the tangential velocities did show variations, as 
noted in Figures 5.1 lb and 5.1 lc, of up to 10% in the worst cases. 

For the data of Figure 5.1 lb, which represents three independent velocity surveys, profiles 
were examined at a series of z/s locations corresponding to locations near the wing surface and 
between the surface and the core location. The u, v, and w component are plotted in Figures 
5.12, 5.13, and 5.14 respectively. Repeatability, including sudden peaks and drops is good. 

The axial velocity profiles show the smallest deviations. A sharp, very repeatable drop from 
positive to negative velocity is recorded in the v component as indicated in Figure 5.13a. This 
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Figure 5.12 Axial Velocity Profile Comparison a) Surface b) Midcore c) Core 
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Figure 5.14 Normal Velocity Profile Comparison a) Surface b) Midcore c) Core 
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was observed at almost every chordwise location. Comparison of this spanwise location of 
approximately y/s = 0.85 with the surface visualization results in Figures 5.1 and 5.2b indicates 
this location to be directly beneath the secondary vortex. The reversed flow direcdon is therefore 
to be expected. 

The largest deviations are seen in the w velocity component at a position near the wing 
surface. Figure 5.14a, and the v component through the core. Figure 5. 13c. Two of the profiles 
in Figure 5.14a give indication of incorrect signs in the region of y/s = 0.4 to 0.6. It should be 
kept in mind that the technique used to reduce the data is most suited to strongly three 
dimensional flows and suffers if one component is much less than the others. 

The chordwise variation of the velocity data profiles is given in Figure 5. 15. The axial and 
tangential components obtained from the core traversals are overlaid for chordwise stations of x/c 
= 0.3 to 0.8. Each profile is plotted against its local spanwise direction. All the profiles exhibit 
similar characteristics, both in shape and magnitude. This indicates the flowfield in this region to 
be scaling in a linear fashion, indicative of a conical behavior. The largest differences occur 
outboard of the core location. For consistency, every spanwise pass for each survey was made 
in the direction from the chord centerline to the leading edge. As the probe encounters the core 
region, it has the tendency to 'push' the core slightly ahead of itself as it is encountered. If the 
probe is exactly on the z/s location of the core axis, this results in a slight displacement of the 
core laterally which may offset the position of the maximum velocity recorded. However, if the 
probe is above or below this height, the spatial location of the core is displaced below or above 
its natural location. In addition, each velocity profile in Figure 5.15 represents the z/s location 
where the maximum u/U°° value was recorded. Because ot resolution, this may not be the 
maximum in the field. The same applies to the tangential velocities. 
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Two other aspects of interest that result directly from the velocity field data can be 
immediately examined. The behavior of the maximum axial velocity in the core, with respect to 
magnitude as it proceeds in a chordwise direction, has been a source of discussion for several 
investigators, as noted in Chapter 1. The majority of the data was acquired for the 75° sweep 
planform at 20° angle of attack and the maximum axial velocity components are compiled in 
Figure 5.16. No breakdown was occurring on the wing in this configuration and the data falls 
roughly in a band of u/U°° = 2.0 to 2.5. The values crest in magnitude at about x/c = 0.6 and 
become less to either end of the wing. The grid resolution for this data was y/s = z/s = 0.03 at 
each station. Therefore the grid increment size in actual units is increasing towards the apex. If 
the subcore region behaves in a cylindrical manner, as has been proposed by other investigators, 
the resolution of the subcore flow should also increase as the apex is approached. The direct 
implication of this is that the decrease in the axial velocity towards the apex represents a real 



*, ** See below 

Figure 5.16 Chordwise Distribution of Maximum Axial Velocity at A = 75° a = 20 



147 


phenomenon. Since breakdown was not present on the planform, the decrease in axial velocity 
towards the trailing edge could be a direct result in the decrease in the actual grid size. It should 
be remembered, however, that the actual resolution ot the probe is increasing towards the trailing 
edge, since the probe size is fixed and the vortex size is increasing. The data marked with an 
asterisk and the double asterisk represent increases on the resolution of the survey grid by two 
and three times, respectively. Similar trends are seen for both configurations and although the 
y/s = 0.015 case, *, at x/c = 0.7 is greater in magnitude, it too falls onto the y/s = 0.03 data by 
x/c = 0.8. 


The A = 75° at a = 30° data, in Figure 5. 1 7, shows a large drop in the axial velocity aft of 
the x/c = 0.4 mark. Breakdown was occurring at approximately x/c = 0.55 to 0.60 and no data 
was taken rearward of this position. The velocity does drop over a distance of about 10% of the 
chord before breakdown. Data for the 70° planform is also included. Breakdown for the 70 
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Figure 5.17 Chordwise Distribution of Maximum Axial Velocity with Breakdown Present 
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wing data was observed at about the x/c = 0.5 to 0.55 station. Only a very slight indication of an 
imminent drop in the velocity was observed, however, at the x/c = 0.5 position. No data was 
available after the onset of breakdown. Both sets of data in Figure 5.17 suggest that the axial 
velocity does not vary appreciably in a region some distance forward of the breakdown zone. 

For comparison, SFIP data from Payne ( 1 987) for an 85° sweep wing at 40° angle of attack is 
also presented in Figure 5.17. There is only a small indication of the upcoming drop in the 
velocity due to the breakdown. Stations more than 5% forward of the breakdown location 
exhibit no major change of the maximum axial velocity. Upstream of the wing, the axial velocity 
component must begin with a freestream speed and then accelerate over the apex and along the 
core as the vortex develops. From the present data, the axial velocity continues to increase up to 
some distance before breakdown. A position is eventually reached where the velocity does not 
increase further, but reaches some maximum, remaining at that magnitude for some distance, 
until falling abruptly through the breakdown region. 


The last topic of interest dealing directly with the velocity data is the intensity of the 
fluctuations when the velocity signals were acquired from the hot wires. Since the derivation of 
the velocities required data that was not taken simultaneously, a direct determination of the 
turbulence intensity of the velocity can not be obtained. The mean voltage signal and the 
associated standard deviation for each wire, however, can be used as an indicator of the local 
fluctuation intensity of the flow. By taking the ratio of the standard deviation to the mean, and 
multiplying by 100, a percent fluctuation intensity was defined. Intensities for the 75° sweep 
planform at a = 20° based on the data from wires 1 and 2 of the first spatial pass are represented 

in Figures 5.18a and 5.18b respectively. With the probe in position two, two more sets of 
information were obtained. These correspond to Figures 5.18c and 5.18d. 

Since the orientation of the wire to oncoming flow is different for each configuration, as 
noted in each figure, the sensitivity to the local fluctuation is different and hence each plot shows 
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variations in the levels of intensity. The largest fluctuations occur in the region of the shear layer 
as it leaves the wing. The fluctuation intensity can actually be seen to mark the wrapping up of 
this layer above the wing as well, especially in Figure 5.18a and 5.18c. A maximum intensity of 
almost 14% was recorded with the wire parallel to the wing in Figure 5.18b. The orientation of 
the wire can actually be determined by closer observation of the figures themselves. In Figure 
5.18a, there is a drop in the intensity in the area of y/s = 0.5 and z/s = 0.4, the region where the 



y/s 


Figure 5. 19 Fluctuation Intensity at x/c = 0.5 for A = 75° at a - 30° 
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flow is oriented along the axis of the wire. This would be true in the flow coming over the wing 
as well, however the effect is masked by a much stronger leading edge flow. The wire parallel to 
the planform in Figure 5.18b shows its insensitivity in the regions of flow parallel to its axis as 
well, namely at y/s = 0.7 and z/s = 0.2 and 0.6. Similarly, Figures 5.18c and 5.18d reveal 
fluctuation deficits along lines at 45° to the plane of the wing. Even stronger fluctuations were 
recorded above the planform at an angle of attack of 30°. A sample of this data is given in Figure 
5.19. Intensities as large as 16% were noted. 

The contours suggest that another region of strong fluctuations exists in the secondary 
vortex. To facilitate observation of this. Figure 5.18b was enlarged and is presented with a color 
contour of the flow in Figure 5.20. Directly inboard of the leading edge lie pockets of intensity 
on the order of that in the shear layer. Whether this is due to the local influence of the shear layer 
or not is unclear. As indicated in Figure 5.6a and Figure 5.9 the velocity vector in this region is 
considerably lower in the secondary vortex than that in the shear layer. Hence, any local 
fluctuations register as larger values when normalized by their average, whereas on an absolute 
scale they would in fact be smaller. 

The intensity data reveals some additional noteworthy characteristics of the vortex flow. 
Strong fluctuations are present in the secondary vortex flow regime as well as the shear layer. 
This may be an indication of the naturally turbulent nature of the flow in these regions. The 
boundary layer induced on the wing by the primary vortex is laminar and any transition to 
turbulence would have to occur outboard of the secondary separation point. It has also been 
suggested that some type of periodic phenomena may exist in the vicinity of the leading edge 
flow, such as discrete vortex shedding by Gad-el-Hak and Blackwelder (1985). Since a periodic 
phenomena would increase the magnitude of the measured fluctuations, it could be surmised that 
there may also be some type of periodic phenomena occurring in the secondary vortex. Finally, 
it is interesting to observe that as the shear layer is followed as it wraps up into the core region, 
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the level of the intensity of the fluctuations drops by over an order of magnitude to a virtually 
negligible value along the axis of the vortex. 

Based on the overall consistency of the velocity data above, and the good correlation with 
the SHP data, is was determined that the velocity profiles measured with the cross wire technique 
satisfactorily represented the vortex velocity field. Axial velocities indicated that for no 
breakdown flows the variation in the magnitude about the mean over the majority of the chord is 
± 5%, dropping by about 15% near the apex region. Both the magnitudes and locations of the 
structures, including the shear layer region and core location, are well represented and compare 
favorably with previously obtained data. Hence, the data could be manipulated further to 
examine other derived properties of the flow, such as vorticity and circulation. 


5.2.2 Vorticity 


The velocity field data was centrally differenced spatially to obtain the axial vorticity 
component: 


_ oby 3v 
i2x “ 57 ' dz 


(5.1) 


The resulting vorticity field for the A = 75° a = 20° data of Figure 5.6a is contoured in Figure 
5.21a. A repeat test appears in Figure 5.21b for comparison. The majority of the axial vorticity 
is concentrated in the region immediately around the core of the the primary vortex. The extent 
of this region is approximately y/s = 0.55 to 0.75 and z/s = 0.3 to 0.5. Outside of this regime, 
the flow presents itself as essentially free of the axial vorticity component. A smaller region of 
flow with vorticity of the opposite sign and a much lower magnitude is located in the vicinity of 
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the secondary vortex. The recorded maximum and minimum values differ by about 7% and 14% 
respectively. The calculated axial vorticity for the seven hole probe data of Figure 5.6b is 
presented in Figure 5.22. Substantially smaller magnitudes of for both the positive and negative 
maximums, over 30% and 50% respectively, are evident when compared to the cross wire data. 
This is suspected to be a direct result of the increase in the grid resolution of the cross wire data, 
which is 30% finer than that used by the SHP. As previously emphasized in Table 1 of Chapter 
2, the aspect of grid resolution should be strongly addressed when calculation of point properties 
in flow regimes containing large gradients is performed. 


A = 75° a = 20° x/c = 0.5 



Figure 5.22. Axial Vorticity at A=75°, a=20°, x/c=0.5 SHP Data 
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The most obvious difference is in the form of the contours themselves. The cross wire data 
displays a significantly more ragged appearance away from the core region than the SHP data. 
The nature of this jaggedness is believed to arise from two sources. Despite the repeatability of 
the velocity data, the error associated with maintaining the exact same spatial locations increases 
with each of the four successive cross wire passes required for a complete survey. Since the 
vorticity depends on this spatial derivative, a slight discrepancy in the actual location will be 
magnified in the resulting vorticity values. An estimate of the uncertainty in the vorticity due to 
the error in the spatial measurement depends on the chordwise station measured. At x/c = 0.3 
this uncertainty is less than 1 % using the resolution of the measuring system as the maximum 
error. Aft of this location, the uncertainty decreases further. Note that this uncertainty is based 
solely on spatial error and does not account for the error assiciated with the values of the 
velocity. Appendix C contains further details. The second source of error stems from the 
sensitivity of the technique to direction in regions of lower speeds. This seems to correspond 
with the contour shapes, for as distance from the core regions is increased, the contours become 
more distorted.. 


A hypothesis presented in Chapter 2 was that the vortex structure may be unable to exist in 

the cohesive, pre-breakdown nature if a maximum local value of vorticity in the vortex, say 
^ x max’ * s reac hed. The usual convention is to nondimensionalize the calculated vorticity by the 

ratio of the root chord to the freestream velocity to provide an overall view of the absolute 
vorticity in the entire flowfield. As mentioned previously, the use of the local spanwise distance 
provides a means to scale each cross sectional flow plane in order to examine the effect of the 
local geometry on the flow characteristics. It the flow scales in this manner, it demonstrates a 
conical behavior. 


The maximum value of vorticity was obtained at each station as well as the minimum. 
These values are presented in Figure 5.23 for the 75° sweep wing at 20° and 30° angle of attack. 


A 



0 

* A 

% 

\ O 

Angle of Attack 
20° 30° 

o maximum a 


d. 

A 

• minimum A 


* o 0 

r - - _ 

o “ 



o * 

6 * ‘ 



• - 



/ • 




Trend of Breakdown Free Row 


Chordwise Location (x/c) 


A A A 
A 


Angle of Attack 

20° 

30° 

o 

maximum a 

• 

minimum A 


O 


. 


• r • 


— , Trend of Breakdown Free Row 


Chordwise Location (x/c) 

(b) 

Figure 5.23. Chordwise Variation of Maximum and Minimum Vorticity 
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Note that breakdown is occurring over the wing at the angle of attack of 30° between the 0.50 
and 0.55 x/c location. No breakdown occurred over the wing for the 20° angle of attack case. 
The data is nondimensionalized using the two conventions mentioned above, the local semispan 
and the chord. From Figure 5.23a, in which the vorticity is multiplied by the ratio of the root 
chord to the ffeestream velocity, the magnitudes of both the positive and negative maximum 
values are seen to decrease with increasing distance from the apex. The positive values of 
vorticity in Figure 5.23a show a slight drop in magnitude at this station, however the negative 
values indicate no change at all. If the vorticity is now nondimensionalized by the local semispan, 
s*, as shown in Figure 5.23b, the data indicates a constant value behavior in the chordwise 
direction for both the positive and negative values. The previously noted drop in the positive 
value for the 30° case is now more evident. 

The scatter in the data can be attributed to several factors. The vorticity is a measure of the 
smallest scales of the flow. Details of the velocity gradients will be lost if the grid is not 
sufficiently fine. The grid size increment for the data in Figure 5.23 is 3 % of the local semispan. 
Finer grids were taken for certain configurations and are plotted in Figure 5.24. Data acquired 
on the larger chord halfspan model is included as well. The magnitude of the positive axial 
vorticity increases with increasing grid resolution, as might be expected since smaller flow scales 
can be resolved. The effect on the maximum negative magnitude was not of the same order, 
however, possibly due to the small size of the secondary vortex. It should be kept in mind that 
these values of vorticity were obtained by differentiating discrete data which has the effect of 
increasing the error associated with such measurements. As the grid sizing becomes smaller, 
though, the relative error in each spatial measurement will increase the scatter of the derived 
vorticity, since the absolute error of the traversing system is fixed. Hence, simply examining the 
maximum value of the vorticity does not make a definitive statement of the status of the vortex at 


that chord station. 
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The leading edge sheet continually feeds vorticity into the vortex flowfield along the entire 
length of the delta wing. In conjunction with the earlier proposal of a maximum vorticity criteria, 
was the hypothesis that the vortex can only exist in the pre-breakdown state if the amount of 
vorticity is below a certain level . The axial vorticity profiles throughout the vortex were 
therefore examined for possible maximum amounts of this distributed vorticity at a chordwise 
location. Integrating the positive and negative values of the axial vorticity separately over their 
respective areas leads to the results shown in Figure 5.25. 

Figure 5.25 contains data for several different configurations and methods of data 
acquisition. It is important to keep in mind that only the 75° sweep planform at oc = 20° is 
maintaining a breakdown free flow above the wing. Both of the other planforms have a 
breakdown region slightly aft of the x/c = 0.5 station. From Figure 5.25a, the amount of 
negative vorticity present above the wing is observed to be relatively similar for all three 
configurations, the A = 70°, a = 20° case having the largest magnitude at each respective chord 
location. For all three planforms, the integrated positive values show a marked increase with 
distance from the apex. The magnitude of the 70° sweep data at 20° angle of attack is seen to be 
approximately the same as the values associated with the 75° sweep at an angle of attack of 30°. 
The 70° wing will generate a stronger vortex than the 75° wing if both are at the same angle of 
attack, but since the 75° wing is at a =30° while the 70° is at 20°, the strength of the 75° wing 
vortex is increased. Scaling the data by s* causes the positively integrated values to remain 
relatively constant or slightly decrease with increasing chord location. In addition, the A = 70 , 
a = 20° case is now seen to fall in the same region as A = 75°, cx = 20°. The negative vorticity 
values remain together at approximately a constant value for all three cases. 
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Figure 5.25. Integrated Vorticity a) Scaled by chord b) Scaled by local semispan 
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In an effort to examine the validity of the theoretical ideas proposed by Brown and Lopez 
(1988) values of azimuthal and axial vorticities, Q 0 and Q x respectively, were examined. Recall 
that the basic premise of Brown and Lopez was that a negative azimuthal vorticity 
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will induce a negative axial velocity and initiate the breakdown process. The calculation of £2 0 
and Q x was based on the grid traverse through the core center utilizing polar coordinates. 
Under the assumption that the radial velocity and its gradient were negligible compared to the 
other terms along the traverse through the core: 


and 


Q 0 = - 


dV x _ du 
dr - dr 



w 3w 

7 + dr 


(5.3) 


(5.4) 


The nondimensional calculated values of Q x and Q 0 are given in Figures 5.26a and 5.26b. 
The axial voracity is seen to increase to a maximum at the center of the core. A slight rise also 
occurs at r/s = 0.30, probably due to the shear layer in the feeding sheet. The earlier maximum, 
presented in Figure 5.21 and based on the sparial derivative in Cartesian coordinates, shows a 
comparable value. Values were also derived from the seven hole probe data of Payne (1987). 
The SHP data follows the cross wire data probe, but does not exhibit the maximums present in 
the latter at the core center or at the leading edge shear layer, r/s = 0.3. The majority of the axial 
vorticity component is contained within 20% of the semispan, r/s = ± 0.10 about the core center 
location. 

In a similar manner, Q 0 also increases as it nears the core centerline, however it drops to a 
value of zero at r/s = 0. Since du/dr is always negative with increasing distance from the core, as 



Figure 5.26 Core Vorticity Profiles for sweep = 75°, alpha = 20°, x/c = 0.5 a) Axial Vorticity 
b) Azimuthal Vorticity c) Helicity 







165 


indicated in Figure 5.15, Q 0 is always positive except at r/s = 0 where u is a maximum and the 
gradient becomes zero. The azimuthal vorticity is seen to have a wider, more erratic distribution 
than the axial component. It also appears to have a rise at r/s = 0.30 due to the feeding sheet. 

The SHP data has also been overlaid on this figure. The scatter of both curves is believed to be 
based mostly on the lack of resolution. The most important thing to note is that for the majority 
of the profile, the values are positive. This would be in agreement with the theory of Brown and 
Lopez as these locations are pre-breakdown. The values do become negative in one region, 
however, just outboard of the leading edge, r/s = 0.3. The axial velocity suffers a deficit in this 

region, leading to a positive 0u/dr. The axial helicity profile is also included in Figure 5.26c. 

The deficit in the SHP data on the axis is due primarily to the lower vorticity magnitudes. 

If the axial vorticity profiles are now overlaid for each of the measured x/c locations, Figure 
5.27a, similar profiles are seen to exist at each station. Both A = 75°, a = 30 and A - 70 , a - 
20° cases exhibited similar features. It was hoped that some indication of the imminent 
breakdown, especially with the A = 75°, a = 30° case, would be evident from these distributions. 
The A = 75°, a = 20° vorticity data scales in the same manner as the axial velocity profiles, that is 
linearly with distance from the apex which indicates a conical behavior of the flowfield. The data 
for the A = 75°, a = 30° configuration does not give any indication of the upcoming flow 

transition, however the probe was still upstream of breakdown. 

The vorticity profiles were then nondimensionally integrated along each respective 
semispan to derive a set of values indicative of the local vorticity density distribution. These 
values are given in Figure 5.27b. As would be expected from Figure 5.27a, the A = 75°, a = 

20° no breakdown case revealed a relatively constant value in the chordwise direction. The single 
SHP point for the A = 75°, a = 20° configuration falls at the lower edge of the band of this data. 
The A = 70°, a = 20° values also fall in this range, while the data for A = 75°, a = 20° maintains 
slightly higher integrated values, possibly indicating a local maximum at the 40% station. 
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The calculated axial vorticity profiles of Iwanski LDV data presented earlier in Figure 2.3 
were also integrated in the same spanwise manner and are plotted against their chord locations in 
Figure 5.27c. These profiles were taken at stations in the breakdown region unobtainable using 
the present cross wire technique. Both the integrated vorticity and the helicity, defined as the 
product of the axial velocity and vorticity components, appear to reach some maximum at a 
distance of about 10 to 15% of the chord upstream of the breakdown region. Whether this 
indicates that the vortex reaches a saturated or critical condition is still open to question, due to 
the scarcity of data upstream of breakdown. A rapid decline through the breakdown region and 
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beyond, is then observed. The axial component of helicity, obtained by taking the dot product of 
the vorticity field with the velocity field, also appear to indicate a minimum at the 55% chord 
position. Helicity can be viewed as a correlation of the velocity and vorticity components in any 
one direction and the scaler parameter is a measure of the inhibiting process of the normal energy 
cascade from large to small scales (Davies-Jones (1982)). 

The illustrative nature of the LDV results above, in contrast to the cross wire data, 
prompted additional derivations of the azimuthal vorticity component as a means of verifying the 
conjectures of Brown and Lopez (1988). Azimuthal vorticity values derived from Iwanski's 
LDV data are presented in Figure 5.28a. This component of vorticity was observed to become 
negative at certain spanwise locations, corresponding to the beginning of the breakdown region, 
thus supporting the conjectures of Brown and Lopez (1988). The azimuthal component appears 
to reach a maximum negative value, however as this data represents time averaged values of a 
highly fluctuating region and care should be taken when drawing conclusions. A further point 
should be made about the derivation of these values. As the core expands at the onset of 
breakdown, a radial velocity component must certainly exist. There is a real possibility that this 
now becomes a non-negligible quantity and would then act to provide a positive contribution to 
£2 0 , serving to offset the negative component. This is not considered in the evaluation of 

Iwanski's data due to a lack of sufficient information. Pagan and Solignac (1986) also 
determined that the value of Q 0 becomes negative as the breakdown region is encountered. Their 
spatial results indicate local minimums shortly after the breakdown. They conclude that the entire 
vorticity vector, Q, rotates from an essentially longitudinal direction to a generally tangential one 
in the breakdown zone as Brown and Lopez's analytical hypothesis surmise. 

Physically there is no meaning to integrating the azimuthal vorticity distributions and thus 
the maximum values of Q 0 were plotted instead in Figure 5.28b. These values can also be 
interpreted as the maximum gradient of the axial velocity in the radial direction and are seen to 
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increase in the axial direction for the present A = 75° tests, especially for the cc — 30 case. The A 
= 70° reflects the same local peak about the 40% station. The data of Iwanski shows the change 
from a positive to a negative sense, as would be expected from Figure 5.28a. Unless there is a 
local axial velocity deficit present, -Ou/dr will always be positive. Other investigators have 
shown this deficit to occur in the breakdown region, but the present method was found to be 
inaccurate for regions of the flow where u is much less than v or w. The low axial velocity 
regions present at breakdown gave results which could not be considered reliable. In addition 
negative values of u could not be measured using this method. 

The behavior of the axial vorticity strongly support the conical nature of the delta wing 
vortex in the region preceding breakdown. Scaling both the maximum vorticity values and the 
profiles based on a traverse though the core by the local spanwise geometry indicates that the 
distribution of the axial vorticity through the core is similar at each chordwise station. The 
majority of this axial vorticity component is seen to be concentrated in the region immediately 
around the core of the the primary vortex. The onset of breakdown, however, cannot be simply 
characterized by an experimentally derived point vorticity value. Grid resolution and the locally 
steep gradients deter this type of quantification. A relative decline can be observed, as Figure 
5.23 suggests, but nothing absolute. The positive components of vorticity also seem much more 
sensitive than their negative counter-parts to chordwise location. The positive axial vorticity 
profiles also strongly indicate a conical flow behavior prior to breakdown. The integration of the 
vorticity profiles derived from the LDV data of Iwanski, indicates a maximum plateau region in 
the distribution of vorticity upstream of the breakdown region. Despite the scattered trends of the 
azimuthal component derived from the LDV data, the change from a positive to a negative sign, 
through the breakdown zone, appears to support the propositions ot Brown and Lopez. Since all 
the vorticity components above were arrived at using calculations based on differentiating the 
experimental data, a possible reduction in associated errors might be obtained by analyzing the 
data after it has been integrated. This then leads to the following section. 
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5.2.3 Circulation 

The circulation, T, was calculated in the plane of the survey grids and compared to the 
vorticity field over the area normal to it. The relationship between each is detailed below: 

T = j V • dr = | ( V X V ) • dA = f Q x dA (5.5) 

r A A 

In order to evaluate the reliability of differentiating the discrete velocity field, the vorticity 
field distribution was spatially integrated and is presented for the A = 75°, a = 20° case, in Figure 
5.29 along with the equivalent line integral. The values have been nondimensionalized by the 
freestream velocity and the root chord. They are plotted outward from the core center (r = 0), 
where the radial distance has been scaled by the local semispan. Each curve represents a 
chordwise location and the circulation is seen to grow in a chordwise manner. This is what one 
would expect, as more of the feeding sheet is being wrapped into the vortex as distance increases 
from the apex. The circulation increases at a decreasing rate from the center of the vortex and 
reaches a maximum inside of the planforms leading edge. The profiles experience an additional 
rise in circulation as the integration path encloses the shear layer outboard of the leading edge. 

The circulation and integrated vorticity values are seen to correspond quite well, which 
would indicate that differentiating the velocity fields did not substantially increase the error. The 
circulation was not calculated along paths corresponding to a constant radius from the core center. 
Instead, square integration paths were followed, as the data fields for the v and w velocity 
components lie tangent to these lines of integration. The enclosed path formed a square enclosing a 
circle of radius equal to the r/s location. The uncertainty in the calculated circulation values, based 
on the error in spatial measurement, was estimated to be approximately half that of the vorticty, 
resulting in less than 0.25% at x/c = 0.5. As with the estimate in the uncertainty of the vorticity, 
no account is made of the error associated with the velocity. Payne (1987) stated that the overall 






173 


positional accuracy of the system is i 1 mm. If this is taken to be the actual error in position, the 
uncertainty in the circulation rises to over 9%. The corresponding increase in the vorticity 
uncertainty is over 20%. 

The circulation values for the A = 75° planforms were then scaled by the local semispan, 
s*, and plotted in Figure 5.30, to examine the dependence on distance from the apex. This has 
the effect of bringing the curves close together. The agreement of the resulting scaled profiles 
implies that the circulation at a constant angle of attack is increasing in direct proportion to the 
increase in local semispan. This linear increase of circulation in a chordwise direction for any 
radius re-enforces the conical nature flowlield assumption. Since Figure 5.30a represents two 
angles of attack for the same planform, a possible way to account for the difference in attack 
angle would be to further incorporate some function of a as a scaling parameter. As an example, 
the circulation values for three x/c locations are divided through by sin(ot) and plotted in Figure 
5.30b. The data for A = 70°, a = 20° has also been added for comparison. The data does not 
collapse to a line, but rather a band. A universal curve of the form suggested in Figure 5.30b 
would be very useful, however it is difficult to ascertain the error present in each profile. For 
comparison, circulation profiles from the seven hole probe data of Payne (1987) were calculated 
and are given in Figure 5.31. A similar behavior is observed in Figure 5.31b when the values 
are scaled by the sine of the incidence angle and the resulting values fall in the range indicated in 
Figure 5.30b. 

A closer look at Figures 5.31 and 5.32 indicates the growth of the circulation with radius to 
follow a definite pattern. That is, there is a region of increasing magnitude in a nonlinear 
manner, followed by a leveling off, and then the subsequent increase due to the shear layer. 

Thus the flow exhibits rotational qualities outside of the subcore. As outlined in Chapter 1, 
Figure 1.7, the archetypical vortex is divided into three regions: the solid body center and the 
rotational Euler region, which comprise the core', and the irrotational flow outside of this. Each 
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Figure 5.31 Radial Variation of SHP Circulation Data Scaled by a) Local Semipan b) Angle of Attack 
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of these regions will exhibit a distinct behavior in the radial variation of the circulation and the 
velocity profiles. A solid body core typically possesses a velocity distribution of the form v = 
cor, for some radius r < a where a is the radius of the maximum tangential velocity. The 
tangential velocity profiles presented earlier in Figure 5.15 suggest a linear velocity distribution 
through the core region. The resulting circulation profile would therefore increases parabolically 
with distance from the core center. None of the circulation profiles in Figures 5.29 to 5.31, 
however, indicate this upward concave behavior. The curves all possess negative second 
derivatives with respect to r for r < leading edge. Since the subcore diameter is on the order of 
10 to 15% of the local semispan, as will be discussed later, the upwards curvature may not be 
present simply due to lack of resolution. 

As illustrated in Figure 5.15, the tangential velocity was seen to decrease with increasing 



Figure 5.32 Theoretical Radial Circulation Distributions 
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radius for r > a. An irrotational flow, of the form V 0 = 1/r, will demonstrate this behavior, but 
there will be no accompanying increase of circulation with radius as illustrated in Figure 5.32. 
Two other theoretical circulation profiles derived trom velocity distributions based on an inverse 
relation of the radius were plotted in Figure 5.32 for comparison. An increase in the circulation 
similar to that in Figure 5.29 is seen for one of the distributions, possibly indicating a 1/r" 
velocity distribution for n <1 exists at radial distances greater than the point of maximum 
tangential velocity. 

Further insight on the nature of this rotational region can be gained from the work of 
Hoffmann and Joubert (1963) on turbulent line vortices. Based on their analysis of the equations 
of motion, utilizing the assumption that the 'inertia' terms are negligible compared to the 
Reynolds stress terms, a circulation distribution of the form below was derived 



where A is a constant. This is also plotted in Figure 5.32 with the corresponding velocity profile 
equation. 

The data in Figure 5.30a and 5.31a was replotted in 5.33a and 5.33b using a logarithmic 
scale for the radius. The cross wire data, in Figure 5.33a, shows evidence of the suggested 
logarithmic dependence with radius for a turbulent vortex before the shear layer flow is 
encountered. Hoffmann and Joubert (1963) also obtained data at a series of stations downstream 
a vortex generated by a differential airfoil. They suggest a universal circulation distribution 
exists of the form 


p 

= 2.14 logjo (~) + 1 

r(a) a 


(5.7) 
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Figure 5.33 Dependence of Circulation on Logarithmic Radial Position a) Present Tests b) SHP data 
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This suggests a constant slope which is not apparent in the present cross wire data. Figure 5.33a 
indicates the slope for the planform at a = 30° to be greater than at a = 20°. The slope of the 
curves follows a fit of 2sin(a)/sin(A) which would lead to a circulation relation of the form 


r 

Ha) 


2 sin(a) 
T(a) sin(A) 


In(^) 


+ 1 


(5.8) 


Interestingly, a consistent logarithmic behavior is not evident from the SHP data as shown 
in Figure 5.33b, the most deviation occurring at the larger radii. The reason for the apparent 
discrepancy is unknown. Referring back to the linear scale plot in Figure 5.31, it can be seen 
that the radial increase of circulation shows less curvature than the cross wire profiles, even to 
the point of being linear with radius. 

The previous circulation profiles in Figures 5.29 to 5.31 illustrate the distribution of the 
circulation in the radial direction. As noted, the profile magnitudes increase with radial distance 
to approximately r/s = 0.25, representing the extent of the primary vortex. A further increase in 
magnitude with increasing radius corresponds to contributions from the shear layer at the leading 
edge. The maximum value of the circulation, scaled by the chord, is plotted for the A = 75°, a = 
20° planform at each chord station in Figure 5.34a to examine the axial variation of the vortex 
strength. Although there is some scatter in the data, a linear trend corresponding to the maximum 
radius of integration, r/s = 0.6 is to be seen . This may be attributable to the location of the lower 
boundary of the integration path, that being the upper surface of the wing for r/s distances greater 
than the z/s location of the vortex axis. The exact proximity of the probe to the wing surface will 
result in measurements of the spanwise velocity components that could vary in intensity. This 
would reduce the magnitude of the circulation in a different manner at each x/c station, leading to 
a the variation in the maximum values. 
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Figure 5.34 Chordwise Variation of Circulation Scaled by a) Chord b) Local Semispan 
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A better representation of the behavior of the vortex is obtained if the circulation values at 
the radial location where the profile initially levels off are plotted. The values for this radial 
distance of r/s = 0.25 are also plotted in Figure 5.34a. A much smaller scatter is present in the 
data. The values tend to follow a near linear distribution, except near the aft of the wing surface, 
whereupon a leveling off is observed. This further supports the arguments pertaining to the 
conical behavior of the delta flowfield, since a conical flowfield demonstrates the linear growth 
of circulation in a chordwise direction, as noted by McCune and Tavares (1988). 

Scaling the values in Figure 5.34a by s* instead of the chord produces the distributions of 
Figure 5.34b. Both the r/s = 0.6 and r/s = 0.25 cases indicate values which are dropping slightly 
in the axial direction. The SHP data falls below the cross wire for both the r/s = 0.6 and the 0.25 
radii. This may be due to the resolution increase of the cross wire data. 

The data from the other wing configurations is presented in Figure 5.35a, along with that 
of Figure 5.34a, for r/s = 0.25 . The increase in vortex strength due to angle of attack is clear 
for the 75° sweep case, as one would expect. The A = 70°, a = 20° planform also generates a 

larger magnitude in circulation than that of the 75° sweep wing for the same angle of attack. 
Hence, as is already well documented in the literature, increasing the angle of attack for a given 
wing geometry, or decreasing the sweep angle for a planform at a given angle of incidence 
causes breakdown to occur closer to the apex. Both of these results therefore support the 
argument that, all things being equal, a stronger vortex will breakdown earlier, that is closer to 
the apex. One can surmise from this that there must be a limit to the strength, or amount of 
vorticity, beyond which the vortex can not sustain the pre-breakdown state. 

An interesting observation is made, however, of the data if s* is again used as the scaling 
parameter in Figure 5.35b. The data tor A = 70° is seen to come much closer to, or even below, 
that of A = 75° for a = 20°This would imply that there is a similar, or slightly lower, total 
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amount of distributed vorticity over the 70° wing per unit span at any chordwise station than the 
75° wing, despite the larger absolute values. Since breakdown occurs on the 70° sweep wing at 
about the x/c = 0.5 chord station and the A = 75° for a = 20° planform does experience 
breakdown at all, it could be construed that knowing the local strength of the vortex is not 
enough to quantify the vortex state with regards to breakdown at any given chordwise station. 

To further illustrate the chordwise behavior of the circulation, consider the SHP data 
obtained by Payne (1987) above a A = 85° for a = 40° planform at several stations directly in the 
breakdown region. Derived circulation values at r/s = 0.25 and 0.7 are presented in Figure 5.36. 

The circulation values are seen to fall through the breakdown region. Since the measured 
circulation is directly related to the axial vorticity component, either the axial vorticity has been 
redistributed or it has become less. Both reasons are correct possibilities. Brown and Lopez 



Figure 5.36 Chordwise Dependence of SHP Derived Circulation for A = 85° and a =40° 
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(1988) postulated that the predominantly axial vorticity vector rotates at breakdown, away from 
the axial direction to a predominantly tangential direction. This would explain the drop in the 
values of calculated circulation. Alternatively, the expansion of the core causes the vorticity to be 
spread over a larger area and the integration path would therefore result in lower circulation 
values. The total circulation would continue to rise, however, as distance increased from the 
apex, since vorticity is still being added to the vortex structure though the shear layer. 

Hence, it can be concluded that at breakdown the expansion of the core must be 
accompanied by an redistribution of circulation over the entire vortex, possibly in a combination 
of expansion and turning of the vorticity vector. This supports a hypothesis that it is not just the 
total amount of circulation present which determines whether the vortex breaks down or not. It 
is the concentration of that circulation which gives rise to the breakdown. Just as vorticity is 
integrated over the area it acts to give a value of the circulation, the circulation can 
be multiplied over the area it encompasses to give an indication of the density of the circulation. 
This would also serve to explain the nature of the circulation values discussed with respect to 
Figure 5.35b above. 

The range of angles of incidence and sweep angles investigated by Payne allows the 
dependence of circulation on these variables to be examined. The derived values are presented in 
Figure 5.37. All the circulation values correspond to a chordwise location of x/c = 0.5. The 
circulation at every configuration grows in a linear manner with angle of attack. As the 
breakdown moves ahead of the x/c = 0.5 station for the higher angles of attack, the circulation 
values drop, in the same manner as was pointed out in Figure 5.36. This again seems to support 
the Brown and Lopez (1988) argument the the longitudinal, or axial vorticity component loses 
some intensity to the tangential direction. For the three lowest sweeps the circulation is also seen 
to decrease in some constant fashion with an increase in sweep. This becomes more evident in 
Figure 5.37b when the data is scaled by the local span. 
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Figure 5.37 Dependence of Circulation on Sweep and Angle of Attack 
a) Scaled by chord b) Scaled by local semispan 
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The behavior of the circulation, as examined in the preceding discussion, supports the axial 
vorticity and velocity distributions trends in indicating that a near conical flowfield exists for the 
the majority of the pre-breakdown leading edge vortex. The flow scales with the local geometry, 
and to a lesser extent with the sine of the angle of incidence. The radial circulation profiles for 
the present tests exhibit a logarithmic dependence on radius. The strength of the vortex increases 
in a linear manner with distance from the apex and with angle of attack for a fixed chord location. 
This increase in strength is inevitably followed by breakdown which reduces the value of 
circulation in the v-w plane, or axial vorticity density, about the vortex axis. Breakdown cannot 
eliminate the circulation present and thus the onset of breakdown cannot be solely attributable to 
the total vortex strength in absolute terms, for the total circulation ol the vortex is still increasing 
with x/c, even after breakdown . Instead it is believed that the local circulation taken about the 
core region plays a most significant role in the onset of breakdown. There is a sharp decrease in 
this value after breakdown, presumably in a non-reversible manner for the flow never transitions 
to the pre-breakdown state. The local strength of the vortex is not believed to be the sole initiator 
of breakdown. The circulation is determined only from the velocity flowfield and no indication is 
given about the pressure forces present. No account is taken of the the axial flow component.. 
To this end, the investigation turned to the correlation of the flow properties from a parameter 
perspective. The most widely used of these is the swirl angle. 


5.2.4 Swirl Angles and Other Correlation Parameters 

The swirl angle defined as the inverse tangent of the tangential velocity component to the 
axial component. The magnitude of this angle reaches a maximum of approximately 42° to 50 in 
the flow directly preceding breakdown. Values computed from the current tests based on 
traversals through the vortex core are plotted in Figure 5.38 alongside of those derived from the 
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SHP data. The cross wire data gave larger swirl angle values by about 5° both inboard and 
outboard of the vortex core when compared with the SHP data. Even larger angles were 
recorded outboard of the leading edge, reaching as high as 52° in the 70° sweep shear layer. The 
reason for these discrepancies lie in the measurement of the axial component as noted in Figure 
5.10. The SHP data gave velocities which were larger in magnitude, resulting in lower swirl 
angle values than the cross wire data. 

The change in the swirl angle with distance from the apex is displayed for all three of the 
present test planforms in Figure 5.39. The maximum values are seen to change little for each 
wing. The 20° angle of attack cases in Figures 5.39a and 5.39b are observed to maintain swirl 
angles of about 40° despite the difference in sweep. The increase in attack angle to 30° in Figure 
5.39c increased the maximum swirl by another 10°. This change indicates that the swirl angle is 
sensitive to the flow, for at x/c of 0.5 the 75° sweep wing at a = 30° is much closer to breaking 
down than at a =20°. Yet the 70° sweep wing is also closer to breaking down and gives no 
indication of a larger swirl angle. Since swirl angle is only a measure the ratio of rotational 
velocity to the axial convective velocity, it does not account for any local shearing in the fluid. 
The swirl angle gives an indication of the direction of the flow and of how it is being turned 
locally, but does not indicate whether there is a change in the local fluid rotational rate. 


A more physically descriptive parameter can be determined using the both local vorticity 
and velocity vectors. Brown and Lopez (1988) use the ratio of these two vectors in their study 
as an precursor to breakdown. Since only the axial vorticity component is available in the 
present study, the swirl parameter derivation outlined in Chapter 2 from the tornado studies, 
which is an indication of the convection of vorticity, is examined here. The parameter is simply 
the ratio of the local axial vorticity to the local axial velocity as defined by 




i in Vorex Core Swirl Angle a) sweep = 75 alpha = 20 
alpha = 20 c) sweep = 75 alpha = 30 
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Computed fields for the 75° sweep wing at a =20° and x/c = 0.4 and 0.7 are given in 
Figure 5.40. Both stations show comparable values and area of concentration. Three areas of 
intensity are noted: the core of the primary vortex, the secondary vortex and the shear layer. The 
vortex core reveals large magnitudes ot this swirl parameter, despite being scaled by the locally 
large axial velocities, because the majority of the positive vorticity is concentrated there. The 
secondary vortex region also indicates a local concentration of this parameter, albeit in the 
negative sense. Both the secondary vortex and the shear layer region are areas where the axial 
velocity component is less than freestream and this serves to increase the local convection 
parameter. The maximum and minimum values at each station are directly comparable. 

The SHP data for this wing is presented in Figure 5.41a. Although there is a qualitative 
comparison in terms of the associated regions of intensity, the maximum and minimum values 
are lower by as much as 40% and 30% respectively. The grid resolution is again believed to be 
the primary cause of this. The areas do seem to be ol a comparable size, however, indicating a 
comparison may be possible in terms of integrated values. The SHP data at a = 30 is also 
included as Figure 5.41b. The maximum and minimum values are slightly higher in magnitude. 
It also appears that the region of the negative convection values, bounded by the zero contour line 
has increased in size, relative to the a = 20° case. 

The maximum and minimum convection parameter values for the present tests are plotted in 
Figure 5.42. If the chord is used as the length scale, the values follow a behavior similar to the 
vorticity values in Figure 5.23. The positive and negative magnitudes are decreasing with 
distance from the apex. It could be argued, however, that the presence of the probe in the 
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smaller vortex structure near the apex, while not causing breakdown, may serve to locally 
accelerate the flow artificially and measure erroneous gradients. This could result in larger 
vorticity values. The A = 75° a = 20° planform data marked with an asterisk indicates data taken 
at a higher grid resolution. A single asterisk stands for a y/s = z/s increment of 0.015 while two 
asterisks denote an increment of 0.010. These increased resolution values indicate a similar 
range of convection values. The finest grid results are up to 100% greater than the nearest larger 
grid values. This directly results from the observed vorticity values in Figure 5.24. 

The convection parameter was integrated over the area of the survey, as was done to the 
vorticity for comparison with the circulation. When scaled by the chord, Figure 5.43a, the 
integrated values for the survey plane decrease in a chordwise direction, quite steeply. The 
relative state at each station, indicated in Figure 5.43b, show a similar distribution at each x/c 
location. It might be expected that the integrated values would show a relative decrease as 
breakdown was approached for the A = 70° a = 20° and the A = 75° a = 30° configurations near 
x/c = 0.5. No indication of this is apparent. Seven hole probe data is include for comparison. 
The values fall below their cross wire counterparts, presumably due to resolution. 


As a final approach to the type of convection considerations discussed, a Rossby number 
was calculated for the above flows to examine its significance, if any. The parameter was 
defined as 

Average s ubcore axial velocit y . (5.10' 

— Average subcore axial vorticity * radius of the subcore 

Figure 5.44 illustrates the values of this parameter for the present study. The majority of the 
values fall to either side of R 0 = 1 . Repeat tests for the 75° planform at 20° angle of attack and 
x/c = 0.5 indicate a substantial variation. The proposed nature of this parameter is to decrease as 
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the breakdown region is encountered. Very slight evidence of this is indicated from the A 75 
a = 20° or the A = 70° a = 30° tests. The finer grids, marked with an asterisk, seem to approach 
a value of 1.1. The SHP data is considerably larger in magnitude, again owing to the lower 
vorticity values. Overall the results of this parameter evaluation with respect to the onset of 
breakdown, and that of the convection parameter above, are inconclusive. 


The series of parameters evaluated above were based on the local and integrated properties 
of the flowfield. In a sense, these can be regarded as conditions that result from changes in the 
planform configuration or external flowfield. It is also possible to set up parameters to 
incorporate these external conditions, such the sweep angle and angle of attack, along with the 



measured flowfield properties. This has already been done, to a certain extent, by scaling the 
field properties by the local semispan and even the sine ot the angle of attack. As outlined in 




Chapter 1, this line of reasoning was incorporated by Hemsch and Luckring (1990) in their 
correlation results for the measured vortex circulation. Their derivation, as a function of the 
Sychev parameter and expressed in terms of the apex half angle e is: 
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where 


g = 


U c tan^e cosa 


= AK n 


K = 


tan a _ 1 
tan e ^1 


(5.11) 


for some value n. Hemsch and Luckring noted that if g and K are plotted in a log-log format, a 
fit of the form g = AK^ “ was seen for data obtained from Wentz and MacMahon (1967), A = 
62° and Delery, Pagan, and Solignac (1987). Values of g an K ranged from 0.5 to 10.0 and 
0.2 to 2.0 respectively. Theory by Smith (1971) proposed a relation of the form g = 4.63K n . 
The g and K parameters were derived from the seven hole probe data of Payne (1987) and 
plotted in Figure 5.45a with that of Wentz and MacMahon as well as a curve of Smith for n = 

1.2. Payne's data represents sweep angles of 70°, 75°, 80°, and 85° at various chord locations. 

It is seen to extend the line of Smith to a g of 100 and a K of 10. 

The data from this investigation, representing three K values, was also plotted in Figure 
5.45a, however the expected collapse of data was not immediately evident. Hemsch and 
Luckring used data that was acquired in the wake of the models. Since the present tests were 
conducted at locations above the wing surface, a further scaling of g by the local chord ratio, x/c, 
was found to bring the data into line with that of Smith as shown in Figure 5.45b. The data of 
Payne was also scaled by x/c and is shown in Figure 5.45b. Thus it would appear that this 
relation strongly correlates the vortex strength with the angle of attack and the wing geometry. 







It is interesting to see how closely this corresponds to the function used in Figure 5.31b. If n is 
set equal to 1 in the above expression, they would be exactly identical. 

The derivation of Hemsch and Lucking correlates the data quite well, however, the __ 

expression g = f(K) is such that the independent variables of a and £ appear on both sides of the 
equation. Physically, the circulation is not inversely proportional to some function of the semi — 
apex angle, but increases as £ increases. Perhaps more representative evaluation of the data can 
be done by simply plotting the dependent versus the independent variables in the manner: 

= f(K’) (5.13) 

Note that the local chord is now employed instead of the local semispan as done in the previous _ 
section, removing the implicit dependence on A from the left hand side of the equation. This 
relation using K' is plotted in Figure 5.46 for a sample of the present data. The data behaves in a 
linear manner with what appears to be a larger spread in the data than was evident in the 
correlauon of Hemsch and Lucking. If the axis of the plot are changed to accommodate a larger 


r 

U M c (x/c) 


= f(tanatan£) - f 
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linear maximum or are replotted in a logarithmic format, as in Figure 5.46b, a trend similar to 
that in Figure 5.44 can be observed. The linear nature of this parameter is further demonstrated 
by scaling the circulation by K' and shown in Figure 5.47. The data for the 75° planform at 20° 
angle of attack indicates a downward trend. If breakdown was imminent, the curve would drop 
off, as would be expected from the SHP data in Figure 5.36. No indication of this trend is seen 
for the other two test configurations. 

A relation of the form examined by Hemsch and Lucking appears to provide one of the best 
correlations of the present data and that in the literature over a wide range of test conditions for 
conditions upstream of the breakdown region. The behavior of the convection parameters 
correlate the flow upstream of the breakdown as well, but seem to be more sensitive to the nature 
of the flow than the Hemsch and Lucking relationship, perhaps because of the logarithmic 



Figure 5.47 Chordwise Dependence of Scaled Circulation Parameter 
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format. No concrete evidence of an impending breakdown is present in terms of a distinct 
chordwise change in the convection parameter representation or the Hemsch and Lucking 
relationship. This would indicate that the breakdown phenomena is very localized spatially and 
does not transmit substantial time averaged information upstream. Since the change in any type 
of these time averaged parameters in a streamwise manner is not indicative of breakdown, it 
appears for the moment that data must be acquired directly in the breakdown region, as was done 
with the SHP, to directly observe any direct change on any of these relations. It may be that to 
measure some indication of the onset of breakdown in the upstream flow a time dependent signal 
correlation will be required. 

5.2.5 Core Spectral Behavior 

Several chordwise stations on the 75° sweep planform, at 20 and 30 degrees angle of 
attack, were examined to determine it there was any characteristic frequencies, or possible 
standing waves, associated with the axial velocity component. A single hot wire was positioned 
in the core and the voltage signature was analyzed in terms of its frequency content given by a 
power spectrum distribution. The following spectral data are ensemble averages of twenty 
spectra as the power spectrum routine was constrained to a maximum of 1024 points. The 
sampling frequency was 6000 Hz as higher rates indicated no dominant frequencies above the 

3000 Hz range. 

Before the spectrums are examined, a typical pre-breakdown signal is given in Figure 5.48, 
followed by a signature of the post breakdown flowfield. The passage of more turbulent flow 
structures is evident in Figure 5.48b. The resulting power spectrum in Figure 5.49a, however, 
did not reveal any dominant frequency for the post breakdown signals. The post breakdown 
spectrum indicates a much greater magnitude of the entire spectrum at the lower frequencies than 
the x/c = 0.4 pre- breakdown. The x/c = 0.4 spectrum is magnified in Figure 5.49b and reveals a 





Amplitude 
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small peak at about the 300 Hz range. The x/c = 0.6 spectrum rises sharply in magnitude as zero 
Hz is approached. If enlarged there appears to be a peak between 5 and 10 Hz. Since the 
frequency resolution for these tests was only about 6.1 Hz, however, the peak could actually be 
an indication of low frequency noise. 

This peaking of the spectrum between 250 to 500 Hz appeared to be characteristic of all the 
pre-breakdown signals. The magnitude ot the signal was greater at 30° than 20° for the x/c = 0.3 
location as seen in Figure 5.50a although the profiles had similar shapes. As the Reynolds 
number was reduced, the difference in magnitude of the signals was seen to increase. 

At x/c = 0.5, in Figure 5.50b, close to the breakdown location for a = 30°, the power 
spectra were seen to differ considerably from Figure 5.50a. There is a resemblance to the x/c = 
0.6 profile of Figure 5.49a in terms of shape for a = 30. It is possible that the a = 30° spectrum 
in Figure 5.50b represents the flow in the breakdown region. The magnitude of the fluctuations 
is seen to have dropped considerably compared to the spectrum in Figure 5.50a. The a = 20° 
spectrum is seen to exhibit a larger magnitude in Figure 5.50b compared to Figure 5.50a. 


The peak present in the power spectrum was also observed to grow in magnitude as x/c 
was increased. Spectra at chordwise stations ot x/c = 0.3, 0.5, 0.7, and 0.9 are presented in 
Figure 5.51 for an angle of attack of 20 degrees. An abrupt change in the spectrum profile from 
x/c = 0.7 to 0.9 is evident. The transition resembles the pre to post breakdown transition noted 
in Figure 5.49a. The magnitude is considerably lower at x/c = 0.9 than 0.7, similar to Figure 
5.50b. Whether this is a result of probe induced breakdown or of proximity to the trailing edge 
of the wing is unknown. 

The freestream flow was tested with no model in the tunnel at the various Reynolds 
numbers to observe any frequency contributions from the wind tunnel. The values presented in 
Figure 5.52 lie three orders of magnitude below the amplitudes shown previous. There is a peak 


Amplitude 


Re = 250000 sweep =75 x/c = 0.3 


a = 20 
a = 30 


1.0e-5 • 


0.0e+0 


1000 1500 2000 2500 3000 

Frequency (Hz) 


Re = 250000 sweep = 75 x/c = 0.5 


a = 20 
a= 30 


o 500 1000 1500 2000 2500 3000 

Frequency (Hz) 
b) 

Figure 5.50 Effect of Angle of Attack on Power Spectrum a) x/c = 0.3 b) x/c = 0.5 









208 


between 5 and 10 Hz, but this maximum value is smaller than any of the previous values by a 
factor of about 100. On the otherhand it may be that the flow field above the delta wing is 
sensitive to these low frequency components and that they serve to create the peaks evident on 
the lower end of the spectrums in Figures 5.49b and 5.50b. 


5.2.6 Core Dimensions 


Flow visualization has indicated another interesting phenomena that would substantiate 
arguments of a cylindrical behavior of the subcore region and a conical nature of the outer portion 
of the core. Although the outer flow region of the vortex scales with the local geometry of the 
wing, the subcore maintains what appears to be a constant diameter as seen in Figure 2.5. If 
titanium tetrachloride smoke is introduced into the core. Figures 2.5a, the traced diameter does 
not widen at the same rate as the outer flow in the axial direction. In what resembles an inverse 
image of the above in Figure 2.5b, the void produced in the core region when kerosene smoke is 
introduced upstream of the model also remains constant in diameter with an increase of distance 
from the apex. The nature of the velocity profile data leads to the possibility of two core size 
definitions: a jet core due to the presence of the axial velocity component and a much smaller core 
based on the distance between the maximum and minimum values of w or Vp, usually referred to 
as the subcore. If a minimum value of u/U,^ is used to set a threshold cutoff, a jet core can be 

defined and the growth rate compared to that of the subcore defined above. Verhaagen has 
defined a third type of core based on the vorticity profile. This rotational core is based on the 
radius to which the vorticity falls to a certain level, in the same manner as the jet core outlined 
above. 

The core diameters evaluated from the A = 75° a = 20° data are presented in Figure 5.53. 
The jet core is defined as the diameter of the axial velocity profile above a threshold of u/U 00 = 
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1.5. This was felt to be the minimum u/U^ cutoff, since the core tends to widen appreciably as 
u/U^ is decreased further. Both the local semispan and absolute units are given in Figure 5.53a 

and 5.53b respectively. The grid resolution of 0.03 y/s leads to more scatter in the determination 
of the size of the subcore than the jet core. Finer grid resolution data corresponding to y/s — 
0.015 and 0.010 are also indicated in Figure 5.53a and 5.53b by the solid black triangles. 

Figure 5.53b indicates that both the jet core and subcore are growing in size, the jet core at a 
slightly larger rate than the subcore. When scaled by the local semispan in Figure 5.53a, the 
relative jet core diameter is seen to decrease in the axial direction as the trailing edge is 
approached. This is confirmed by Figure 5.53b where the absolute size of the jetcore is leveling 
off near the trailing edge. The finer grids also indicate the subcore to be increasing slightly in 
size towards the trailing edge, although a constant diameter can be interpreted over the central 
portion of the wing. The scatter of the data makes interpretation difficult and more data will be 
required at a finer grid resolution before any definitive statements can be made. 

Core size evaluated from the SHP data tor the same geometry is plotted in Figure 5.53 
along with the cross wire data and falls in the same range. Five hole probe data from Verhaagen 
and van Ransbeeck (1990) is also included in Figure 5.53b. They measured a subcore diameter 
from x/c = 0.1 to 0.7 on a 2.22 m, 76° swept delta wing. The constant diameter would indicate a 
cylindrical rather than conical flow. Data obtained from Hawk, Barnett, and O Neil (1990), 
however, points to a linear type of core growth with chordwise distance over a 762 mm chord, 
70° swept wing. These discrepancies again indicate that due to such factors as grid resolution, 
which is usually not compared between experimental data sets, differing statements concerning 
observable properties can exist between investigators and comparisons should be examined in 
light of these factors. 


The core sizes calculated for the other planforms of the present study are presented in 
Figure 5.54 along with the SHP equivalents. Both techniques indicate similar diameters for both 
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the jet and sub core. The A = 75° a = 20° jet core is substantially larger than the 70° planform. 
The subcore exhibits a constant value over the chordwise distance shown. Both planforms in 
Figure 5.54 have breakdown occurring on the wing in the vicinity of the x/c = 0.55 location. 
Since there is only one point at each chordwise station, and this data is taken over the central 
portion of the wing, the trend of the data in the chordwise direction is difficult to ascertain. 

The data above would seem to indicate that the the subcore is growing in diameter in a 
downstream direction, not at the rate of the local wing geometry, perhaps, but growing 
nonetheless. The flow visualization in Chapter 2 suggests that the growth rate, if any, is small 
and not easily seen visually. This indicates that care must be taken when making quantitative 
judgements from flow visualization. On the other hand, the scatter in the data and the indication 
that at times there does appear to be a constant diameter for some x/c range points to the need for 
a more precise measurement technique before coming to a strong conclusion on this subject. In 
addition, there is a probe interfering with the flow when the measurements were taken which is 
not present the flow visualization. The effect of this probe, which was on the order of the 
subcore diameter, has the tendency to displace the core are the result of this on the core diameter 
has not been quantified. 

As mentioned at the beginning of this section, a core diameter can also be defined according 
to the region containing the rotational flow. Based on this definition, Verhaagen and van 
Ransbeeck (1990) report the rotational core to occupy approximately 35% of the local semi span. 
From the cross wire data in Figure 5. 27, based on a threshold vorticity value equal to 90% of the 
maximum, a rotational diameter equal to 30% of the local semispan is measured. If the 
rotational core definition is to contain all the rotational region, that is until the axial vorticity drops 
to zero, 90% of the local semispan would have to be included. 
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5.3 Pressure Measurements 

The effect of pressure on the behavior of the vortex structure and the breakdown location is 
very substantial as has been pointed out by several investigators in Chapter 1 and discussed in 
Chapter 2. Payne (1987) measured a change in the slope of the total pressure coefficient on the 
vortex axis versus chordwise position curve from negative to positive at the 60% chord mark for a 
85° swept wing at a 45 degree angle of attack. This is then followed by a significant rise in the 
value of Cpt. These measurements were taken through the breakdown region and agree physically 
with what would accompany the local stagnation of the axial velocity component. 

It was therefore decided to observe the changes in the pressure distributions on and above a 
delta wing to examine if the pressure distribution gave an indication that breakdown was evident. 
Whether this came as a result of some maximum value of the local surface pressure coefficient or a 
critical change in the pressure gradient, dp/d\, this could then be incorporated into a parameter for 
the breakdown prediction. 


5.3.1 Surface Pressures 

The tests conducted with the surface pressure models centered on three factors: the 
influence of the planform size on the measured pressures, the effect of the sting mount location, 
and whether or not it was possible to trace the path ol the breakdown by observing the pressure 
signature on the wing surface. 

The values of the surface pressure coefficient, Cp, along a ray extending from the apex, 
gave no obvious indication that a breakdown existed at some chord location over the wing. The 
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data along the y/s =0.6 location exhibited the lowest values of Cp, indicating they were closet to 
being under the core, and thus will be the only values which will be discussed here. 

The 16 inch model Cp values are presented in Figure 5.55a. The value of dp/dx is always 
seen to be positive. Increasing the angle of attack only serves to increase the nonlinearity of the 
curve, but no indication of where the breakdown position could be determined. This is also 
observed from the data presented by Roos and Kegelman (1990) although it was not explicitly 
stated by them. This would seem to indicate that the use of this variable in developing a 
breakdown parameter would be minimal. No plateau is observed and no drastic sign change in 
the slope occurs. Kirkpatrick (1970) noted a definite step type of change in the value of Cn 
verses angle of attack as the breakdown moved forward on a 68° wing, if the spanwise pressures 
are integrated. Roos and Kegelman comment that although the lift aft of breakdown decreases 
with incidence, the overall lift continues to increase until the flow is fully separated from the 
wing. In any event, simply observing the surface pressure signatures is not enough to determine 
a breakdown location which was the intent here. 

There was an interesting change discovered in the pressure measurements when they were 
repeated on the 16” chord wing. This second set of data, presented along with the initial data in 
Figure 5.55b, was taken with the wing centered in the tunnel. The first set of tests were 
conducted with the wing much closer to the tunnel floor and severe differences appear at angles 
of attack of 35° and greater . The close proximity to the wall appeared to delay the stall of the 
wing in the first series of tests until over 45 degrees, as compared to the typical angle of attack 
around 35 degrees. The second set of data shows and indication of the onset of stall at 35° at the 
x/c = 0.7 location. The a = 40° data shows a rise in the pressure at all stations on the centered 

wing as compared to the wing near the tunnel floor. The latter was then seen to experience a 
similar pressure rise once an angle of attack of 50° was reached. 
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The effect of increasing and decreasing the chord length relative to the 16 inch model is 

shown in Figure 5.56. Each plot represents an angle of attack ranging from 25 to 45 degrees. 

Both 16 inch cases are included. The 12 and 20 inch models have the sting mounted in the same 

absolute location to simulate a reduction and extension of the trailing edge respectively, relative to 

the 16 inch model. At angles of attack below 40°, the models show similar profiles, although the 

20" model exhibits a consistently lower pressure field. At 35° the 12" model indicated a 
considerable pressure rise which would corresponds to a stall according to the Q a curve. The 

centered 16" model also exhibited the verge of this transition as mentioned previous. The 20" 
model indicates a still lower pressure profile, which continues to drop as the angle of attack is 
increased to 45° in Figure 5.56d. The profile of the centered 16" model approaches that of the 
12" case by a = 45°, but no indication of stalling is present on the 20" model. 

The effect of the sting location is illustrated in Figure 5.57. On both the 12 and 20 inch 
models, the sting was first positioned at the same absolute location as the 16" model, that being 8 
inches from the apex. The mount was then changed to the local x/c = 0.5 mark which is closer to 
the apex of the 12" by 2" and close to the 20" model trailing edge by the same amount. Overall 
the resulting pressure profile differences are minimal. The 12” model shows slightly lower 
values Cp for the x/c = 0.5 location, while the 20" model gave slightly higher pressures. This 
could be stated such that for both models, movement of the sting rearward decreased the 
measured Cp values. The largest difference, of approximately 9.5%, occurred at a = 30° for the 
12" model. In general the 12" and 16" centered models correlate quite well, while the 20" model 
was found to have consistently lower pressure values. 

The surfaces pressure distributions gave no indication of the onset of breakdown when 
viewed along rays extending from the apex. It was not until total separation occurred over the 
wing that there could be a noticeable difference in the OP/dx behavior at the surface. Because of 
the lack of a definitive trend in this measured parameter, it is ineffective as an indicator of the 
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onset of breakdown and would contribute little if included in a correlation parameter. It is 
apparent that the blockage effects involved with the measured wing configurations are not 
negligible. This factors are already known, however, and blockage correction procedures are 
available. More importantly because of its subtlety, and perhaps not as well known, is that the 
mounted position of the model in the tunnel is critical to representing proper flow conditions. 


5.3.2 Total Pressure and Centripetal Acceleration in the Core 

Total pressure measurements were made, via the total pressure probe, above the 75 sweep 
wing with the primary intent of examining the ideas put forth in section 2.2. The intent was to 
compare the balance of the radial pressure gradient with the centripetal acceleration of the fluid. 
The data obtained from the total pressure probe could be combined with the velocity data of the 
hot wire surveys to evaluate these parameters. The total pressure data is presented in Figure 
5.58a in the form of total pressure coefficient, Cpt , against chordwise location. Two angles of 
attack were tested, 20° and 30° Previous flow visualization with the hot wire on the a = 20° 
tests showed no breakdown occurring, either naturally or induced by the hot wire probe. The a 
= 30° case indicated that breakdown was induced by the hot wire probe to occur at about 55% of 
the chord, where it would naturally occur in the wake. 


From Figure 5.58a, it can be seen that the pressure measurements indicate breakdown is 
not occurring on the 30° case until perhaps at about the x/c = 0.8 station. There is no probe 
induced breakdown at the x/c = .55 location. The total pressure probe is considerably smaller 
than the hotwire configuration as seen in Figures 4.4 and 4.5. Whether breakdown was actually 
occurring at the 80% station was not observed visually, however it was certainly not occurring 
where it does with the x-wire probe. Thus any comparison, such as that in section 2.2, must be 
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made forward of the breakdown induced by the hotwire probe. This serves to indicate how 
sensitive the vortex flow is to the intrusion of a probe and most especially the size of the probe. 


The 20° case reveals an increase in the the total core pressure with distance, up to about x/c 
= 0.6, followed by a subsequent decrease, but careful interpretation of the data is required. As 
the probe is moved rearward on the delta wing, the local flow increases in size accordingly and 
the probe size becomes relatively smaller when compared to the local vortex station. Thus it may 
be that an area encompassed by the probe tip as a single point measurement further near the apex, 
may now be resolved into several points. In essence the probe can now measure a lower 
pressure more accurately, if it exists, since it is not being smeared into a single point. This 
implies that the actual core pressures nearer the apex may be lower than measured. In other 
words, if the pressure is decreasing in the axial direction, the relative probe size may be the 
cause, whereas if it is increasing as x/c increases, the phenomena would not be a result of the 
probe resolution. Conversely, the oc = 30° measurements show a rise in the core pressure, after 
an initial drop, up to the x/c = 0.8 station and it can be surmised that the pressure is actually 
rising. Error bars represent the maximum error associated with the largest measure Cpt values . 

Recall, from Chapter 1, that Lambourne and Bryer (1961) noted that the total pressure was 
seen to have a nearly constant value along the length of the vortex at roughly Cpt = -5.0 for a 65° 
sweep wing at a = 15. Kegelman and Roos (1990) report a similar constant total pressure with 

axial distance for 60° and 70° flat plate delta wings. Naarding and Verhaagen (1988) reported a 
drop in the axial total pressure with increasing x/c as well as a sharp drop near the apex for their 
76° sweep wing at ot = 20. Their data is given in Figure 5.58b for comparison and is seen to 
have Cpt values lower than the present values by up to 100%. They explain the difference 
between their data and that of Lambourne and Bryer is because Naarding and Verhaagen used a 
probe size on the order of the viscous core size. The present tests used a probe on the order of 
1/3 to 1/4 the size of the subcore. The better resolution of Lambourne and Bryer would explain 



223 


their lower pressure values, but does not explain the reduction of values as Naarding and 
Verhaagen move closer to the apex. Smearing ot the pressure profile due to a locally large probe 
size, would result in higher pressures, not lower. Kegelman and Roos (1990) also point out that 
the area where this significant total pressure loss occurs becomes proportionately smaller with 
increasing x/c. This could indicate that the subcore region, where the majority of the pressure 
loss occurs, is maintaining an constant absolute size. 

Lamboume and Bryer mention that three factors on which the pressure distribution along the 
axis depends. The first of these is that increasing vortex strength tends to provide a falling axial 
static pressure. One might also propose that a lower axial pressure or a greater radial pressure 
gradient would allow for the existence ot a stronger vortex. In either case, the strength of the 
vortex would seem to be directly dependent on the local radial pressure gradient, however it is 
established or influenced by external factors. Theretore the proposed ratio of forces was 
introduced at the end of Chapter 2 and will now be examined in more detail. 


Four ratios were suggested to characterize the vortex state at a given axial station. The first 
two were based on the ratio of the difference in pressure between the core axis and the 
ffeestream, either static or total, to the maximum centripetal acceleration measured at the edge of 
the subcore. The latter two required the local pressure at the edge of the subcore. Unfortunately 
it is not possible to obtain these two latter ratios, even with the total pressure data available from 
the SHP surveys. If the total pressure probe is not aligned directly with the local flow vector, 
the total pressure is not being measured. One might argue that, given the velocity vector, the 
dynamic components could be added to the measured pressure at the required location to obtain 
the total pressure. This would be true in a field which did not have any losses. Since the 
leading edge vortex indicates a drop in the total pressure at the axis relative to the external flow, 
it is impossible to know the total pressure at a given radius without actually measuring it. 
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The two pressure/acceleration ratio, PAR, terms were therefore calculated in the following 
manner based on the ratio: 
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The ratios are defined explicitly as: 
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and is obtained from the hotwire data at that station. 

An example of the SHP total pressure data at the core axis and and the associated velocity 
profiles is given in Figure 5.59. The maximum value of the tangential velocity occurs at the edge 
of the subcore, indicated by the two solid vertical lines. The pressure coefficient used is the 
maximum recorded and the axial component required for ratio 1 was also taken as the maximum. 
The calculated ratio 1 and ratio 2 PAR values are plotted in Figure 5.60a against their respective 
angle of attack. Since all the stations were taken at an x/c location of 0.5, changes in the axial 
direction are unobtainable. 


225 


The ratios based on the static pressure distribution, ratio 1, are up to four times greater than 
the ratio 2 values based on the total pressure. The concept of the PAR itself implies that a value 
greater than one indicates the pressure gradient is sufficient to keep the vortex together. In this 
respect all the calculated values support this for they were taken in the pre-breakdown flow. The 
ratio 2 values are closer to a value of 1 and might be thought of as to the more representative of 
the flow. It is the local static pressure, however, that determines the local fluid forces and this is 
the more proper term to consider in this evaluation. As stated previously, a more indicative ratio 



Figure 5.59 SHP Total Core Pressure and Velocity Profiles A = 75°, a = 20° x/c = 0.5 


would be obtained by knowing the actual pressure gradient at the radius of the maximum 
tangential velocity, but this is unavailable. 
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Seven hole probe data was available for a chordwise evaluation of the 85° sweep wing at 
45° angle of incidence. The computed values are presented in Figure 5.60b. Interestingly, the 
profiles peak in the region just before breakdown. Since the core expands, it would be expected 
that the ratio would drop below one, rather than increase. After breakdown, the vortex is again 
in a state of balance and a ratio of one or greater would be expected. Close examination of the 
data revealed that at the station where the peak occurs in Figure 5.60b, the tangential velocity 
had dropped, but the minimum Cpt did not indicate any rise until the next x/c station. Care must 
be taken with these values since they represent time average quantities in a fluctuating field. 

Time averaging may have smoothed out the instantaneous drops in the core pressure which 
would send the PAR to below 1. One would expect this ratio to be greater than 1 in the post 
breakdown flow as well, for it too is a stable rotating vortex. 

Lastly, the data from the present tests and the corresponding similar geometry planform 
results from other investigators is presented in Figure 5.61. The majority of the ratio 2 PAR 
values are slightly above 1. Values derived from the five hole probe data of Verhaagen and van 
Ransbeeck (1990) registers larger values of the ratios due to their larger magnitude Cpt values. 
The cross wire data shows similar trends tor both the ratios, a relatively constant value, followed 
by a sharp increase after x/c = 0.7. The ot = 30° cross wire data indicates a similar trend. 

No strong trends are readily evident from the results of this pressure acceleration ratio. The 
values are all seen to be greater than unity which w'ould be expected in the stable vortex flow, but 
the expected drop towards one as breakdown approached was not observed. More data is 
required, especially near the breakdown zone, and instantaneous velocities and pressures are 
required for the breakdown region itself. In addition the pressure and velocity should be 
recorded simultaneously. The unsteady nature of the flow in the breakdown region precludes 
any evaluation as to the instantaneous, local state of the flow, based on time averaged 


measurement techniques. 
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Figure 5.61 Pressure Acceleration Ratio at A = 75° a - 20 

Although the above correlations do not give strong indications of the impending breakdown 
downstream of their position, this is not to say that they should be disregarded in evaluating the 
flowfield. The data used to derive these correlations is based on averaged quantities which were 
obtained at different times. In the broadest sense, they do indicate that the flow is in a stable 
condition and not about to breakdown, which is true since the data was acquired in a region of 
the flow upstream of the breakdown region. More data is required to properly evaluate the 

effectiveness of such parameters. 

The results of the present study have led to several conclusions, mostly dealing with the 
vortex flowfield upstream of the breakdown region. The velocity fields have been compared and 
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discussed. The values of vorticity and circulation, as well as their spatial distributions, are good 
indicators of the flow variation in the chordwise direction. The breakdown region is seen to alter 
these values and distributions in a significant way. Dimensional scaling has been demonstrated to 
provide a means for correlating the vortex behavior . The conclusions of this study are 
summarized in a more complete form in the following chapter. In addition, several comments are 
made concerning the role of the vorticity and circulation in creating the flow conditions necessary 
for breakdown. 


CHAPTER VI 


CONCLUSIONS 


The preceding discussion strongly supports the belief held by many investigators, both 
experimental and theoretical, that the flow over a delta wing, upstream of the breakdown regions 
and away from the apex and trailing edge regions, behaves in a conical manner. That is, properties 
remain constant along rays emanating from the apex. Dimensional scaling of the radial circulation 
distribution by the local semispan indicates the flow properties to be similar at each chordwise 
location for regions of pre-breakdown vortex flow. Spanwise vorticity and velocity distributions 
based on a single traverse through the core of the vortex also scale with the local geometry in the 
pre-breakdown state. This is a good indication that the local semispan is a representative length 
scale which should be included with any nondimensionalization scheme used to examine flows 
generated by swept leading edge geometries. 

Further evidence of the conicity of the flow field is supplied by the nearly linear increase in 
strength of the vortex, at a given radius from the core axis, with distance from the apex. For a 
fixed chordwise station, the dependence of circulation on the angle of incidence was also found to 
increase linearly. The relation of the form examined by Hemsch and Luckmg appears to provide 
one of the best correlations of the present data, as weU as that found in the literature, over a wide 
range of test conditions. The correlation parameters indicate the sensitivity of the flow to the 
slenderness of the planform and incorporates both the effects of sweep and angle of attack. 
Inclusion of the relative chord location, x/c, as shown in the present study, accounts for regions 
where the vortex size and strength are increasing, such as over a delta wing. 
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This increase in strength is inevitably followed by breakdown which was seen to reduce the 
circulation about the vortex axis. The local concentration of axial vorticity is reduced while the 
azimuthal vorticity deceases throughout the breakdown zone. This supports the theory of Brown 
and Lopez (1988), who postulate that the physical mechanisms involved in breakdown rely on the 
production of negative azimuthal vorticity resulting from a tilting and stretching of the 
predominantly axial vorticity vector. 

The nature of the flow to follow a conical behavior has been an underlying assumption in 
many theoretical derivations, including Stewartson and Hall (1963) and Mangier and Weber 
(1967). Some of these analysis, such as Stewartson and Hall, also employ a cylindrical subcore 
region. Flow visualization and to a certain extent the measured core diameters of the present tests, 
indicate that this cylindrical assumption is a realistic approximation. The best agreement for either 
of these geometrical flow considerations applies over the central portion of the planform, that is 
away from the apex and trailing edge regions. This was also shown to be the case with the 
experimental data of by Verhaagen and van Ransbeeck (1990). Near the apex, the geometry of the 
flat plate delta wing is such that the span of the wing is on the order of the thickness. The shape 
that the flow senses resembles that of a blunt body rather than a slender,thin planform and 
deviations from a conical flow assumption can be expected. In the same manner, the flow is 
altered as it approaches the trailing edge by the pressure recovery of the flow. The reversed flow 
regions, indicated by the surface visualization, support the presence of these effects. Overall, 
however, the majority of the flowfield properties are of a conical form, outside of the subcore 
region, and the use of this approximation in a theoretical analysis is valid. It should be kept in 
mind that this applies to the primary vortex structure. The flow near the surface is subject to 
changes depending on the nature of the boundary layer. The secondary vortex also appears to 
scale with the flow geometry, but flow visualization pictures have revealed some form of transition 
at a given chordwise location. Whether this is a form of breakdown of the secondary or not could 
not be determined with the present measurement techniques. 
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Several other conclusions were reached and are summarized below: 

a) The use of the present cross wire technique, employing four spatial passes and two 
different probe configurations for a single survey, was satisfactory in measuring the 
flow conditions and establishing trends of the pre-breakdown flow. Velocity profiles 
obtained with the cross wires compared well with previous seven hole probe and laser 
doppler anemometry data taken using the same planforms. 

b) The grid resolution at which experimental data is acquired must be accounted for when 
comparisons of data are made. This is especially so if point values are to be reported which 
result form discrete differentiation of the data. In the present study, the grid resolution of the 
cross wire data was 30% finer than the seven hole probe data and this resulted in an increase in 
the maximum vorticity of over 60%. 

c) The circulation profiles for the present tests exhibited a logarithmic dependence on radial 

distance from the core axis. The slopes of the different configurations appeared to follow the 

form 2 S * n ^ — for the planforms investigated. 
sin(A) 

d) Although the axial vorticity distribution is severely altered as the breakdown region is 
encountered, there appears to be some maximum in the spanwise vorticity distribution that 
occurs in the region just preceding breakdown. 

e) The majority of the positive axial vorticity is confined to the subcore region of the vortex in 
the pre-breakdown vortex state. The use of the maximum value of axial vorticity is deceptive in 
determining the local strength of the vortex structure, because it is a point property and 
dependent on the grid resolution. 

f) The trends in the negative axial vorticity distributions, which are concentrated in the 
secondary vortex region, remained similar for the different configurations investigated in the 
present study, despite geometric and angle of attack differences. Adjustment of the flowfield to 
changing test conditions, such as angle of attack occurred for the most part in the positive axial 
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vorticity regions. In addition, the magnitude of the negative axial vorticity was similar for all 
three test cases. 

g) Two vortex core definitions were examined using the present cross wire data: a jet core, 
based on the diameter of the vortex within which the axial velocity was greater than u/U°° = 1.5, 
and subcore, defined as the distance between the maximum and minimum tangential velocities in 
the core. In general, both cores gave indications of an increase in size with distance from the 
apex and the growths rates appeared to be scaling with the local semispan. The jet core size 
remained constant aft of approximately 70% of the chord to the trailing edge. Finer grid 
resolution data and flow visualization results suggested that the subcore maintains a constant 
diameter over the central regions of the planform as reported by Verhaagen and van Ransbeeck 
(1990). 

h) The surface pressure distribution, as measured along rays underneath the vortex core, give 
virtually no indication of the onset or passage of breakdown. 

i) The surface flow directly inboard of the leading edge and outboard of secondary attachment 
line indicates a flow direction inwards from the trailing edge, towards the wing centerline, and 
in the axial direction. If the flow is in a direction towards the leading edge, the axial component 
of the flow would need to be in a direction opposite to that of the oncoming stream. 

j) The maximum axial velocity on the 75° sweep planform at 20° angle of incidence 
mai ntains a value between u/U°° = 2.0 to 2.5 over the majority of the wing, dropping 
off in the vicinity of the apex and the trailing edge . The maximum value of the axial 
velocity was also seen to maintain a constant value for some distance upstream of the 
breakdown region for the planforms exhibiting breakdown. 

k) The power spectrum revealed a distinctive frequency peak between 250 and 500 Hz in 
all the pre-breakdown axial flow signals. 


The results of this study and others, along with theoretical discussions in the literature, 
suggest that there are two significant conditions in the flow which will cause breakdown to occur: 
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the amount of local circulation present and the distribution of the external pressure field. Both of 
these are equally important, however it is felt that the initiation of breakdown is most sensitive to 
the latter. Whether this is the local adverse gradient at the trailing edge, the gradient in the external 
flow, or the result of an object in the downstream, this is felt to be the single most important factor 
for initiating breakdown, provided the vortex is near the point of criticality in terms of strength. 

All breakdown location versus chordwise station plots are very flat as they near the trailing edge. 
The breakdown rarely maintains a position at the x/c = 1.0 position. Instead, it will move to 
approximately x/c = 0.8 an then continue upstream. It must overcome some type of bamer 1 at the 
trailing edge, which is surmised to be the pressure gradient at the trailing edge as the flow recovers 
to ffeestream at the surface of the wing. This is also one of the primary reasons for the hysteresis 
effects in breakdown location that are seen with regard to breakdown location pitching up or 
pitching down. Furthermore, this sensitivity to the pressure of the flowfield is felt to be the main 
cause for discrepancies that appear in the literature for locating the breakdown. Models are 
generally constructed with great accuracy and angles of incidence are set with ±0.1 degree in some 
instances. Yet locations can vary by as much as 25% of the chord as compiled by Kegelman and 
Roos (1989). Hall (1966) has shown that pressure gradients at the edge of the vortex core are 
amplified at the centerline as a result of the substantial swirl velocities. This causes the centerline 
flow to be very sensitive to the local environment. Hence a discrepancy in measurements could 
easily arise if the tunnel pressure gradients vary, especially in smaller tunnels. 


The second flow condition that defines the the state of the vortex, with respect to 
breakdown, is its circulation at a certain chordwise location. This was especially seen to be 
influential for tube vortices. Increasing the circulation eventually caused breakdown. This is also 
the case with respect to delta wings. An increase in angle of attack, or a sweep decrease, causes 
the vortex to strengthen and breakdown follows, proceeding toward the apex with increasing 
incidence. Yet the circulation alone can not be responsible for breakdown. The breakdown 
location occurring on a delta wing planform at a fixed angle of attack can be easily displaced 
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towards the apex by the creation of an adverse pressure gradient downstream of the model itself as 
demonstrated by Delery, Pagan and Solignac (1987). The circulation required for breakdown to 
occur has now been reduced. Thus, even the presence of the trailing edge and the accompanying 
local pressure gradient could be crucial to the location and onset of the breakdown. 

Whether the flow exists in the pre or post-breakdown state, the position of breakdown 
inevitably depends on some equilibrium with regard to the the forces in the vortex. The simplest 
explanation for the dominating factors in this balance of forces are the radial pressure gradient and 
the centripetal acceleration of the fluid. The state of the external environment which governs the 
balance of these forces determines the state of the vortex. On the pressure side, the gradient is 
determined by the pressure field outside the vortex and along the core. An external gradient in the 
chordwise direction could appreciably influence the pressure and of course the state of the vortex. 
A local rise in the pressure along the core axis will also initiate breakdown, whether this happens to 
be the effect of the trailing edge or an object placed in the path of the vortex downstream of the 
trailing edge. The geometry of the delta wing sets up the velocity field above the planform. 
Vorticity is continually being added to the vortex from the leading edge. The circulation is seen to 
increase with angle of attack, distance from the apex, and a decrease in the sweep angle, and the 
accompanying increase in the centripetal acceleration required to keep the fluid on the rotational 
path is second order in nature. Lamboume and Bryer noted in 1961 that an increase in the vortex 
strength reduced the axial pressure. The majority of correlations evaluated to date involve only the 
circulation or swirl of the flow but not the pressure forces directly. Recall that no evidence of the 
onset of breakdown is present in the Hemsch and Lucking relationship and that only the circulation 
is examined. The present attempt to derive a relationship incorporating both the centripetal 
acceleration, or strength of the vortex, and the effects of the pressure distribution only indicated 
that the vortex was stable at the measured locations, but did not indicate that some critical condition 
in the flow was imminent. Yet it is felt that this type of relationship, possibly using instantaneous 
values of velocity and pressure, would indicate that breakdown was forthcoming. 
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Many fundamental questions, arising from the simplest observations, still remain 
unanswered. Although some researchers may disagree, it has been shown that tracer particles 
introduced into the core are not centrifugally 'spun out’ as seems to be indicated by upstream 
smoke injection techniques. Neither is smoke entrained into the subcore region if introduced about 
the delta wing, except at the apex. In fact there appears to be no significant mass diffusion to or 
from the subcore region once the vortex has been established and begins to grow. Does this imply 
that all the fluid in the subcore region arrives from a streamline impinging on the apex? Although 
the Reynolds number effects are reportedly small on the overall flow, above a certain value, there 
is an effect, which has not been thoroughly quantified, on the reduction of the subcore size as 
noted by Leibovich (1984). Since this is the region of viscous flow it might make more sense to 
define a Reynolds number characteristic of the flow based on the core diameter and not some fixed 
geometry such as the wing chord. Additional questions on the flow also need to be addressed. 

Why does the vortex initially assume the intense, pre-breakdown form and not the post breakdown 
form? How is the vorticity transported into the subcore region? Is viscous diffusion at the interface 
of the viscous subcore with the outer flow the sole means for vorticity transport? Regarding the 
last question, it has been shown that the strength of the vortex is growing as it progresses along the 
delta wing, indicating a definite transport of vorticity to the core region The manner in which this 
occurs has generally been assumed to be via the shear layer. Flow visualization indicates, 
however, that streamlines emanating from the leading edge shear layer never spiral directly into the 
core, but follow trajectories about the along a cone projecting from the apex. 

Furthermore, what really occurs along a streamline path that impinges on the apex, and is 
subsequently redirected into the core, to cause so great a loss in the total pressure? A simple 
evaluation of the drop in total pressure associated with a solid body rotation yields A Cpt - - 
(V0/Uoo) 2 . For a V0 of 1.5 the A Cp t = -2.25. This is not enough to account for the measured 
drop in the pressure data. The drop in the total pressure implies substantial losses somewhere. 
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The apex seems to be the only logical location. If this is so, the importance of the geometry of the 
apex takes on a greater significance. If it is possible to tailor such geometry to the point where the 
drop in the total pressure is maximized, the static pressure in the core will drop accordingly, 
enabling larger lift. 

With respect to further investigations, several suggestions are offered. The first pertains to 
the measurement procedures. The number of spatial passes required for a measurement plane must 
be reduced, preferably to one. Although techniques such as a three component LDV are desirable, 
they have their own problems and intrusive probes have indicated the ability to make comparable 
measurements. The size of the probe is crucial and efforts should be aimed at the smallest probe 
possible that can measure a completely unknown velocity vector. To incorporate these constraints 
may even require the design of a new form of velocity detection measurement probe. 

Instantaneous measurements appear to be essential in understanding completely what is 
occurring in the pre and post-breakdown flow. The flow upstream of breakdown exhibits a quasi- 
steady nature and many investigators have used time averaged measurement techniques as a valid 
approach to understanding the flow. The breakdown phenomena does not appear to transmit 
substantial time averaged information upstream, at least based on the parameter evaluations 
considered here. The oscillatory nature of the breakdown zone remains unresolved. Conditional 
sampling and spatial correlation techniques could be used to determine the reasons for this 
fluctuation in breakdown location. The indication that the radial circulation distribution follows a 
logarithmic behavior points the way to an evaluation of the Reynolds stress field. Instantaneous 
pressures as well as velocities are a must to properly examine the concept of pressure acceleration 
ratios . 


Secondly, comparative experiments should us as big a model as is possible under the 
circumstances to increase the resolution of the probe itself. The half span measurements of 
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Verhaagen and van Ransbeeck (1990) use a model 2.22 m in length compared to the present test 
size of under one half meter. Yet they report a subcore diameter of under 9 mm, on the same order 
as the present tests. 

Finally, the condition of the streamwise pressure gradient bears further study in light of the 
possible sensitivity of this flow to this condition. A study using artificially imposed gradients 
would indicate the extent of influence that this parameter has on the vortex flow and help to 
establish bounds for actual test conditions. 



APPENDIX A 


DELTA WING VORTEX DATA 


Anders, K., ” LDV Measurements of the Velocity Field of a Leading Edge Vortex Over a Delta 
Wing Before and After Vortex Breakdown" Von Karman Institute for Fluid Dynamics, 
Technical Note 142, March 1982 

- wind tunnel. Re = 170,000, Uoo = 18 - 19.5 m/s 

- flat plate, 137.5 mm chord, 1.5 mm thick, 8.5° stnd bevel 

- aspect ratio = 1.6, A = 68.2; 

- LDA core traversals: 

a = 19.3 ; x/c = 0.67, 0.74, 0.82, 0.93, 1.0 
a = 28.9 ; x/c = 0.55, 0.65, 0.73, 0.88, 1.0 


Anderson, M.W., Beran, P.S., and McCann, M.K. "Vortex Breakdown Over Delta Wings” 
Graduate Aeronautical Laboratories, California Institute of Technology, Pasedena, CA 91125, 
1983 

- water tunnel, dye injection, Uoo = 0.5 ft/s 

- flat plate, 7.8 - 15" chord, 1/16 - 1/4" thick, 26.6 - 33.7° stnd bevel 

- A = 64.8, 75, 80, 85 ; a = 20, 25, 30, 35, 40, 42, 45 


Carcaillet, R., Manie, F., Pagan, D. and Solignac, J.L. "Leading Edge Vortex Flow Over a 75 
Degree-Sweep Delta Wing- Experimental and Computational Results." ICAS 86. 

- water tunnel, dye injection, Uoo = 40 m/s, Rec = 4e06 

- flat plate, 0.5 m, 1.45 m chord, 20° stnd bevel 

- A = 75 ; a = 20 

- LDV, five hole probe, surface pressure measurements 

Chigier, Norman A., "Measurement of Vortex Breakdown Over a Delta Wing Using a Laser 
Anemometer" NEAR TR 62, Nielsen Engineering and Research Inc., Mountain View, Ca. 
June 5, 1974. 

- wind tunnel, Re = unknown, Uoo = unknown 

- flat plate, 137.5 mm chord, 1.5 mm thick, symmetrical bevel 

- aspect ratio = 1.46, A = 70; a = 20, 25, 35 

- LDA core traversals: 

x/c = 0.71, 0.789, 0.868, 1.507 

Eamshaw, P.B. and Lawford, J.A. "Low Speed Wind Tunnel Experiments on a Series of 
Sharp Edged Delta Wings" ARC R & M, No 3424, March 1964. 

- wind tunnel. Re = 200,000 - 400,000, Uoo = 80 ft/s 

- plano/convex, 0.589 - 1.178 ft chord. 
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- A = 45, 55, 60, 65, 70, 76 ; a =12- 37 

- force, surface flow vis. 

Elle, B. J. "An Investigation at Low Speed of the Flow near the Apex of Thin Delta Wings with 
Sharp Leading Edges " ARC R & M, No 3176, January 1958. 

- water tunnel, air bubbles, Re = 700,000, Uoo = 12 ft/s 

- flat plate, 1/16" thick, 7° stnd bevel 

- A = 70 half span ; a = varied 


Elsenaar, A., Hjelmberg, L., BUtefisch, K., and Bannik, W.J. "The International Vortex How 
Experiment" 1988 in Portugal 

- water tunnel, air bubbles, Re = 700,000, Uoo = 12 ft/s 

- flat cropped delta w/ NACA 64A005 profile , 1/16" thick, 7° stnd bevel 

- A = 55,65 full span ; a = varied, M= 0.4 to 4.0 

- LDV.five hole probe, force, surface pressure measurements; 

surface flow vis. 


Erickson, Gary E„ "Vortex How Correlation ".Technical Report AFWAL-TR-80-3143 Hight 
Dynamics Laboratory, Wright Patterson Air Force Base, January 1981 

- water and wind tunnels, various Re, A, 

- many different configurations 

Hummel, D., and Srinivasan, P.S. "Vortex Breakdown Effects on the Low-speed Aerodynamic 
Characteristics of Slender Delta Wings in Symmetrical How" Journal Qf the Royal Aeronautical 
Society Vol.71 April. 

- wind tunnel. Re = 1.4e06 - 1.7e06, Uoo = 40m/s 

- flat upper and beveled lower surfaces, triangular cross section 
0.52 and 0.625 m chord 

- A = 60, 68.2 ; a = varied 


Hummel D. "On the Vortex Formation Over a Slender Wing at Large Angles of Incidence 
AGARD CP-247 January 1979, pp.15-1 - 15-7. 

- wind tunnel, Re = 2e06 . , 

- flat upper and beveled lower surfaces, triangular cross section, 750 mm chord 

- aspect ratio = 1 .0, A = 76; a = 20.5 

- surface pressure measurements 

- Five Hole Probe grid surveys: 

in 4 planes behind the wing trailing edge 


Iwanski, Kenneth P., "An Investigation of the Vortex How over a Delta Wing With and 
Without External Jet Blowing" Masters Thesis, University of Notre Dame, April 1988. 


- wind tunnel. Re = 150,000 

- flat plate, 15.85" chord, 1.0" thick, 45° stnd bevel 


- A = 70; a = 30 

- LDA core traversals: 

x/c = 0.15, 0.225, 0.3, 0.375, 0.411, 


0.448, 0.484, 0.521, 0.557, 0.593 
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Kirkpatrick, D.L. "Analysis of the static Pressure Distribution on a Delta Wing in Subsonic 
How" ARC R&M 3619 1970 

- wind tunnel, U°° = 120 ft/s. Re = 3.15e06 

- symmetrical X-section, max t/c = 0.048, 50" chord 

- Aspect Ratio 1.616, A = 68; a = -2° - 26° 

- Surface Pressures 

Kegelman, J.T. and Roos, F.W. "The Howfield of Bursting Vortices over Moderately Swept 
Delta Wings" AIAA Paper 90-0599 Jan 8-11, 1990 Reno, NA. 

Kegelman, J.T. and Roos, F.W. "An Experimental Investigation of Sweep Angle 
Influence on Delta Wing Hows" AIAA Paper 90-0383 Jan 8-1 1, 1990 Reno, NA. 

- wind tunnel. Re = 0.3e06 - 2.0 e06, Uoo = 10 m/s 

- A =60, a = 12 - 20; 70, a = 25 - 33 

- flat plates, 24" (60°) and 30" (75°), 0.5 " thick, 25° stnd bevel 

- LDA, SHP, surface pressures, flow visualization 

Kjelgaard, S. and Sellers, W.L. IQ, "Detailed Howfield Measurements over a 75° Delta Wing 
For Code Validation" NASA Langley Research Center, Hampton, Virginia NASA TN 2997 

- wind tunnel, Re = 0.5e06, 1.0e06, 1.5e06,Uoo = 30 ft/s 

- flat plate, 22.392" chord, 0.3" thick, 10° stnd bevel 

- A = 75; a = 20.5; 

- LDV, surface vis, total pressure, five hole probe, 

x/c = 0.7, 0.9, 1 . 1 

Lamboume, N. C., and Bryer, D. W., "The Bursting of Leading Edge Vortices-Some 
Observations and Discussion of the Phenomena," Aeronautical Research Council. Reports and 
Memoranda. No. 3282, 1962. 

- wind and water tunnel, Re = 10,000 - 4.6e06, Uoo = 80 ft/s 

- flat and cambered plates, 3.1, 8, 8.5, 47.5 " chord, 16° bevel 

- A = 50,65; a =16 -28 

McKeman, J.F. and Nelson R.C. "An Investigation of the Breakdown of the Leading Edge 
Vortices on a Delta Wing at High Angles of Attack" Masters Thesis, University of Notre Dame, 
January 1983. 

- wind tunnel. Re = 225,000 

- flat plate, 16" chord, 0.75" thick 

- A = 70 ; a = 20, 25, 30, 35, 40 ; P = 0 - 12 in 1° increments 

Naarding S.H.J. and Verhaagen, N. G. "Experimental and Numerical Investigation of the 
Vortex How Over a Sharp Edged Delta Wing with and without Sideslip" Delft University 
Report LR-573, December 1988. 

- wind tunnel. Re = 2.5e06, U^, = 44 m/s 

- biconvex, 850 mm chord 

- aspect ratio =1.0, A = 75.96 ; a = 21.1 ; p = -20 - 20 

- force, surface pressure, five Hole Probe, Laser light sheet, surface vis 
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Pagan, D. and Solignac, J.L. "Experimental Study of the Breakdown of a Vortex Generated By 
a Delta Wing" La Recherche Aerospatiale, No 3, May-June 1986. 

- wind tunnel. Re = 580,000, Uoo = 14.5 m/s 

- flat plate, 560 mm chord 

- A = 75; a = 19.3 

- LDA and Five Hole Probe in the far wake and breakdown in wake 

Payne, F. M. "The Structure of Leading Edge Vortex Flows Including Vortex Breakdown" 

PhD Dissertation, University of Notre Dame, May 1987. 

- wind tunnel, Re = 255, 000, Uoo = 30 ft/s 

- flat plate, 16" chord, 0.25" thick, 25° stnd bevel 

- a = 10, 20, 30, 40 

- LDA and SHP core traversals: 

A = 70, 85, x/c =0.5 a = 10, 20, 30, 40 
A = 85, a = 40, x/c = 0.55, 0.6, 0.65, 0.675, 0.7, 0.75 

- SHP grid surveys (normal to model): 

A = 70, 75, 80, 85, x/c = 0.5 a = 10, 20, 30, 40 
A = 85, a = 40, x/c = 0.55, 0.6, 0.65, 0.675, 0.7, 0.75 

Peckham, D.H. "Low-Speed Wind Tunnel Tests on a Series of Uncambered Slender Pointed 
Wings with Sharp Edges" ARC R&M 3186 N62-10253, Report Aero 2613, December 1958. 

- wind tunnel. Re = 2.3eo6 - 8.6 e06,Uoo = 80 - 300 ft/s 

- flat plates, biconvex, thick wings 100" chord, 0.25" thick, 25° stnd bevel 

- A = 76, gothic ; a = 10, 20, 30, 40 

- pressure, force balance, surface flow visualization 

Schrader, K. F., Reynolds, G. A., and Novak, C. J. "Effects of Mach Number and Reynolds 
Number on Leading-Edge Vortices at High Angle of Attack AIAA Paper 88-0122 January 11- 
14, 1988. 

- wind tunnel. Re = 0.25e06 - 10 e06,Uoo = 80-300 ft/s 

- flat plate, sharp edge, 20.3 cm chord 

- A = 63; a = 0 to 33 

- LDV, force balance, surface pressure, surface flow visualization 
Sforza Data?? 

Thompson, D.H. "A Water Tunnel Study of Vortex Breakdown Over Wings with Highly 
Swept Leading Edges" Australian Defence Scientific Service Note ARL/A. 356 May 1975 

- water tunnel, dye injection, Re = 9,800,Uoo = 74 mm/s 

- flat plate, 150 mm chord, 1.8 mm thick, symmetrical 30° included bevel 

- A = 60, 65, 70, 75, 80; a = prebreakdown - 45 

Verhaagen, N.G.. "An Experimental Investigation of the vortex Flow Over Delta and Double 
Delta Wings at Low Speed" Delft University of Technology Report LR-372 Sept 1983, Delft, 
The Netherlands AGARD-CP-342 Aerodynamics of Vortical Type Flows in Three Dimensions, 
April 1983. 

- wind tunnel. Re = 1 .4e06 ,Uoo = 30 m/s 
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- flat plate, sharp edge, 20.3 cm chord 

- A = 76; 76/60, 76/40 double delta wings, a = 5 to 25 

- force balance, , surface oil flow visualization, laser light sheet 

Vorropoulos,G. and Wendt, J.F. ."Preliminary Results of LDV surveys in the Compressible 
Leading Edge Vortex of a Delta Wing" Von Karman Institute for Fluid Dynamics, Technical 
Note 137, August 1982 

- wind tunnel. Moo = 0.4, 0.62, 0.80 

- flat plate, 150 mm chord 

- aspect ratio = 2.0, A = 63.4; a = 10 

- LDA grid surveys: 

x/c = 0.8, 0.6, 1.0 

Wentz, W.H. "Wind Tunnel Investigations of Vortex Breakdown on Slender Sharp Edged 
Wings" PhD Thesis, University of Kansas, 1968 

- wind tunnel. Re = le06 

- flat plate, 13 models, 15° stnd bevel 

- A = 50 - 85 ; a = varied 
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FullSpan sweep: 


75 

406.4 mm 


Ay*: 
Az* : 


0.03 

ao3 


alpha 

X* 

filename 

u/U- 

v/U- 

w/U- 

Mu u/U Core 



Min 

Mu 

Min 

Mu 

Min 

Mu 


@ z* 

20 

0.20 

O2Z05.dat 

056 

2.05 

-1.14 

1.38 

1.38 

1.21 

072 

038 

20 

030 

2004 .dat 

035 

2.27 

-1.02 

1.48 

-1.36 

1.47 

066 

036 

20 

030 

1 1,05-dal 

038 

2.25 

-084 

1.42 

-1.24 

1.17 

066 

036 

20 

0.40 

12/04.dat 

052 

2.33 

-093 

1.48 

-1.25 

1.28 

069 

037 

20 

041 

13/04. dal 

041 

2.28 

090 

1.39 

-1.27 

1.25 

066 

037 

20 

042 

14/04 .dal 

031 

2.30 

-091 

1.40 

-1.25 

1.23 

066 

040 

20 

043 

16/04. dal 

046 

2.29 

-088 

1.48 

-1.34 

1.29 

066 

036 

20 

043 

04 Z05.dat 

050 

2.26 

-1.09 

1.44 

-1.20 

1.62 

066 

036 

20 

043 

07/OS.dat 

041 

2.29 

-081 

1.44 

-1.23 

1.19 

066 

036 

20 

044 

17/04.4x1 

039 

2.29 

-1.04 

1.40 

-1.22 

1.21 

066 


20 

0.45 

18/04 .dal 

044 

2.42 

-081 

1.45 

-MS 

1.36 

066 

036 

20 

047 

23/04 .dat 

051 

Z47 

-080 

1.51 

1.33 

1.25 

066 

036 

20 

047 

17/05.dat 

052 

146 

-091 

1.41 

-1.10 

1.17 

066 

036 

20 

048 

03/05,d*t 

038 

135 

-096 

1.44 

-1.36 

1.30 

066 

039 

20 

049 

2324-dat 

045 

135 

-084 

1.41 

-1.16 

1.20 

066 

036 

20 

050 

11/04 .dat 

037 

133 

-092 

1.38 

-1.25 

1.19 

069 

040 

20 

050 

19/04. dal 

042 

142 

-1.08 

1.46 

-1.43 

1.37 

069 

039 

20 

060 

24Z04.dat 

046 

130 

-085 

1.43 

-1.19 

1.25 

066 

035 

20 

0.63 

23/05.dal 

044 

150 

-082 

1.37 

-1.18 

1.32 

066 

037 

20 

067 

24/05.dal 

039 

141 

-098 

1.54 

1.25 

1.43 

066 

037 

20 

070 

3 0/04. dal 

045 

134 

-080 

1.34 

-1.18 

1.26 

066 

037 

20 

070 

21Z05.dat 

039 

132 

-091 

1.43 

-1.51 

1.38 

066 

037 

20 

075 

25Z05.dat 

028 

133 

-088 

1.34 

-1.29 

1.25 

066 

037 

20 

080 

Ol/DSdai 

040 

136 

-079 

1.34 

-1.02 

1.19 

066 

037 

20 

0.90 

29A)5.dal 

024 

126 

-098 

1.23 

-MO 

1.19 

063 

037 

20 

1.10 

30(05. dat 

029 

110 

-075 

1.27 

-087 

1.53 

060 

036 

30 

030 

l4/05.dat 

039 

3.38 

-1.22 

2.07 

-1.84 

1.93 

066 

042 

30 

040 

14/05.dat2 

044 

3.49 

- 1.11 

2.12 

-201 

1.81 

066 

043 

30 

050 

15/0 5.dat 

031 

3.20 

-1.17 

200 

-1.87 

1.83 

063 

045 

30 

055 

01/06.dal 

025 

182 

-1.05 

1.10 

• 1.11 

1.18 

066 

045 

30 

060 

27/OS.dal 

023 

132 

-1.80 

230 

-1.92 

1.77 

060 

047 

30 

051 

1 106. dal 

0.00 

3.22 

-1.16 

278 

-1.91 

3.26 

063 

041 

30 

049 

18/06 

034 

187 

-1.27 

1.96 

1.83 

231 

060 

039 

FullSpan 

sweep : 

75 
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250000 


Ay*: 

0015 

y/r 



chord : 

406.4 mm 




Az* : 

0015 

i * : 
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x* 

filename 

u/U- 

v 

/U- 

w/U- 

Mu u/U Cora 



Mm 

Mu 
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Mu 
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Mu 

@y* 

<3> i* 

20 

040 

2208. dal 

« 

. 

. 






20 

040 

2308. dal 

1.03 
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-096 

1.56 

-1 JO 

M3 
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039 

20 

050 

240Rdal 

1.08 

154 

-0.93 

1.64 

-1.25 

US 

070 

039 

20 

060 

2808. dal 

1.02 
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-1.05 

1.49 

-133 

134 

069 

037 

20 

050 

2908.4a! 

098 

151 

-0.94 

1.49 

-135 

1.40 

070 

039 

20 

070 

3008.dat 

098 

154 

-0.90 

1.36 

-1.19 

1.25 

066 

037 


y*: 0 to 1.2 

z* : surface k> 1 


Axial Voracity (spn/U**' 
Mia Mu 


1AJ —chord l/U '-s pan 


Rjbyl Rxbyv2 Raby mx, Iuegrttod Convection d subccrc d jetcare Comments: 
Am rpn itim crd y* mm y ft mm 


-48.08 
-42.92 
-15.39 
-3ft 88 
35.84 
29.40 
3050 
-42.93 
-40.47 
39.79 
32.03 
35.21 
16.62 
-31.74 
-41.09 
-27.61 
-31.62 
-16.08 
-19.88 
-30.17 
18.56 
-17.98 
-25.89 
-19.00 
-22.42 
-13.00 

- 22.20 

-3185 

22.50 

-15.80 

-22.55 

-18.39 

-19.35 


6102 

75.48 

64.75 
65.67 
7128 
6190 

75.11 

51.58 
66.17 
7ft 10 
68.10 
7189 
63.37 

74.84 
70.82 
69.31 
73.87 
66.95 
69.66 

80.48 

60.58 
85.79 

70.11 
64.33 
67.27 
69.66 

111.20 

11100 

98.90 

57.70 

53.84 
120.14 

109.75 


z*-06 

r* ■ 06 

on 

205 

016 

1.97 

016 

200 

020 

1.91 

021 

1.89 

024 

210 

024 

204 

024 

210 

022 

1.87 

023 

1.99 

024 

1.95 

023 

1.80 

025 

202 

027 

206 

023 

1.76 

029 

213 

027 

205 

031 

1.92 

034 

1.82 

038 

1.76 


038 

1.07 

073 

1677.09 

31294.82 

0.090 

1.960 

0.231 

5.024 

037 

0.78 

067 



0.090 

2940 

0.222 

7.246 

136 

1.29 

077 

1693.71 

2107010 

0.090 

2940 

0.262 

8.544 

1.60 

1.18 

079 

1596.76 

14888.79 

0.090 

3.920 

024 

10.47 

136 

1.05 

053 



0.120 

5.358 

0.248 

11.1 

0.94 

039 

049 



0.150 

6.860 

0.268 

12.27 

136 

1.02 

1.02 



0.060 

2809 

025 

11.71 

1.27 

0.99 

073 



0.120 

5.619 

024 

11.24 

1.48 

1.16 

077 



0.090 

4.214 

022 

1031 

0.96 

0.86 

059 



0.120 

5.880 

022 

10.77 

139 

1.13 

075 



0.090 

4.606 

0.237 

12.11 

0.96 

038 

086 



0.090 

4.606 

0.24 9 

12.77 

1.06 

0.95 

042 



0.150 

7.840 

0.263 

13.72 

0.90 

031 

074 



0.090 

4.802 

0.234 

12.47 

1.18 

037 

056 

184631 

13781.78 

0.120 

6.534 

0.213 

1U8 

0.93 

0.60 

036 

1719.28 

1283838 

0180 

9.801 

0.245 

1331 

037 

031 

076 

1644.66 

10229.90 

0.090 

5.880 

0.205 

1339 

0.97 

0.82 

060 



0.120 

8.232 

0.228 

15.67 

1.01 

034 

040 



0.150 

10944 

0.228 

16.64 

233 

134 

064 

155338 

828135 

0120 

9.147 

0.191 

1438 

0.96 

030 

036 

1789.61 

9541.29 

0150 

11.434 

0.204 

1535 

0.91 

032 

074 



0.090 

7.350 

0.188 

1532 

0.96 

035 

061 

1494.77 

6973.18 

0120 

10454 

0.176 

1532 

037 

031 

056 

1438JJ5 

5963.20 

0.120 

11.761 

0.146 

1433 

131 

1.02 

1.01 



0.060 

7.187 

0.096 

11-52 

0.93 

033 

067 

2178.10 

27095.97 

0.090 

2940 

0.419 

13.68 

1.61 

1.04 

1.04 

2281.44 

21286.08 

0.060 

2613 

0.415 

18.08 

034 

0.80 

036 

224265 

16739.37 

0.180 

9.801 

0391 

2137 

_ 

-78.54 

1.08 



0.090 

5.390 

0327 

1937 

_ 


036 



0.240 

15.681 

0.18 

11.78 

032 

039 

089 


- 

0.060 

3.332 

0331 

1838 



035 



0.150 

8.004 

0336 

17.93 


windy 
breakdown 
breakdown 
zero axial vel 


23.60 

-24.18 

-15.60 

-15.90 

-11.83 
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055 
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034 

0.12 
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0.99 
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0.93 

0.83 
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167138 
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0.060 

4.438 
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Mu 
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Mu 
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Mu 
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20 
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2.18 

-0.82 

1.27 

1.05 

1.42 

066 

033 

-17.16 

98.83 

20 

050 151 l.Zdat 

1.16 

2.28 

-0.82 

1.07 

-1.00 

1.15 

067 

032 

-6.48 

1 1 1.75 

20 

050 161 l.dat 

1.00 

2.35 

-0.92 

1.39 

-1.09 

1.13 

065 

032 

-1064 

112.10 

20 

050 1711.dat 

092 

2.09 

0.83 

1.23 

-1.13 

1.18 

065 

032 

-15.85 

96.20 
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N/A 
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068 
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d aubcore 

d jetcare 
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y/» 

mm 

y* mm 

align errors 



N/A 


d aubcore 

d jetcare 
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9-63 
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Rel6" erd/oe er 
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638 
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ORIGINAL PAGE IS 

OF POOR QUALITY 



APPENDIX C 


ASSOCIATED ERRORS 


Overheat Error 


Using the error of the anemometer as ± 0005 or 0.002 ohms, an estimate of the error on a 
required overheat ratio can be determined. For the purposes here, an overheat of 1.8 will be used. 
Assume that the nulling of the probe to obtain the cold wire resistance has resulted in a value of 
4.000 ohms. This would then represent a range of 3.998 to 4.002 Q. For an overheat of 1.8 the 

resulting operating resistance would be set to 

(1.8 * 4.000 Q) = 7.200 Q. 


but the range of this value is 7. 198 Q to 7.202 Q. Thus the actual overheats could range from 


7.198 

4.002 


1.799 


7.202 
t0 3.998 


1.801 


or 


± 0.08% 


Velocity Measurement Error 

The nature of the equations used to determine the velocity field data makes error analysis 
difficult using standard methods, such as that outlined by Kline (1985). Thus, a comparison with 
data obtained by other previous investigators was used to qualify the present data. An indication of 
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the sensitivity of the output of the velocity equations to changes in the measured input voltages can. 
however, be indicated using a sensitivity analysis. 

Two spatial locations were considered. Both were examined using the data at the x/c = 0.5 
chord station for the 75° sweep wing at 20° angle of attack. The first spatial location was on the 
core axis, cotresponding to y/s = 0.69 and z/s = 0.39. Listed below are the measured voltages for 
the three wires at that location, along with the standard deviations and the resulting velocities using 

the equations (3.13 a-c). 


VI = 3.202 volts a = 0.87% 
V2 = 3.114 volts a = 0.77% 
V3 = 3.211 volts a = 1.6% 


u/U°° = 2.318 
v/U°° = 0.2917 
w/U°° = -0.4762 


The input voltages were then altered, one at a time, to examine the effect on the calculated 
velocities. Each voltage was first changed by ± its standard deviation. Then a 2% change was 
applied to all the cases. The resulting change in the velocities, compared to the baseline values are 

listed below. 


AVI 

AV2 

AV2 

+0.87% 

0% 

0% 

-0.87% 

0% 

0% 

0% 

+0.77% 

0% 

0% 

-0.77% 

0% 

0% 

0% 

+1.6% 

0% 

0% 

-1.6% 

+2.0% 

0% 

0% 

-2.0% 

0% 

0% 

0% 

+2.0% 

0% 

0% 

-2.0% 

0% 

0% 

0% 

+2.0% 

0% 

0% 

-2.0% 


Au % 

Av % 

Aw % 

+0.8 

-81.8 

+25.8 

-2.2 

+80.6 

-3.3 

+2.2 

- 7.4 

-63.7 

-2.5 

+18.3 

+39.2 

-3.1 

+99.4 

+45.5 

+0.3 

-13.4 

-4.9 

+0.4 

-53.0 

+70.7 

-5.6 

+ 157. 

+3.5 

+2.0 

+ 110. 

-63.9 

-6.3 

+43.0 

+84.0 

-6.1 

-160. 

+78.1 

-2.2 

+76.9 

+34.2 


It is immediately obvious that, although the axial velocity changed by no more than about 
3% for all the standard deviation perturbations, the other components demonstrated exceedingly 
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high differences when compared to the baseline values. The v velocity component is seen to be 
most sensitive to the change from V 1 . Less than a 1 % change can cause over 80% change in the 
value of v. The w component shows similarly large deviations for changes in V2. 


The above case consists of one large component, u/U,*,, and two smaller components, 

roughly an order of magnitude below the axial velocity. Therefore a second spatial location was 
chosen from the same survey plane where the velocities were of similar magnitude. The location 
selected was y/s = 0.51 and z/s = 0.24. The resulting voltages and velocities were: 

VI = 2.434 volts a = 3.5% u/U^ = 0.899 

V2 = 2.488 volts a = 2.3% v/U^ = 0.937 

V3 = 2.705 volts a= 1.6% w/U^ = -0.764 


Again the input voltages were perturbed by their standard deviations, which can be seen to 
be larger than the previous case, and the resulting change in the output velocity was recorded. 
Considerably more uniform variation in the output velocities is observed. The axial component 

shows larger changes than before, however none of the other components show changes greater 
than those of u/U^. 


AVI 

AV2 

AV3 

Au % 

Av % 

Aw % 

+3.5% 

0% 

0% 

+17.4 

-15.5 

-1.3 

-3.5% 

0% 

0% 

-18.8 

+ 12.2 

+ 1.9 

0% 

+2.3% 

0% 

+11.0 

+ 1.7 

-14.4 

0% 

-2.3% 

0% 

-12.2 

-0.7 

+12.6 

0% 

0% 

+1.6% 

-15.4 

+9.6 

+14.2 

0% 

0% 

-1.6% 

+12.3 

-9.6 

-14.9 


Vorticity Error 


An estimate of the error in the axial vorticity based on the error in the spatial location can be 
determined using the standard method of Kline (1985). The vorticity is defined as 
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Ox = 


8w dv 
" 5z 


or in terms of discrete values 


_ Aw A^ 
x — Ay Az 


(Cl) 


(C2) 


The uncertainty in Ox , Uq^, is determined from 


f 7" " 80x^2 Z" 8Qx S^ 77r dQx S^ , f TT dO x \2 
u Ox = y ^Ay 8 Ay J + ( U Az 8 Az J + ( U Aw 9 AwJ ( Av 9 Av J 


(C3) 


or 


U ftx = 


AV 

Av Az I 


(C4) 


Disregarding the velocity errors for the moment, the uncertainty based only on spatial error 
can be determined. Hence 


un, = -\/(v0) 2 + ( u a*i?y 


(C5) 


If the uncertainty in the spatial location is set equal to the resolution of the system, 

U Ay = 0.0254 mm and = 0.0208 mm 

The calculated uncertainty can be expressed as a fractional uncertainty of the value of O x itself by 
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eliminating the need to calculate an actual value of Q x . 


Three stations were selected along the 75° sweep planform, x/c = 0.3, 0.5, and 0.8. The 
values of Ay/s = Az/s = 0.03 for the data acquired, however in order to maintain the same spatial 
location for the vorticity values the spatial distance used in equation (C2) is twice this. The actual 
values of Ay and Az depend on the chordwise station at which they are determined, as shown 
below, along with the calculated fractional uncertainties for these using equation (C6). 


x/c 

Av=Az (ram) 

2Av (mm) 

Uq^/Qx % 

U n x /Q x % (lmm) 

0.3 

0.98 

1.96 

0.84 

36.1 

0.5 

1.63 

3.26 

0.50 

21.7 

0.8 

2.61 

5.22 

0.31 

13.5 


Payne (1987) states that the overall positional accuracy of the system is ± 1 mm. If this is 
interpreted as the uncertainty in Ay and Az, the fractional uncertainty of the vorticity increases 
dramatically to the percentages noted in the last column. 

To include the effect of the velocity error, equation (C4) must be employed. An error 

estimate for the velocity is required as well as some values for Aw and Av. Using the values of Ay 
= Az = 3.26 mm at the x/c = 0.5 station, and selecting values of Aw = -Av = O.IU^ m/s, a value 

of Q x = 61.3 1/s is determined. If the error in the velocities is only taken to be 15% of the 

values obtained in the core, the uncertainty from (C4) is calculated to be over 80% of the value of 

Qx- 

Circulation Error 

The error in the calculated circulation can be calculated in a similar manner as the vorticity. 


The circulation, defined as 
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r = I V • dr (C7) 

r 

can be expressed as a summation of terms, each of which has an error associated with it. 

r = £(VAy) n (C8) 

i=l 


where Ay = Az. Thus, ignoring the error due to velocity, the uncertainty based on the error in the 
spatial location can be derived as 


U r = 



ar 

Ayl 3Ayi 



ar Y 

Ay 2 3Ay2 J 


+ 


(C9) 


where 


ar = ar 

aAyi ~ 3Ay n 


If the values of V n are assumed to be the same along a path of integration, or for the worst case 
scenario V n is set equal to the largest value along the: 


u r 


= V n V 


U 


Ay 1 


+ U 


Ay2 


or ‘ 



(CIO) 
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As an example, the station of x/c = 0.5 is used again. If the integration path is taken to be out to a 
radius of r/s = 0.25, the value of V n can be set to the maximum of V n = 1.3 from Figure 
5.15b. The circulation at this station and radius was calculated to be approximately 0.2* (cU^) 
from Figure 5.29. At this radius approximately 34 summations are required based on the local grid 
spacing of y/s = z/s = 0.03. Using the uncertainty in the grid spacing of 0.0254 mm for both y and 
z directions yields a value of 


Up = 0.193 e-03 U^ 


Up 

— = 0.24% 

r 


about half that of the vorticity for the same case. If the maximum uncertainty in Ay of 1 mm is 
used, however, an error of 9.3% is determined, again about half that of the vorticity. This, of 
course, assumes a negligible error in the velocity. 



APPENDIX D 


FURTHER CROSS WIRE DERIVATIONS 

The technique of Sherif and Pletcher (1987) uses two positions of an X-wire array with the 
wires are in a plane perpendicular to the probe axis. Use of this method gave success in terms of 
magnitudes that were comparable to other measurements of the same delta wing flowfield. These 
equations were rewritten for the probe configuration of the two wires parallel to the oncoming 
stream and at ±45° and are detailed below. The resulting equations were similar to that obtained 
by Sherif and Pletcher. They also consisted of a set of three nonlinear equations with three 
unknowns plus a fourth equation. The direction of V and W could also be determined from the 
voltages, as was shown earlier. Hence if these nonlinear equations could be solved, the velocity 
vector would be known. This derivation is now explained in more detail. 


Recall the general geomertry for a single hot wire introduced in Chapter 3 and based on the 
figure below 



The effective velocity measured by the wire was expressed by Jprgensen (1971) as. 
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U e ff^ “ U N^ + k i^U T ^ + k 2 ^U B 2 

As outlined in Chapter 3, the present study required four spatial passes to uniquely determine the 
unknown flowfield. It was desirable to minimize the number of grid sweeps required to reduce 
both the data acquisition time and the potential for error. Several schemes were tried, including 
iterative techniques based on simplified forms of the Jorgensen equation, however none were 
successful. The original derivation of Sherif and Pletcher (1987) was then applied to the 
configuration with the wires lying in a plane parallel to the probe axis, using the geometry below. 



Figure D2 Single Slant Wire Geometry 


The velocity is decomposed into the components below 
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U N = u sin0 + v cos0 cos0 + w sin0 cos0 
U T = - u cos0 + v cos0 sin0 - w sin0 sin0 
U B = v sin0 + w cos0 

where 0 is the angle of the wire relative to the x-axis and 0 is the rotation angle of the probe, 0 
being equal to zero degrees in the x-y plane. 

For an x-wire in the x-y plane, the two wires lie at 0 = 45° and 0 = 135° with 0= 0°. This is 
also equivalent to 0 = 45° and 0 = 0° and 1 80°. Hence 

a) 0 = 0°, 0 = 45°: U N 

U T 

U B 

Therefore: \ij = + k ,2 + ^2 *2 

Similarly 

b) 0 =0°, 0 = 135“: U efr 2 = + k) 2 + ^2 w 2 

2 2 

Hence: a) U eff 2 = ^-(l+kj 2 ) + y (l+k, 2 ) + k 2 2 w 2 -MivO-kj 2 ) 

b) U eff 2 = ^-(1+k! 2 ) + y (1+kj 2 ) + k 2 2 w 2 - uv(l- kj 2 ) 


= U a 2 (Dl) 

= U 5 2 (D2) 


U V 

f2 f2 

V u 

V2 ' V2 


f <U + V) 
f(v-u) 


= w 


Rotation of the probe to 0 = 90° yields 
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0 = 45° 


t j 2 _ (u+w) z 2 (-u-w T + k 2 V 2 

°) u eff ” 2 + K 1 2 k 2 


2 2 

U e ff 2 = yO+kj 2 ) + k 2 2 v 2 + (l+k^ + uwCl+k^) = U ( 


2 (D3) 


0=135° 


2 ( u * w > , k 2 ( u ~ w ) + k 2 V 2 

d > u eff = ~^r + k l 2 + k 2 v 


U eff 2 = yO+k! 2 )+ k 2 2 v 2 + (1+kj 2 ) - uw(l+ kj 2 ) = U d 2 (D4) 


Thus, four nonlinear equations in three unknowns result with the effective velocities U a , U b , 
U c , and Ud determined from the calibration curves. The directions for v and w are simply found 
in the same manner as the present method. The angles are defined from the x-axis as 


a = 45 - tan* * 




u 

V a 2 


angle of the velocity component in the xy plane 


Y = 


45 - tan" 1 


l U cJ 


angle of the velocity component in the xz plane 


Further manipulation of equations 1 to 4 yields 

U a 2 - U b 2 = 2uv(l-k 1 2 ) 

U c 2 - Uj 2 = 2uw(l+kj 2 ) 

U a 2 + U b 2 = u 2 (l+ kj 2 ) + v 2 (l+ k^ 2 ) + 2 k 2 2 w 2 


(D5) 

(D6) 

(D7) 
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U c 2 + U d 2 = u 2 (l+ k x 2 ) + 2k 2 2 v 2 + w(l+k! 2 ) 

Subtracting (D8) from (D7) and defining the result as £ 

£ = U a 2 + U b 2 - U c 2 - U d 2 = v 2 (1+ k^ 2 - 2k 2 2 ) + w 2 (2k 2 2 - 1 - kj 

Dividing equation (D5) by (D6) and denoting this as (3. 

U a 2 ' u b 2 = v < 1 - k l 2 > - B 
u c 2 - U d 2 wfl+kj 2 ) 


Rearranging 


v = w (3 


(l+kj 2 ) 

(1- k! 2 ) 


or 


■> •> „•> < i+ k i 2 > 2 

' - w w 


= w 2 o 2 


(1+k! 2 ) 


where a equals (5 multiplied by — ^2) 


Finally, substitution of (DIO) into (D9) 

£ = w 2 a 2 (1+ k t 2 - 2k 2 2 ) + w 2 (2k 2 2 - 1 - kj 2 ) 
or 

£ = w 2 (k 1 2 (o 2 - 1) + 2 k 2 2 (l - a 2 ) + a 2 - 1) 

and solving for w 


(D8) 


) (D9) 


(DIO) 


(Dll) 
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w 


ki 2 (c 2 
V 1 


c 

1) + 2 k 2 2 (l - ct 2 ) + a 2 - 


\0.5 

1 

J 


CD 12) 


Thus it would appear that w, and susequently u and v, could be solved explicitly. 


Unfortunately, lines of very large or small magnitudes were observed in the velocity fields 
calculated from equation (D12) and the resulting expressions for u and v. In addition, the values 
of u determined were not representitive at all topographically, although the magnitudes were 
comparable. From equation (D12) it can be seen that when a approached 1, the value of w 
becomes infinite. Equation (DIO) indicates that for a = 1, w = v, and thus it may be that a separate 

solution must be determined when w = v. To determine this, v was set equal to w in equations 
(Dl) to (D4) 

2 2 

U a 2 = ^-(1+kj 2 ) + ^-(l+kj 2 )+ k 2 2 w 2 +uw(l-k! 2 ) 

2 2 

Ujj 2 = ^-(1+kj 2 ) + ^-(l+kj 2 )+ k 2 2 w 2 - uw(l-kj 2 ) 

U c 2 = ^-<l+k 1 2 )+ k 2 2 w 2 + (1+kj 2 ) + uw(l+ kj 2 ) 

U d 2 = y(l+kj 2 )+ k 2 2 w 2 + ^-(1+kj 2 ) - uw(l+ kj 2 ) 

Now 

U a 2+ u b 2= u 2 (l+ k x 2 ) + w 2 (l+k 1 2 +2k 2 2 ) (D13) 

and 

u a 2 ' u b 2= 2uw(l-k 1 2 ) 
or 

»a 2 - Ub! 

2w(l- kj 2 ) 


u 


(D14) 
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Substituion of (D14) into (D13): 


U a 2 * U b 2 = 


('U a 2 -U b 2 '| 


2w(l- 


(l+k 1 2 )+ w 2 (l+k 1 2 +2k 2 i ) 


or 


w 


w 


2 U a 2 ^ u b 2 




1+ ki 


1+ kj 2 + 2k 2 2 + 4(1- kj 2 ) 2 l+k 1 2 +2k 2 i 


= 0 


(D15) 


Solutions to equation (D 15) also failed to produce velocity fields without regions of non- 
convergence. Newtons iteration scheme helped, however constraints (as in P ± .2) were required 
to get rid of the singularities. Thus, for the present tests, the method outlined in Chapter 3 was 


employed. 
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